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ABSTRACT 


Larralde,  Jesus.  Ph.D.,  Purdue  University,  December  1984. 
Structural  Analysis  of  Rigid  Pavements  with  Pumping.  Major 
Professor:  Dr.  Wai  Fah  Chen. 


Pumping  has  been  considered  a  major  contributor  to 
premature  failure  in  rigid  pavements.  The  traditional 
design  procedures,  based  on  a  stress-fatigue  criterion, 
have  been  unsuccessful  in  providing  pavements  which  accom- 
plish their  intended  service  lives.  In  addition  to  the 
stress-fatigue  criteria,  an  erosion  criteria  must  be  con- 
sidered. Concrete  pavements  have  been  analyzed  using  the 
Westergaard's  equations,  or  charts  developed  from  those 
equations.  These  methods  assume  full  contact  between  the 
slab  and  subgrade,  as  well  as  continuity  in  the  whole 
pavement  structure.  Recently,  with  finite  element  tech- 
niques, the  consideration  of  discontinuities  in  the  pave- 
ment structure  has  been  possible.  Nevertheless,  in  the 
traditional  analysis  methods,  closed  form  solutions  or 
finite  element  methods,  the  reduction  of  pavement  strength 
resulting  from  the  damage  caused  by  traffic  has  not  been 
considered.  Fatigue  damage  caused  by  repetitive  load 
applications  reduces  the  strength  and  stiffness  of  the 
pavement . 


XV 

In  this  thesis,  a  nonlinear  finite  element  method  for 
analysis  of  rigid  pavements  which  considers  the  decay  of 
the  their  stiffness  properties  due  to  fatigue  damage  is 
presented.  The  damage  caused  in  the  concrete  slab,  the 
load  transfer  devices,  and  the  pavement  foundation  are 
considered.  The  computer  scheme  has  been  developed  from 
the  previously  existing  ILLISLAB  code.  The  computer 
implementation  presented  herein,  called  PMARP,  has  the  fol- 
lowing main  features:  nonlinear  finite  element  method; 
orthotropic  plate;  consideration  of  reinforcement  in  con- 
crete slab;  resilient  modulus  of  subgrade  reaction;  pump- 
ing development;  determination  of  states  of  stresses  and 
deformations  at  all  nodal  points;  and  determination  of 
state  of  pavement  damage.  This  method  can  be  used  for 
design  and  for  prediction  of  maintenance  requirements  of 
rigid  pavements . 


CHAPTER  1   INTRODUCTION 


1.1  Int  roduc t  ion 


Historically,  pavements  have  been  divided  into  two  general 
types  [l]i  flexible  pavements,  which  generally  consist  of 
a  relatively  thin  wearing  surface  built  over  a  base  course 
and  subbase  course,  which  in  turn  rest  upon  the  compacted 
subgrade;  and  rigid  pavements  made  up  of  a  Portland  cement 
concrete  slab  and  sublayers  termed  bases  or  subbases 
between  the  slab  and  subgrade  (Fig.  1.1).  The  research 
reported  in  this  thesis   deals  with  rigid  pavements  only. 


SLAB 


SUBLAYERS 


Figure  1.1  Cross  section  of  rigid  pavement. 


_1_.2_  Statement  of  the  Problem 
Since  the  1940's,  pumping  has  been  considered  as  a  major 
contributor  to  failure  in  concrete  pavements.  It  has  been 
reported  [2,3,4,5]  that  most  of  the  failures  in  concrete 
pavements  have  been  preceded  by  pumping  of  the  material 
underneath  the  concrete  slab.  The  traditional  design  pro- 
cedures, based  on  a  stress-fatigue  criterion,  have  been 
unsuccessful  in  providing  pavements  which  accomplish  their 
intended  service  lives.  Premature  failure  due  to  excessive 
pumping,  erosion  of  sublayers  and  joint  faulting  have 
occurred.  Therefore,  in  addition  to  the  stress-  fatigue 
criteria,   erosion  criteria  must  be  considered. 

As  part  of  the  design  process,  the  structural  analysis  of 
rigid  pavements  has  been  used  to  determine  the  stresses 
produced  by  certain  standard  loads  at  specified  critical 
points  of  the  pavement.  Thickness  design  has  been  based  on 
limiting  these  critical  stresses  at  appropriate  levels  to 
control  fatigue  damage  of  concrete. 


Generally,  concrete  pavements  have  been  analyzed  using  the 
Wes ter gaard 's  equations  [6]  or  charts  developed  from  these 
equations  [7]  .  These  methods  assume  full  contact  between 
the  slab  and  sublayer,  as  well  as  continuity  in  the  slab. 
More  recently,  the  development  of  finite  element  analysis 
has   allowed   the  consideration  of  discontinuities  such  as 


joints,  cracks,  and  loss  of  sublayer  support.  Addition- 
ally, variable  wheel  placement,  variable  material  proper- 
ties, and  variable  load  transfer  at  joints  can  be  con- 
sidered with  this  method. 

Nevertheless,  the  reduction  of  strength  due  to  the  com- 
bined effects  of  load  repetition  and  erosion  has  not  been 
considered.  Packard  and  Tayabji  [8]  have  indicated  that  in 
addition  to  the  failure  of  pavements  caused  by  excessive 
fatigue  cracking  of  the  concrete  slab,  pavements  can  also 
fail  due  to  excessive  pumping  and  erosion  of  the  sub- 
layers . 

The  conventional  methods  for  structural  analysis  of  con- 
crete pavements  do  not  consider  the  variation  in  the  pave- 
ment properties  caused  by  repetitive  application  of  loads. 
These  methods  are  conservative  in  the  sense  that  they 
assume  the  geometric  and  mechanical  properties  to  remain 
at  the  same  values  through  the  life  of  the  pavement. 

The  major  limitation  in  the  analysis  method  developed  by 
Westergaard,  or  those  based  on  his  equations,  is  that  it 
is  applicable  only  if  the  analyzed  structure  is  continu- 
ous. A  rigid  pavement  is  composed  of  several  elements  with 
different  mechanical  behaviors.  In  general,  the  pavement 
can  be  divided  into  a  concrete  slab,  sublayers  and  load 
transfer  devices.  Each  one  of  these  components  behaves  in 
a  completely  different  manner,  therefore,  it  is  impossible 


to  use  a  general  differential  equation  to  represent  the 
behavior  of  the  whole  structure.  Additionally,  the 
existence  of  joints,  cracks,  and  gaps  between  slab  and 
sublayers  is  very  common  which  makes  the  analysis  more 
complex  and  the  solution  in  a  closed  form  impossible. 
Therefore,  these  methods  are  used  only  to  check  results 
obtained  with  other  methods  using  simple,  idealized  cases. 

One  of   the   main   causes   of   pavement   failures   is   the 
deterioration   and   gradual  weakening  of  the  pavement  com- 
ponents.  Repetitive  loading  produces  micr o era  eking  in  the 
slab   as  well  as  erosion  and  voids  in  the  sublayers.   This 
reduces  the   resistance   of   the   pavement   to   the   loads 
applied   by   traffic.    The   response  of  the  pavement  to  a 
specific  load  is  different  throughout  its   entire   service 
life.    Cumulative  damage  reduces  the  general  stiffness  of 
the  structure.   Therefore,  in  the  analysis  of  rigid   pave- 
ments, especially  if  erosion  and  pumping  are  expected,  the 
pavement  components   should   be   considered   as   having   a 
stress   dependent   behavior.    The  response  of  the  slab  as 
well  as  that  of  the  sublayers  should  be  dependent  upon  the 
previous   stress  history  applied  on  the  pavement.  Pavement 
properties  should  be  updated  every  few  years  to  take   into 
account   the  reduction  of  its  stiffness  caused  by  exposure 
to  traffic  loads.  In  this  way,  a  more  exact  calculation  of 
stresses   and   deflections  at  different  stages  of  the  life 
of  the  pavement  is  accomplished.  Therefore,  a  more  precise 


estimation   of   the   number   of  allowable  load  repetitions 

(and  thereby  of  service  life)  can  be  achieved.  Then,   the 

designer  has  more  and  better  information   for  an  econom- 
ical ana  lysis . 

1_.3_  Objective 
The  main  objective  of  the  research  presented  here  is  to 
develop  a  method  of  structural  analysis  of  concrete  pave- 
ments which  considers  the  discontinuities  in  the  pavement 
structure,  as  well  as  the  long-term  effects  produced  by 
traffic.  The  research  presented  herein  is  part  of  pro- 
ject DT  FH61-82-C-00035  ,  sponsored  by  the  Federal  Highway 
Administration,  to  develop  a  method  of  design  to  prevent 
pumping  in  rigid  pavements.  This  project  is  divided  into 
four  main  tasks:  (1)  development  of  a  simplified  labora- 
tory method  for  evaluating  pump  susceptibility  of  materi- 
als; (2)  development  of  methods  to  detect  voids  under  con- 
crete pavements;  (3)  development  of  a  theoretical  pro- 
cedure for  structural  analysis  of  concrete  pavements;  and 
(4)  selection  of  methods  of  economic  trade-off  analyses. 
The  research  presented  herein  deals  with  the  development 
of  theoretical  procedures  for  structural  analysis  of  rigid 
pavements. 


To  accomplish  the  objective,  a  non-linear  finite  element 
method  for  analysis  of  rigid  pavements  has  been 
developed.  This  method  of  analysis  takes  into  account   the 


damage  produced  in  the  concrete  slab,  in  the  load  transfer 
devices,  and  in  the  sublayers  by  the  repetitive  applica- 
tion of  traffic  loads.  The  computer  implementation  of  the 
method  is  based  on  the  computer  method  ILLISLAB  previously 
developed  at  the  University  of  Illinois  [9],  In  the 
method  presented  herein,  the  pavement  components  are 
assumed  to  have  a  stress-  or  strain-dependent  behavior, 
and  the  concrete  slab  is  considered  as  an  orthotropic 
plate . 


_1_.A_  Thes  is  Arrangement 
The  present  thesis  is  divided  into  two  parts.  In  the 
first  part,  consisting  of  Chapters  2,3,  and  4,  the  avail- 
able methods  for  analysis  of  concrete  pavements  are  anal- 
ized.  Chapter  2  is  a  literature  review  of  the  methods  of 
structural  analysis  in  which  closed  form  solutions, 
influence  charts,  and  finite  element  methods  are 
presented.  In  Chapter  3  a  theoretical  evaluation  of  these 
methods  is  presented,  while  Chapter  4  deals  with  a  numeri- 
cal evaluation  of  them.  In  the  second  part,  which  consists 
of  Chapters  5,  6,  and  7,  a  finite  element  method  for 
analysis  of  rigid  pavements  with  fatigue  in  the  concrete 
slab,  fatigue  in  the  load  transfer  devices,  and  pumping  is 
presented.  In  Chapter  5,  the  method  of  analysis  and  the 
computer  Implementation  are  presented.  Numerical  results 
obtained  using  the  proposed  method  are  also  given. 


CHAPTER  2   LITERATURE  REVIEW 


2.1  Introduction 


For  over  50  years,  efforts  have  been  made  to  develop  reli- 
able analysis  techniques  to  improve  the  design  and  perfor- 
mance of  concrete  pavements.  First  attempts  to  analyze 
pavements  were  made  by  Westergaard  in  the  early  1920's 
[10].  Later  Pickett  and  Ray  in  1951  [7]  used  the  work, 
done  by  Westergaard  to  develop  influence  charts  which  were 
used  by  The  Portland  Cement  Association  for  the  design  of 
highway  and  airport  pavements.  Since  the  1960's  discrete 
elements  and  finite  element  methods  have  been  used.  In 
this  chapter,  the  development  of  the  different  methods  for 
structural  analysis  of  rigid  pavements  is  illustrated. 
Also,  a  brief  description  of  the  most  common  methods  for 
design  of  rigid  pavements  is  presented. 

2.1.1  Closed  Form  Solutions. 
Closed  form  solutions  of    slabs    on   elastic   sublayerss 
were   obtained   by   Westergaard   [10,11,12,13],  Hogg  [14], 
Volterra  [15],  Bergstrom  [16],  and  Tiraoshenko  [17]. 


These  solutions  were  based  on   the   differential   equation 
developed   by   Lagrange  for  elastic  isotropic  thin  plates. 


This  equation  is  as  follows 
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Where       q       is     the    distributed       lateral       load,         D         is       the 
flexural    rigidity    given    by: 


D    = 


Et 
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(2.2) 


w       is     the    deflection    at    the    slab    with    elastic       properties 

E       and       \i    . 

For  a  plate  resting  on  an  elastic  sublayers,  the  differen- 
tial equation  becomes: 


3  w 


„4  *A 

?   3  w  3  w    q-kw 

4       22      4     D 

3x      3x  3y  3y 


(2.3) 


Where  k  is  the  so-called  modulus  of  subgrade  reaction. 
Note  that  the  term  "subgrade"  is  used  regardless  of  the 
proper  designation  of  the  sublayer  in  the  pavement  sec- 
tion. 

Even  though  Westergaard  did  not  provide  details  of  the 
way  this  equation  was  solved,  an  approximate  expression 
for  the  displacement  in  the  form  of  a  exponential  series 
was  assumed.  Using  the  first  two  terms  of  these  series, 
the  following  equation  was  obtained 


w  =  ZJL  (l.le'ex  -  O.88c.e'20X) 

k  1 


(2.4) 


1  Ik 
Where,   3  =   (k/D)         ,   and   P   is   the   total   load. 

Westergaard   obtained   expressions   for  stresses  caused  by 

three  loading  cases:  1)  load  applied  near  the  corner  of   a 

large  rectangular  slab;  2)  load  applied  near  the  edge  at  a 

considerable  distance   from   the   corner   ;   and   3)   load 

applied  at  the  interior  of  a  large  slab.  These  expressions 

are  as  follows: 
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where  a  is  the  radius  of  area  of  load  contact  (in.). 
This  area  is  considered  circular,  in  the  cases  of  corner 
and  interior  loadings,  and  semicircular  in  the  case  of 
edge     load.     The    value    of     b     is    given    by: 


2       2 

b    =    \|l.2a    +t         -0.675t     .    The    quantity       1  ,       radius       of 

relative    stiffness,     is    given    by: 


1    =    (         Et    2       )1/4 
12(l-u    )k 


(2.8) 


The  above  equations  were   developed   using   the   following 
as  sump t  i  ons  : 


The  slab  acts  as  a  homogeneous   isotropic   linear   elastic 
solid  with  full  support  from  the  sublayer. 
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The  reaction  of  the  sublayer  at  a  given  point  is  vertical 
and  proportional  to  the  slab  deflection  at  that  specific 
point.  This  is  called  a  Winkler  foundation,  in  which  the 
constant  of  proportionality  is  called  the  modulus  of 
subgrade  reaction  k  .  The  load  is  distributed  uniformly 
over  a  circular  or  semicircular  area. 

The  Winkler  foundation  assumption  was  modified  in  1938  by 
Hogg  [14],  and  Holl  [18].  They  considered  the  sublayer 
acting  as  an  elastic  solid  instead  of  the  Winkler  founda- 
tion. These  two  researchers  solved  the  problem  of  a  thin 
slab  of  infinite  size  supported  by  a  semi-infinite  elastic 
solid.  However,  they  only  considered  symmetrical  loading. 
Pickett  et  al  [19]  later  obtained  influence  charts  to 
determine  stresses  and  deflections  for  any  distribution  of 
load.  These  charts  are  given  in  reference  [19], 


Even  though  this  method  can  be  used  for  verification  pur- 
poses in  simple  problems,  the  major  limitation  is  that 
discontinuities  and  variable  properties  cannot  be  con- 
sidered. It  assumes  full  support  conditions;  therefore, 
loss  of  support  or  voids  cannot  be  modeled.  Variations  in 
the  thickness  of  the  slab  or  in  its  elastic  properties 
cannot  be  modeled  either.  Due  to  these  limitations  and 
with  the  development  of  computers,  other  techniques  based 
on  numerical  methods  began  to  be  studied.  Among  these, 
the  finite  element  method  has  had  great  acceptance. 
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2.1.2  Analysis  of  Joints 
Joints  in  rigid  pavements  are  necessary  to  facilitate  con- 
struction and  to  control  the  formation  of  cracks.  The  need 
for  joints  has  been  recognized  since  the  early  1900's. 
Initially,  joints  were  simple  openings  between  slabs,  load 
transfer  devices  appeared  around  1917,  and  from  1912  to 
1936  different  types  of  joints  with  and  without  load 
transfer  devices  were  developed  and  even  patented. 


One  of  the  first  attempts  to  analyze  joints  was  carried 
out  by  Timoshenko  and  Lessels  in  1925  [20].  Considering 
the  case  of  a  semi-infinite  beam  supported  on  an  elastic 
sublayers,  they  obtained  mathematical  expressions  for 
loads  and  deflections  for  a  bar  encased  in  concrete. 
Teller  and  Sutherland  [21],  Friberg  [22],  Sutherland  and 
Cashell  [23],  Kushing  and  Fremont  [24],  Giggin  [25],  Fin- 
ney and  Fremont  [26],  Marcus  [27],  and  Teller  and  Cashell 
[28],  are  among  the  researchers  concerned  with  the 
behavior  of  joints.  Teller  and  Sutherland  [21]  studied 
the  effects  caused  by  variations  in  joint  width,  dowel 
looseness,  and  dowel  stiffness  on  the  performance  of 
joints.  The  importance  of  placing  the  dowel  bars 
correctly,  that  is,  perpendicular  to  the  cross  section  of 
the  pavement,  was  studied  by  Smith  and  Benham  [29],  They 
concluded  that  the  bars  may  be  distorted  or  the  concrete 
damaged  in  the  vicinity  of  the  bars  if  they  are  not, 
within   certain   limits,  perpendicular  to  the  joints. 
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Friberg  [22] s tudied  the  bearing  stress  caused  by  dowel  bars, 
using  the  equations  developed  by  Timoshenko  and  Lessels 
[20]  and  a  Winkler  foundation.  With  that  model,  stresses 
produced  by  the  bar  on  the  surrounding  concrete  of  the 
order  of  68.9MPa  (10  000  psi)  were  obtained.  Later, 
Marcus  [27],  using  the  same  analysis  and  considering  a 
three  dimensional  state  of  stresses,  obtained  lower 
stresses.  He  stated  that  the  primary  cause  of  radial 
cracks  is  the  principal  tensile  stress  resulting  from  the 
tensile  and  shear  stresses  present  at  the  bar-concrete 
interface . 

The  effect  of  repetitive  loading  on  dowel  joints  was  stu- 
died by  Teller  and  Cashell  [28].  They  indicated  that  per- 
formance under  single  loading  is  no  measure  of  performance 
under  repeated  loading.  With  a  test  in  which  an  alter- 
nating load  is  applied  to  either  side  of  the  joint  of  a 
slab  specimen,  the  effects  of  dowel  diameter,  dowel  embed- 
ment, and  width  of  joint  opening  were  studied.  It  was 
indicated  that  in  order  to  transfer  certain  percentages  of 
the  load  across  the  joint,  the  dowel  diameter  should  be 
specified  proportionally  to  the  slab  thickness.  Repeated 
loading  developed  dowel  looseness  which  has  an  important 
effect  on  load  transfer  efficiency. 


A  clear  understanding  of  the  behaviour  of  joints,  as  well 
as  of  the  whole  pavement  structure,  is  important  in  the 
structural  analysis.   A  mechanistic  approach,  in  which  the 


13 

pavement  structure  is  divided  into  pieces  which  then  are 
modeled  by  appropriate  elements,  appears  to  be  the  way  to 
proceed.  Finite  element  methods  allow  such  an  approach. 

2.1.3  Finite  Element  Methods. 
The  structural  analysis  of  rigid  pavements  is  a  complex 
problem.  Before  the  development  of  the  numerical  methods 
and  the  availability  of  high  speed  computers,  pavements 
were  analyzed  by  the  methods  developed  by  Westergaard  or 
derivations  of  such  methods.  These  so  called  "closed  form 
solutions"  are  subject  to  severely  limiting  assumptions. 
With  these  assumptions,  pavements  are  considered  as  con- 
tinuous with  infinite  dimensions,  and  with  full  contact 
between  slab  and  subbase.  Therefore,  the  analysis  is 
unrealistic  since  it  is  not  possible  to  analyze  pavements 
with  cracks,  joints,  or  variable  support  conditions. 


The  limitations  imposed  by  these  solutions  prompted  the 
development  of  different  analysis  techniques.  Rather  than 
representing  pavements  as  idealized  structures  they  should 
be  represented  in  a  mechanistic  way.  Newmark  [30]  proposed 
a  numerical  technique  for  analyzing  elastic  bodies. 
Clough  [31]  and  Turner  [32]  developed  a  finite  element 
approach  to  the  analysis  of  complex  structures,  including 
plane  stress  analysis.  Finite  element  analysis  of  plates 
and  slabs  resting  on  elastic  sublayerss  were  developed  by 
Zienkiewicz  and  Cheung  [33],  Hudson  and  Matlock  [34], 
Huang  and  Wang  [35],  and  Huang  [36,37], 
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In  the  method  developed  by  Zienkiewicz  and  Cheung,  the 
concrete  slab  is  divided  into  rectangular  finite  elements 
treated  as  elastic  isotropic  plates.  These  elements  are 
connected  to  each  other  at  their  corners.  Each  element  has 
three  degrees  of  freedom  at  each  corner;  a  lateral  deflec- 
tion in  the  Z  direction,  and  two  rotations  about  perpen- 
dicular axes  X  and  Y.  These  12  quantities  constitute  the 
nodal  displacements  for  the  given  element  ;  u.  Similarly, 
at  every  node  there  are  three  forces;  two  concentrated 
moments  and  a  lateral  force;  F  .  The  nodal  forces  and 
displacements  can  be  related  as: 


(2.9) 


The  matrix  k  ,  composed  of  12x12  elements,  represents 
the  stiffness  matrix  for  the  element  e  .  Putting 
together  the  stiffness  matrices  of  all  the  elements, 
equilibrium  equations  at  the  joints  can  be  written,  then: 
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e  =  l 


(2.13) 


e-1,2,3,... 


Where  f.  f  ~  f-  ...  are  all  external  forces;  a.  a?  ...   are 
all  nodal  displacements;  And  n  is  the  number  of  elements. 


If  in  the  set  of  simultaneous  equations,  <  f  /=  [K] 

the  boundary  conditions  in  terms  of  forces  and  displace- 
ments are  included,  the  solution  of  the  problem  reduces 
to  the  solution  of   such  a  system  of  equations. 

g 
To  determine   k    ,  Zienkiewicz  and  Cheung  used  the   prin- 
ciple  of     virtual  work,   where  for  a   virtual  displace- 
ment the  external  work  must  be  equal  to  the  internal  work. 


ext 


(2.14) 


Where    6a        is  a  virtual  displacement   which   can   be 
taken  as  I,  and 
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To  compute  the  internal  work,   the   following   expressions 
developed  by  Timoshenko  are  used: 
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M  =  [D]  x     x  =  vector-of -curvature s 
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Then: 


T 
"in.  =  ^{6x}  M  dxdy 


(2.20) 


Whe  re 
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[B][C]  relates    the    curvature    to       the      nodal       displace- 

ments,    therefore: 


=    //  |tB] [C]"1}T[D 


] [B] [C]_1ue    dxdy 


(2.21) 


17 


[,.,"]•[ 


//    tB]T[D][B]     dxdy       [C]_1<a 


Then: 


[-1T[ 


//    [B]T[D][B]    dxdy 


-1    lei"1 


(2.22) 


This  matrix  for  a  rectangular  element  is  given  in  the 
paper  by  Zienkiewicz  and  Cheung  [33], 

Nodal  forces  are  obtained  by  assigning  forces  equivalent 
to  the  lateral  distributed  or  concentrated  loads.  This 
assignment  can  be  done  by  simple  statics,  or  in  such  a 
way  that  during  any  virtual  displacement,  the  work  done  by 
the  actual  forces  and  that  done  by  the  equivalent  concen- 
trated forces  is  the  same. 

To  analyze  slabs  on  elastic  foundations,  an  additional 
element  was  introduced  to  represent  the  support.  In  a 
paper  published  in  1965,  Cheung  and  Zienkiewicz  [38] 
showed  that  the  stiffness  matrix  coefficients  of  the  foun- 
dation are  simply  added  to  those  of  the  slab  element.  They 
obtained  expressions  for  these  coefficients  considering 
the  sublayers  as  a  Winkler  type  or  as  an  isotropic  elastic 
solid.  For  the  Winkler  case,  nodal  forces  and  displace- 
ments are  related  as: 


{,}   -  a,  k[„]   {w} 


(2.23) 
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In  which  a  and  b  are  the  sides  of  the  element,  k  is 
the  modulus  of  subgrade  reaction,  w  represents  the 
deflections  at  the  nodes,  and  a  values  are  coefficients 
to  account  for  the  area  contributing  to  the  nodal  forces. 


For  the  isotropic  elastic  solid,  Cheung  and  Zienkiewicz 
used  the  Boussinesq  equation,  arriving  at  the  following 
expression: 
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where  v    and  E  are  the  elastic  properties  of  the   founda- 


tion,  and    * 


is  the  flexibility  matrix  of  the  foun- 


dation which  is  obtained  from: 
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values    are  coefficients  which  depend  on  — . 
i  i  a 


The  sublayers  matrix  and  the  plate  matrix  give  a  set  of 
simultaneous  equations  for  the  whole  pavement  system.  The 
solution  of  these  equations  provides  the  nodal  displace- 
ments and  with  these  moments,  stresses  and  reaction  pres- 
sure can  be  determined. 

In  the  method  developed  by  Hudson  and  Matlock  [34],  in 
1966,  the  pavement  slab  is  defined  by  a  finite-element 
model  consisting  of  bars,  springs,  elastic  blocks  and  tor- 
sion  bars.   Combining   beams  in  the  X  and  Y  directions,  a 
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grid-beam  system  is  formed  which,  coupled  with  torsional 
bars,  represent  the  concrete  slab.  The  beams  are  con- 
nected with  elastic  joints  and  at  these  joints  the  lateral 
loads  are  applied.  The  sublayer  is  represented  by  spring 
elements  connected  with  the  slab  at  the  beam  joints. 


Expressing  the  equilibrium  equations  for  forces  and 
moments  at  the  joints,  a  set  of  simultaneous  equations  for 
the  whole  pavement  is  formed.  The  solution  of  these  equa- 
tions provides  the  deflections  at  the  nodes,  and  with 
these  the  complete  structure  can  be  solved. 

Huang  and  Wang  [35]  considered  the  presence  of  doweled 
joints  in  a  finite  element  method  essentially  the  same  as 
the  one  developed  by  Zienkiewicz  and  Cheung.  For  the 
analysis  of  doweled  concrete  pavements,  the  case  of  a  two 
slab  system  connected  with  dowel  bars  at  the  joint  was 
considered.  Then,  the  equations  for  the  system  were 
expressed  following  the  method  of  Zienkiewicz  and  Cheung 
and  considering  the  joint  as  a  discontinuity.  This  modif- 
ication destroys  the  symmetry  of  the  stiffness  matrix,  and 
therefore  the  total  matrix  has  to  be  stored.  In  this 
method  by  Huang  and  Wang,  the  sublayers  is  considered  as  a 
Winkler  ty pe  . 

Later  in  1974,  Huang  [36,37]  presented  an  iterative  method 
for  determining  the  stresses  and  deflections  in  pavements 
in  which  the  sublayer  behaves  as  an  elastic  solid* 
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Morever,     the  effect  of  temperature  gradients  in  the  slab 
was  considered. 

The  basic  element  is  the  plate  element  presented  by  Zienk- 
iewicz   and   Cheung.   Temperature  warping  is  introduced  by 

replacing  the  nodal  displacements,  ^a   ,  by  ^a>  =  <A  >  ,  in 
whi  ch: 


(2.26) 


And 


c.  =  ait   27 


(2.27) 


Where,  c  is  the  amount  of  warping  at  node  i,  a  is  the 
coefficient  of  thermal  expansion  of  concrete,  At  is 
the  temperature  differential,  between  top  and  bottom  of 
slab  d  is  the  distance  from  node  i  to  center  of  slab, 
and    t   is  the  slab  thickness. 

In  this  model,  the  resultant  set  of  simultaneous  equations 
is  not  symmetric.  To  solve  this  system,  Huang  developed  an 
iterative  scheme  to  make  the  stiffness  matrix  banded, 
facilitating  its  storage  and  solution. 

In  the  method  presented  by  Huang  and  Wang  ,  [36,37]  the 
effect  of  partial  contact  between  the  slab  and  the  sub- 
layer is  introduced  by  simply  deleting  the  reactive 
forces  at   the   nodes  assumed  not  in  contact.  However,  in 
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actual  conditions  a  node  not  in  contact  with  the  subbase 
may  come  in  contact  after  loading,  and  vice  versa.  To  take 
this  into  account,  the  same  investigators,  later  in  1974, 
presented  a  general  method  to  consider  both  full  and  par- 
tial contact.  Assuming  a  Winkler  foundation,  two  cases  of 
partial  contact  were  considered.  In  the  first  case,  there 
is  no  pumping  or  plastic  deformation  in  the  sublayer. 
Therefore,  the  springs  for  the  sublayer  elements  are 
assumed  of  equal  length.  In  the  second  case,  pumping  or 
plastic  deformation  can  be  introduced  by  using  springs  of 
unequal  length.  In  this  way,  voids  can  be  preassigned  to 
areas  where  excessive  deformation  in  the  sublayer  occurs. 
Temperature  is  considered  in  the  same  way  as  in  the  method 
presented  by  Huang   [36,37], 


Another  improvement  in  the  finite  element  modeling  of 
rigid  pavements  was  introduced  by  Tabatabaie  and  Barenberg 
[9]  .  They  developed  a  two  dimensional  finite  element 
model,  called  ILLISLAB,  based  on  a  plate  element  similar 
to  that  developed  by  Zienkiewicz  and  Cheung.  The  concrete 
slab  is  considered  in  the  same  way  as  in  Zienki ewi cz ' s 
model,  but  additionally  a  second  sublayer  can  be 
integrated  into  the  system.  The  slab  and  top  sublayer  can 
be  considered  to  be  perfectly  bonded  or  not  bonded.  per- 
fect bond  is  assumed  between  the  concrete  slab  and  the  top 
sublayer,  then  an  equivalent  transformed  section  is  used. 
This   transformed  section  is  treated  as  the  original  plate 
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element.  In  the  case  of  no  bond,  the  stiffness  matrix  for 
the  bottom  sublayer  and  slab  are  simply  added  to  each 
o the  r . 

Dowel  bars  are  modeled  by  a  bar  element  with  two  degrees 
of  freedom  at  each  node.  Thereby  the  dowels  can  transfer  a 
vertical  force  and  a  moment.  Aggregate  interlock  and  key- 
way   are  modeled  by  spring  elements  able  to  transfer  shear 

The  overall  structural  stiffness  matrix  is  obtained  by 
superimposing  the  matrix  of  the  individual  elements.  Gen- 
eralized deflections  and  stresses  are  calculated  by  solv- 
ing the  set  of  simultaneous  equations  : 


{■}  ■  WW 


(2.28) 


where  {  f  }  represents  the  equivalent  nodal  forces  for  a 
uniformly  distributed  load  over  a  rectangular  section  of 
the  concrete  slab.  Values  of  {a}  are  the  resultant  nodal 
displacements,  and  [K]  is  the  overall  structural  stiffness 
matrix.  In  this  method,  temperature  effects  are  not  con- 
sidered. 


Similar  models  have  been  presented  by  Tayabji  and  Colley 
[39],  and  Chou  [40],  As  in  the  case  of  the  model 
developed  by  Tabatabaie  and  Barenberg,  the  basic  element 
is  the  rectangular  thin  plate  element  given  by  Zienkiewicz 
and  Cheung.   The  slab  is  divided  into   rectangular   finite 
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elements  with  three  degrees  of  freedom  per  node.  The  sub- 
layer reactive  forces  act  at  the  nodes.  The  stiffness 
matrix  for  the  plate  and  sublayer  are  simply  added  to 
each  other.  In  the  model  by  Tayabji  and  Colley,  called 
JSLAB,  the  sublayer  is  assumed  as  a  Winkler  model  acting 
as  a  series  of  uniformly  distributed  spring  elements.  In 
this  model,  dowels  are  represented  by  beam  elements  with 
two  degrees  of  freedom  per  node:  vertical  displacement  and 
rotation.  Aggregate  interlock  and  keyways  are  represented 
by  spring  elements  by  means  of  which  only  vertical  load  is 
transferred  across  the  joint. 

In  the  JSLAB  model,  a  linear  temperature  gradient  is  con- 
sidered by  application  of  moment  along  slab  edges.  The 
moment  per  unit  width  is: 

3 
Ea  AT   Yj  (2.29) 

where  t  is  the  slab  thickness,  a  is  the  coefficient  of 
thermal   expansion,    AT   is  the  temperature  gradient,  and 

E    is  the  modulus  of  elasticity  of  concrete.   An   itera- 

c  J 

tive  analysis  is  used  to  establish  support  conditions.  For 
each  iteration,  a  check  of  contact  between  slab  and  sub- 
layer is  accomplished,  making  the  subgrade  modulus  equal 
to  zero  where  contact  is  lost. 

The  models  developed  by  Chou  [40],  called  WESLIQID  and 
WESLAYER,  are  also  based  on  the  plate  element  developed  by 
Zienkiewicz   and   Cheung.    WESLAYER   assumes    that    the 
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sublayer  acts  as  a  Winkler  foundation,  while  the  WESLAYER 
considers  it  as  an  elastic  layered  solid.  The  Winkler 
foundation  is  modeled  by  spring  elements.  In  the  case  of 
the  elastic  foundation  (WESLAYER),  Burmister's  equation  is 
used  to  compute  the  stiffness  matrix  for  the  sublayer.  In 
this  case,  the  combination  of  the  plate  and  sublayer 
matrices  gives  a  non  banded  matrix.  To  solve  this,  Chou 
applies  the  same  iterative  scheme  presented  by  Huang 
[36,37].  The  two  programs  proceed  in  an  iterative  manner, 
assuming  at  the  beginning  full  contact  between  slab  and 
sublayer,  except  at  the  nodes  where  gaps  are  preassigned. 
Sublayer  contact  conditions  and  convergence  of  deflec- 
tions along  joints  are  checked  after  each  iteration.  Con- 
nections can  be  considered  as  doweled,  aggregate  inter- 
lock, or  keyway.  For  dowels,  bar  elements  with  shear  and 
moment  are  used.  For  aggregate  interlock  and  keyways,  only 
shear  is  transferred  using  spring  elements  for  modeling. 


2_-2_  Existing  Computer  Models 
Several  computer  models  for  analyzing  concrete  pavements 
have  been  developed  since  the  1960's.  The  models 
ILLISLAB,  by  Tabatabaie  and  Barenberg  [9]  ;  JSLAB,  by  Tay- 
abji  and  Colley  [39]  ;  and  WESLIQID  and  WESLAYER,  by  Chou 
[40]  are  principal  examples.  The  assumptions,  modeling 
techniques,  and  applications  are  briefly  illustrated  in 
this  chapter.  The  models  ILLISLAB  and  JSLAB  are  opera- 
tional at  Purdue  University  at  the  present  time. 
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2.2.1  ILLISLAB 
Assumptions.-  In  this  model,  the  pavement  system  is 
divided  into  the  following  finite  elements:  rectangular 
plates  for  the  slab,  base  or  subbase,  and  overlay;  spring 
elements  for  the  subgrade;  spring  elements  for  joints  with 
aggregate  interlock  and  keyway;  and  bar  elements  for  the 
dowel  bars.  The  plate  element  uses  the  classical  theory 
of  medium  thick  plates,  in  which  lateral  forces  are  sup- 
ported by  the  flexural  strength  of  the  plate.  Lateral 
stresses  and  strains,  as  well  as  transverse  shear  deforma- 
tions are  neglected.  Both  perfect  bond  and  no  bond  between 
concrete,  base  or  subbase,  and  overlay  can  be  assumed. 
Kirchhoff  hypotheses  are  met  for  the   plate  element. 

The  subgrade  is  considered  a  Winkler  type.  It  behaves  as  a 
dense  liquid  for  which  the  deformation  at  a  certain  point 
depends  only  on  the  stress  at  the  point.  The  stress- 
deformation  relation  in  the  subgrade  is  represented  by  k 
,  the  modulus  of  subgrade  reaction.  The  values  for  this 
parameter  can  be  varied  from  node  to  node  and  can  be  taken 
as  zero  at  points  with  no  support . 


The  dowel  bars  behave  as  linearly  elastic  bars  located  at 
the  neutral  axis  of  the  plates.  These  bars  transfer  shear 
force  and  bending  moment  from  the  slab  on  one  side  of  the 
joint  to  the  slab  on  the  opposite  side. 
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For  the  case  of  aggregate  interlock,  and  keyway,  the  load 
Is  transferred  by  shear.  This  transfer  is  accomplished  by 
a  linearly  elastic  spring  element. 

The  slab  is  modeled  by  a  plate  element  with  three  degrees 
of  freedom  per  node.  These  are  a  vertical  force,  a  rota- 
tion around  the  x  axis,  and  a  rotation  around  the  y  axis. 
Corresponding  to  these  degrees  of  freedom,  there  are  three 
nodal  forces:  a  vertical  force,  and  two  moments  around  x 
and  y  axes.  For  the  case  of  perfect  bond  between  layers, 
complete  strain  compatibility  between  layers  is  assumed. 
The  plate  considered  in  the  analysis  has  a  transverse  sec- 
tion equal  to  the  transformed  section  resulting  from  the 
base,  concrete,  and  overlay  sections.  For  the  case  of  no 
bond  between  layers,  the  stiffness  matrices  of  the  con- 
crete, the  base  or  subbase  and  the  overlay  are  simply 
superimposed  by  addition.  The  properties  of  the  plate  can 
be  varied  from  element  to  element. 


The  Winkler  foundation  is  modeled  by  a  series  of  spring 
elements  distributed  uniformly  under  the  slab.  The  spring 
elements  have  one  degree  of  freedom  per  node:  a  vertical 
displacement.  Corresponding  to  this,  there  is  a  vertical 
force.  The  force  at  each  element  is  given  by  the  deforma- 
tion of  the  element  times  the  spring  constant.  This 
spring  constant  is  the  modulus  of  subgrade  reaction  for 
the  subgrade  at  the  location  of  the  element.  The  value  for 
this  parameter  can  be  changed  from  node  to  node. 
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Dowels  are  modeled  by  bars  with  two  degrees  of  freedom  per 
node.  These  degrees  of  freedom  are  a  vertical  displacement 
and  a  rotation  around  a  horizontal  transverse  axis.  There- 
fore, the  bar  is  able  to  transfer  a  vertical  shear  force 
and  a  moment.  In  addition,  some  looseness  between  the 
dowel  bar  and  the  surrounding  concrete  is  modeled  by  a 
spring  element  whose  constant  is  the  so  called  concrete- 
dowel  interaction  factor. 

Aggregate  interlock  and  keyway  joints  are  modeled  by 
spring  elements.  These  elements  have  a  vertical  displace- 
ment as  single  degree  of  freedom.  Vertical  forces  are 
transferred  across  the  joints  by  the  springs. 

Applications .-  This  model  can  be  used  in  the  analysis  of 
cracked  and  jointed  concrete  pavements  with  base  or  sub- 
base,  overlay,  and  subgrade.  Up  to  nine  slabs  can  be 
analyzed  with  one  longitudinal  joint  and  two  transverse 
joints.  Connections  can  be  with  aggregate  interlock,  with 
keyway,  or  with  dowel  bars.  Either  perfect  bond  or  no  bond 
can  be  assumed  between  the  concrete,  the  overlay,  or  the 
sub laye  r s . 


It  is  also  possible  to  analyze  concrete  shoulders  with  and 
without  tie  bars,  as  well  as  concrete  pavements  with  vary- 
ing thicknesses  and  modulus  of  elasticity,  and  with  vary- 
ing modulus  of  support  for  the  subgrade. 
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2.2.2  JSLAB 
Assumptions .-  This  model  is  similar  to  the  ILLISLAB  model. 
The  pavement  structure  is  divided  into  the  same  types  of 
elements:  the  slab,  modeled  by  rectangular  plate  elements; 
the  joints,  modeled  by  thick  beams  or  spring  elements;  and 
the  sublayers,   modeled  by  spring  elements. 

The  plate  elements  for  the  slab  can  be  composed  of  one  or 
two  layers  with  or  without  bond  between  them.  The  stiff- 
ness matrix  for  the  plate  is  developed  using  the  classical 
theory  of  plates  with  small  deflections,  assuming  the 
material  to  be  homogeneous  and  elastic.  Strains  in  the 
vertical  direction,  as  well  as  strains  resulting  from 
shear  in  vertical  planes,  are  neglected.  For  the  case  of 
bond  between  layers,  the  transformed  section  concept  is 
used;  for  the  case  of  no  bond,  the  stiffness  matrices  of 
the  layers  are  superimposed. 

With  the  JSLAB  model,  it  is  possible  to  include  the 
effect  of  a  linear  temperature  gradient  in  the  concrete 
slab.  This  effect  is  considered  by  applying,  along  the 
slab  edges,  a  moment  whose  magnitude  depends  on  the  tem- 
perature gradient,  the  modulus  of  elasticity,  and  the 
coefficient  of  thermal  expansion.  The  model  uses  an  itera- 
tive process  to  establish  the  support  conditions  due  to 
curling.  At  locations  where  loss  of  support  occurs, 
the  subgrade  modulus  is  made  zero. 
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Dowel  bars  are  assumed  to  behave  like  thick  beams  able  to 
transfer  shear  and  moment  across  the  joints.  The  stiffness 
matrices  for  these  beams  is  modified  to  account  for  the 
deformation  that  occurs  in  the  concrete  due  to  the  action 
of  the  dowel.  This  is  similar  to  the  concrete-dowel 
interaction  concept  in  the  1LLISLAB.  Aggregate  interlock 
and  keyway  are  represented  by  spring  elements  transferring 
only  shear  force.  These  elements  behave  like  linearly 
elas  tic  springs. 

As  in  the  ILLISLAB  model,  the  subgrade  is  assumed  to 
behave  as  a  Winkler  sublayers,  for  which  stresses  and 
deflections  at  different  locations  are  independent.  The 
subgrade  is  modeled  by  uniformly  distributed  springs  with 
equal  or  different  elasticity  constants.  These  springs, 
with  one  degree  of  freedom,  can  have  a  spring  constant 
equal  to  zero  at  locations  where  no  support  occurs. 


Applications .-  This  model  can  be  used  in  the  analysis  of 
jointed  concrete  pavements  with  up  to  nine  slabs.  Joints 
can  be  longitudinal  or  transverse  with  dowels,  aggregate 
interlock,  or  keyway.  Any  configuration  of  static  vertical 
loads  can  be  applied  as  either  distributed  or  concentrated 
loads.  It  is  also  possible  to  apply  vertical  deflections 
or  moments  at  various  nodes. 
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The  model  can  analyze  pavements  when  temperature  and  mois- 
ture gradients  exist  in  the  concrete  slab.  The  effect  of 
a  moisture  gradient  is  obtained  by  applying  an  equivalent 
temperature  differential.  The  slab  can  be  considered  with 
one  or  two  layers  representing  the  concrete  slab  and  base 
or  subbase.  Bond  or  no  bond  can  be  considered.  The  sub- 
layers can  be  assumed  to  have  different  properties  at  dif- 
ferent locations,  making  it  possible  to  represent  voids 
underneath  the  slab  by  making  the  subgrade  modulus  equal 
to  zero  at  the  nodes  where  no  contact  occurs. 


2.2.3  WESLIQID  and  WESLAYER. 
Assumptions .-  These  finite  element  methods  are  also  based 
on  the  classical  theory  of  thin  plates  with  small  deforma- 
tions. The  pavement  structure  is  divided  into  the  same 
elements  as  those  used  in  1LLISLAB  and  JSLAB:  slab,  sub- 
layer or  sublayers,  and  joints.  The  slab  is  modeled  by 
the  rectangular  plate  element  developed  by  Zienkiewicz  and 
Cheung.  The  basic  difference  is  in  their  modeling  of  the 
sublayers,  which  are  considered  as  an  elastic  layered 
solid  in  the  WESLAYER  model;  however,  the  WESLIQID  model 
assumes  a  Winkler  sublayers 

In  the  WESLIQID  model,  the  sublayers  are  considered  as  a 
dense  liquid  material  modeled  by  a  series  of  uniformly 
distributed  elastic  springs.  The  forces  due  to  the  reac- 
tion at  the  sublayers  are  simply  added  to  the  forces  of 
the  plate  at  any  given  node.  These   reactions  have  only   a 
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vertical  component  which  is  added  to  the  vertical  com- 
ponent of  the  nodal  forces  in  the  plate. 

In  the  WESLAYER  model,  the  sublayers  are  considered  as  a 
layered  elastic  solid,  in  which  the  deflection  at  any 
given  point  depends,  not  only  on  the  forces  at  the  node, 
but  also  on  the  forces  and  deflections  at  other  nodes.  The 
stiffness  matrix  for  the  sublayers  is  obtained  by  invert- 
ing the  flexibility  matrix  formed  using  the  Burmister 
equation. 

The  WESLIQID  and  WESLAYER  models  can  consider  the  effect 
of  a  linear  temperature  gradient  in  the  concrete  slab.  The 
initial  deflection  due  to  the  temperature  gradient  is  com- 
puted using  the  coefficient  of  thermal  expansion,  the  tem- 
perature gradient,  the  thickness,  and  the  distance  to  the 
center  of  the  slab  where  the  deflection  is  zero.  The 
equilibrium  equations  for  the  nodes  are  formed  including 
the  temperature  deflections.  Full  contact  or  partial  con- 
tact between  slab  and  sublayer   can  be  assumed. 

The  analysis  of  the  joints  is  accomplished  by  combining 
different  amounts  of  shear  and  moment  transfers.  The  two 
models  have  three  options  for  specifying  shear  transfer 
and  one  for  moment.  Specification  of  shear  transfer  can  be 
accomplished  by:  a  shear  transfer  efficiency  factor, 
defined  by  the  ratio  of  the  vertical  deflections  between 
loaded  and  unloaded  slabs  along  the  joint;    a    spring 
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constant  modeled  by  imaginary  springs  along  the  joints;  and 
by  diameter  and  spacing  of  dowel  bars,  which  is  applicable 
only  when  steel  bars  are  the  only  means  of  shear  transfer. 
It  is  also  possible  to  specify  the  conditions  of  looseness 
of  the  dowel  bars  by  a  dowel  support  factor. 

Moment  transfer  across  joints  or  cracks  is  accomplished  in 
a  two  step  procedure.  First,  the  moments  at  nodes  along 
the  crack,  or  joint  are  calculated.  Then  these  moments  are 
multipied  by  the  efficiency  of  moment  transfer,  and 
assigned  to  each  slab  as  externally  applied  moments.  In  a 
second  step  the  applied  moments  are  included  in  the 
analysis.  The  definition  of  moment  transfer  is  based  on 
the  rotations  of  the  slabs  adjacent  to  the  joint  or  crack. 

Applicat  ions .-  These  finite  element  models  can  be  used  in 
the  determination  of  stresses  and  deflections  in  concrete 
pavements  with  cracks  and  joints.  Analysis  can  be  carried 
out  for  pavements  with  slabs  made  up  of  two  layers  of 
materials  with  different  properties.  The  concrete  slab  can 
be  either  plain  concrete  or  reinforced  concrete.  The  pave- 
ment sublayers  can  be  considered  either  as  a  Winkler  sub- 
layers or  as  a  layered  elastic  solid.  Full  or  partial 
support  of  the  concrete  slab  can  be  assumed.  Voids  can  be 
introduced  by  assuming  different  values  for  the  properties 
of  the  sublayers.  The  joints  can  be  doweled  or  connected 
with  other  load  transfer  devices  like  keyways  and  aggre- 
gate interlock.   The  WESLAYER  model  is  limited  to  only  two 


33 

slabs  because  of  the  large  computer  storage  required.  The 
analysis  can  be  done  considering  external  distributed  or 
concentrated  loads  as  well  as  thermal  effects.  With  all 
these  options,  analysis  can  be  performed  for  different 
types  of  concrete  pavements  and  concrete  shoulders. 

—•—  Brief  Summary  of  Design  Methods 
Design  procedures  for  concrete  pavements  have  been  based 
on  a  stress  fatigue  criterion.  Recently,  the  inclusion  of 
an  erosion  criteria  has  been  proposed  by  Packard  and  Tay- 
abji  [8].  In  the  following,  an  overview  of  the  stress 
fatigue  design  methods  is  presented.  The  methods  developed 
by  the  Portland  Cement  Association  [41],  AASHTO  [42],  and 
the  ACI  [43]  are  discussed  as  follows: 


2.3.1  Portland  Cement  Association  Method 
The  basic  concept  that  has  been  used  for  design  of  both 
airfield  and  highway  pavements  is  the  concept  of  fatigue 
of  concrete.  In  accordance  with  this  concept,  the  number 
of  load  applications  that  a  concrete  beam  sustains  in 
flexure  depends  on  the  properties  of  the  concrete  and  on 
the  level  of  stresses.  The  Portland  Cement  Association 
obtained  a  table  which  relates  the  stress  level  with  the 
number  of  load  repetitions  needed  to  produce  failure  [41], 
This  table  shows  that  if  the  relative  stress  (stress 
divided  by  modulus  of  rupture)  is  less  than  0.51,  then 
concrete  is  able  to  sustain  an  unlimited  number  of   stress 
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repetitions  of  this  or  lesser  magnitude  without  failure. 
As  the  relative  stress  increases,  the  possible  number  of 
stress  repetitions  decreases. 


In  the  design  method  developed  by  the  Portland  Cement 
Association,  the  thickness  of  the  concrete  slab  is  depen- 
dent upon  the  magnitude  and  number  of  load  repetitions, 
the  modulus  of  rupture  of  concrete,  and  the  modulus  of 
subgrade  reaction.  The  number  of  load  repetitions  that  a 
pavement  sustains  during  its  service  life  is  estimated  and 
categorized  into  axle-load  groups.  Assuming  certain  thick- 
ness and  using  Westergaard's  equations,  the  stresses 
corresponding  to  each  of  the  load  categories  are  calcu- 
lated. Each  of  these  stresses  is  divided  by  the  modulus 
of  rupture  of  the  chosen  concrete  giving  the  relative 
stresses.  Using  the  PCA  table,  the  number  of  allowable 
load  repetitions  is  obtained  for  each  category.  The  per- 
centage of  actual  load  applications  with  respect  to  the 
allowable  number  of  load  applications  are  calculated,  and 
these  percentages  are  summed.  Theoretically,  the  design 
is  adequate  if  the  sum  of  these  percentages  is  less  than 
100  percent.  The  calculation  of  stresses  is  made  by  the 
Westergaard  equations  or  the  Pickett  and  Ray  tables  [7], 
Continuity  and  full  support  conditions  are  assumed. 
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2.3.2  AASHTO  Method 
The  AASHTO  method,  developed  after  the  AASHO  Road  Test,  is 
also  based  on  the  concept  of  concrete  fatigue.  Empirical 
equations  were  developed  from  the  Road  Test  to  represent 
the  relationship  between  loss  in  servi ciabi li ty ,  traffic, 
and  pavement  thickness.  In  this  method,  the  expected 
traffic  mix  is  converted  into  equivalent  18-kip  single 
axle  loads.  Given  these  equivalent  load  applications,  the 
desired  pavement  thickness  is  determined  on  the  basis  of 
the  desired  working  stresses  and  the  value  of  the  modulus 
of  subgrade  reaction,  using  the  empirical  relations 
obtained  from  the  Road  Test. 


The  total  number  of  equivalent  single-axle  load  applica- 
tions is  determined  by  dividing  the  traffic  mix  into  vari- 
ous weight  categories.  Each  category  is  multiplied  by  a 
load  factor  (these  load  factors  were  developed  from  the 
empirical  equations  of  the  Road  Test  to  have  equivalent 
effects  on  the  performance  of  the  pavement).  Originally, 
the  general  Road  Test  equations  were  obtained  for  certain 
specific  values  for  the  modulus  of  elasticity,  modulus  of 
rupture  of  concrete,  modulus  of  subgrade  reaction,  and 
environmental  conditions.  To  account  for  conditions  other 
than  those  that  existed  at  the  Road  Test,  the  equations 
were  modified  by  comparing  the  stresses  measured  on  the 
test  with  the  stresses  calculated  with  the  Westergaard, 
Spangler,  and  Pickett  equations. 
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2.3.3  American  Concrete  Institute  Method 
The  ACI  committee  325  [43]  presented  a  design  method  to 
determine  the  pavement  thickness  and  amount  of  reinforcing 
steel  for  continuously  reinforced  concrete  pavements.  This 
method  is  based  on  some  modifications  of  the  AASHTO 
method.  In  continuously  reinforced  concrete  pavements,  the 
slab  must  resist  stresses  and  deflections  produced  by 
applied  traffic  loads,  and  the  function  of  the  longitudi- 
nal steel  is  to  keep  the  cracks  in  concrete  tightly  closed 
so  that  there  is  an  effective  load  transfer  across  the 
cracks.  In  the  ACI  method  a  factor  for  considering  the 
load  transfer  in  the  thickness  design  is  included.  Addi- 
tionally, equations  to  determine  the  percentage  of  longi- 
tudinal stress  are  provided.  This  longitudinal  steel  is 
provided  to  control  cracking  due  to  shrinkage  and  tempera- 
ture. The  results  obtained  with  both  methods  are  compar- 
able, and  the  difference  depends  on  the  value  selected  for 
the  load  transfer  factor.  The  ACI  committee  325  recommends 
a  value  of  2.2  for  the  load  transfer  factor.  Assuming  this 
value,  the  results  obtained  with  the  AASHTO  and  the  ACI 
methods  are  similar. 


The  traditional  design  methods  have  been  based  on  the 
assumption  of  a  failure  due  to  concrete  fatigue.  However, 
in  fact,  most  of  the  failures  in  concrete  pavements  have 
occurred  due  to  a  different  mechanism  [2,8], 
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Most  of  the  failures  occurred  as  a  consequence  of  erosion 
of  the  sublayers.  Pumping,  faulting,  and  slab  breakup  are 
the  major  problems.  Therefore,  a  realistic  design  pro- 
cedure should  be  developed  considering  the  deterioration 
of  pavements  caused  by  erosion.  The  analysis  method  should 
also  reflect  this  failure  mechanism. 


2_.4_  Summary 
In  the  preceding  paragraphs,  the  development  of  the  tech- 
niques for  structural  analysis  of  concrete  pavements  has 
been  illustrated.  A  review  of  literature  published  since 
the  1900's  has  been  presented.  This  review  was  done  by 
dividing  the  available  information  into  three  parts: 
literature  related  to  the  analysis  with  closed  form 
methods;  literature  related  to  the  analysis  of  joints;  and 
literature  related  to  finite  element  methods.  Addition- 
ally the  computer  models  ILLISLAB,  JSLAB,  WESLIQID,  and 
WESLAYER  were  reviewed,  and  the  assumptions  and  applica- 
tions of  these  models  were  discussed. 
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CHAPTER  3  THEORETICAL  EVALUATION  OF  EXISTING 
METHODS  OF  ANALYSIS 


2»1_  Introduction 
In  order  to  determine  the  abilities  that  the  existing  com- 
puter models  have  for  analyzing  the  structural  performance 
of  rigid  pavements  with  pumping,  theoretical  and  numerical 
evaluations  of  such  models  were  carried  out.  At  the  same 
time,  the  characteristics  required  in  the  development  of 
an  improved  model  were  obtained. 

3  .2  Criteria   of  Model  Evaluation 
The  models  were  studied  in  accordance  with  three  different 
criteria   as  follows: 

1.  Evaluation   in   terms   of   their   suitabilities    for 
des  ign. 

2.  Evaluation  in  terms  of  the   modeling   of   the   actual 
behavior  of  rigid  pavements. 

3.  Numerical  Evaluation. 

The  first  two  criteria   are   discussed   in   this   chapter. 
Numerical  evaluation  is  given  in  Chapter  4. 
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3.2.1  Evaluation  for  Design  Purposes 
The  objective  of  the  design  process  is  to  obtain  a  pave- 
ment with  such  properties,  that  under  given  environmental 
and  traffic  conditions  provides  a  service  life  of  a  speci- 
fied number  of  years.  This  service  life  is  the  period  of 
time  between  the  time  at  which  the  pavement  is  opened  to 
traffic  until  a  functional  failure  occurs.  Yoder  and 
Witczack  [1]  define  functional  failure  as  the  state  of  a 
pavement  which  is  no  longer  capable  of  carrying  its 
intended  function  without  causing  discomfort  to 
passengers,  or  without  causing  high  stresses  in  the  vehi- 
cles that  pass  over  it.  The  generally  accepted  design 
procedures  are  based  on  the  principle  of  limiting  the 
flexural  stresses  in  the  concrete  slab  to  safe  values. 
With  this,  it  is  intended  to  avoid  fatigue  cracking  caused 
by  load  repetitions.  The  PCA  and  ASSHTO  design  methods 
provide  values  for  slab  thicknesses  such  that  for  the 
predicted  traffic,  the  number  of  stress  repetitions  does 
not  exceed  the  fatigue  resistance  of  concrete.  The  stress 
calculations  are  based  upon  the  assumption  of  full  support 
at  all  points,  and  at  all  instants  in  the  service  lives  of 
the  pavements. 


In  reality,  through  the  service  life  of  a  pavement  its 
stiffness  properties  change.  The  passage  of  certain  load- 
ing over  a  pavement  at  the  beginning  of  its  life,  produces 
an   effect  which  is  different  from  the  one  produced  by  the 
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same  loads  applied  after  several  years  of  service.  With 
time,  microcracks  can  develop  in  the  slab,  decreasing  its 
stiffness.  Additionally,  voids  and  loss  of  support  can 
develop  increasing  stresses  and  deflections. 

The  computation  of  stresses  and  deflections  in  the  struc- 
tural analysis  should  be  based  on  actual  conditions.  At 
the  beginning  of  the  service  life,  the  assumption  of  full 
support  is  acceptable.  However,  after  several  years  in 
service,  the  assumptions  of  full  support  and  initial 
strength  and  modulus  of  elasticity  are  not  valid. 

The  present  methods  of  design  have  been  based  on  the 
assumption  that  failure  in  rigid  pavements  occurs  due  to 
fatigue  in  the  concrete  slab.  However,  it  has  been 
observed  [1,2,4]  that  most  of  the  failures  in  rigid  pave- 
ments are  preceded  by  pumping  of  the  material  underneath 
the  slab.  Packard  and  Tayabji  [8]  have  proposed  a  design 
method  which  considers  the  damage  in  the  material  under 
the  slab.  Generally,  an  improved  method  of  design  should 
consider  fatigue  in  the  slab,  as  well  as  the  erosion  in 
the  sublayers . 

As  part  of  the  design  process,  the  structural  analysis 
must  be  capable  of  determining  the  response  (in  terms  of 
stresses,  deflections,  cracking,  and  damage  in  general  ) 
of  the  pavement  structure  to  environmental  and  traffic 
actions.   With  given  environment  and   traffic   conditions, 
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the  designer  proposes  certain  geometry  and  material  pro- 
perties for  the  pavement.  Next,  the  response  of  the  pro- 
posed pavement  is  evaluated  in  terms  of  service  life. 
Calculating  the  initial  and  maintenance  costs  of  various 
alternatives,  the  designer  chooses  the  best  alternative, 
engaging  in  a  "trade-off"  process. 

In  order  to  perform  the  trade-off  process  with  rigid  pave- 
ments, the  designer  must  be  able  to  consider  all  the  prac- 
tical types  of  pavements.  These  types  of  pavements  are: 
1)  plain,  jointed  pavements  (PJP);  2)  continuously  rein- 
forced concrete  pavements  (CRCP);  3)  jointed,  reinforced 
concrete  pavements  (JRCP);  and  prestressed  concrete  pave- 
ments (PCP).  Yoder  and  Witczak  [1]  classify  the  joints 
as:  1)  contraction  joints;  2)  expansion  joints;  3)  con- 
struction joints;  and  4)  warping  joints.  In  accordance 
with  their  performance,  the  joints  can  be:  1)  plain, 
dummy-groove  joints;  2)  doweled,  dummy-groove  joints;  3) 
doweled  joints;  4)  keyed  joints;  and  5)  tied  joints. 


The  section  of  a  rigid  pavement  consists  of  the  concrete 
slab  and  several  sublayers  (Fig.  1.1).  These  layers, 
which  constitute  the  sublayers  of  the  slab,  can  be  com- 
posed of  a  stabilized  or  unstabilized  base  or  subbase,  and 
a  subgrade.  For  these  components  of  the  pavement,  many 
different  combinations  of  thicknesses  and  strengths  can  be 
used.  In  the  concrete  slab,  the  thickness  and  strength  can 
be   varied.   Values   for   highway  slab  thickness  vary  from 
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152mm  (6  in.)  to  305mm  (12  in.),  and  for  the  compressive 
strength  from  20.7  MPa  (3000  psi)  to  34.5  MPa  (  5000  psi) 
or  higher.  For  the  base  or  subbase  ,  the  desired  strength 
can  be  obtained  by  varying  the  gradation  of  the  granular 
material  as  well  as  by  varying  the  compacting  energy. 
Traditionally,  the  stiffness  properties  of  the  sublayers 
are  measured  by  the  value  of  the  modulus  of  subgrade  reac- 
tion. Recently  however,  it  has  been  proposed  [44]  that  a 
resilient  modulus  be  used.  The  traditional  modulus  of 
subgrade  reaction  is  obtained  from  a  static  test.  This 
test  does  not  reproduce  the  behavior  of  the  sublayer  when 
a  moving  load  (traffic  load)  is  applied.  By  contrast,  the 
resilient  modulus,  which  is  obtained  from  an  impact  test 
and  is  strain  dependent,  represents  more  accurately  that 
behavior.  Typical  values  for  the  subgrade  modulus  vary 
from  27.1  MPa/m  (100  pci  )  ,  for  a  fine  grained  material,  to 
135.5  MPa/m  (500  pci)  or  more  for  a  granular  well  graded 
material.  For  the  resilient  modulus,  values  can  vary  from 
81.3  MPa/m  (300  pci)  for  a  soft  fine  material,  to 
271. MPa/m  (1000  pci)  or  more  for  a  stiff  granular 
material.  Most  of  the  time  a  base  or  subbase  is  pro- 
vided. The  modulus  of  elasticity  of  a  stabilized  sublayer 
may  be  from  6900.  MPa  (1.x  106  psi  )  to  13800  MPa  (2.x 
10  psi).  For  nons t ibilized  sublayers  these  values  may 
be  half  of  the  values  for  stibilized  ones. 
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In  addition  to  considering  the  material  properties  and  the 
geometric  dimensions  of  the  pavement,  the  designer  has  to 
consider  the  climatic  and  traffic  conditions.  The  amount 
of  rainfall  and  the  daily  and  yearly  temperature  varia- 
tions should  be  estimated.  Traffic  can  be  quantified  by 
the  number  of  vehicles  of  different  weights  that  are 
expected  to  pass  over  the  pavement  during  its  service 
life.  When  pumping  is  expected,  certain  special  condi- 
tions should  be  considered  as  well.  It  is  important  to 
study  the  erosion  susceptibility  and  potential  for  forma- 
tion of  voids  in  the  sublayers.  Erosion  of  sublayers  and 
subsequent  formation  of  voids  can  modify  the  support  con- 
ditions in  the  slab,  causing  increase  in  the  stresses  and 
deflect  ions  . 


Finally,  for  the  trade-off  process,  the  estimation  of  the 
possible  maintenance  costs  is  important.  The  economical 
analysis  of  the  pavement  alternatives  is  on  the  basis  of 
the  total  cost,  which  includes  the  initial  cost  plus  the 
maintenance  cost.  Therefore,  although  the  estimation  of 
the  initial  cost  is  relatively  simple,  since  the  geometry 
and  material  properties  are  known,  the  estimation  of  the 
maintenance  cost  is  more  complex,  since  the  actual  mainte- 
nance work  is  not  known  in  advance.  The  allocation  of  cer- 
tain annual  amounts  of  money  for  maintenance,  regardless 
of  the  differing  service  conditions  of  the  pavements  is 
not   realistic.  Instead,  cost  should  be  allocated   to  each 
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alternative  pavement  in  accordance  with  the  expected   con- 
ditions  that   the   specific   pavement   will  achieve  after 
different  time  periods.    Any   estimation   of   the   future 
state   of   the   pavements,   given  the  traffic  and  climatic 
conditions,  will  be  of   help   in   the   estimation   of   the 
maintenance   costs.    The  designer  must  be  able  to  predict 
the  service  life   of  each    of   the   proposed   alternative 
pavements.   Moreover,   he   has   to   predict   the  amount  of 
maintenance  that  each  alternative  pavement  will  need  after 
certain   several  years  in  service.  To  do  so,  the  amount  of 
damage  that  can  occur  in  each   alternative   design,   given 
the   expected   traffic,   should  be  predicted.   This  damage 
can  be  evaluated  in  terms  of  permanent  deflections,  crack- 
ing, or  faulting. 

3.2.2  Evaluation  in  Terms  of  Pavement  Behavior 
Rigid  pavements  are  made  up  of  a  Portland  cement  concrete 
slab  overlying  a  series  of  less  stiff  layers.  The  essen- 
tial difference  between  flexible  and  rigid  pavements  is 
the  manner  in  which  they  distribute  the  loads  over  the 
subgrade.  In  rigid  pavements,  due  to  the  relatively  high 
stiffness  of  the  concrete  slab,  loads  are  distributed  over 
a  wide  area  of  the  subgrade.  Therefore,  the  major  struc- 
tural capacity  is  provided  by  the  concrete  slab.  This 
slab  can  be  a  continuous  slab,  for  the  case  of  continu- 
ously reinforced  concrete  pavement,  or  formed  by  a  series 
of  slabs  separated  by  joints  and   interconnected   by   load 
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transfer  devices  [1].  Generally,  a  subbase  is  placed 
between  the  subgrade  and  the  slab  to  control  pumping, 
frost  action,  and  shrink  and  swell  acting  of  the  subgrade. 
Additionally,  it  provides  drainage.  The  subbase  is  built 
beyond  the  edge  of  the  pavement  so  that  loads  applied  at 
the  edge  of  the  pavement  are  distributed  over  a  wider 
area.  To  prevent  pumping,  subbases  are  made  up  of  granular 
materials,  or  materials  stabilized  with  lime,  portland 
cement,  or  asphalt. 

The  subgrade,  usually  constructed  by  compacting  the 
materials  existing  at  the  place  of  construction,  provides 
the  ultimate  support  of  the  road.  Since  the  subgrade  can 
have  relatively  low  stiffness,  the  slab  and  the  subbase 
must  distribute  the  applied  forces  over  a  wider  area  to 
have  low  stresses  at  the  subgrade  level. 

Damage  to  the  pavement  structure  can  be  caused  by  several 
factors.  In  general,  damage  is  caused  by  environmental 
factors  such  as  presence  of  water  and  temperature  varia- 
tions, and  traffic  factors,  such  as  repetitive  stresses 
and  deformations. 


Water,  either  from  rainfall  or  from  underground,  can  enter 
into  the  sublayer  materials  changing  their  strengths  and 
supporting  capacities.  Surface  infiltration  of  water 
through  cracks  or  joints,  combined  with  repetitive  appli- 
cation of  loads,  can  cause  pumping.  Excessive   water    may 
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erode  the  sublayer  materials  when  pressure  applied  upon 
the  slab  by  traffic  forces  water  to  be  ejected  at  high 
speed,  eroding  the  subbase  and  creating  voids.  Further 
load  repetitions,  after  some  loss  of  support  has  occurred, 
produces  cracking  in  the  slab  and  eventually,  functional 
failure. 

Temperature  variations  can  produce  uniform  longitudinal 
stresses  at  the  cross  section  of  the  pavement,  or  flexural 
stresses  caused  by  temperature  differentials  across  the 
thickness  of  the  slab.  Longitudinal  tensile  stresses 
develop  when  the  concrete  slab  is  cooled  down  and  its 
contraction  is  prevented  by  the  friction  with  the  subbase. 
Longitudinal  compressive  stresses  occur  when  the  tempera- 
ture in  the  slab  increases  and  the  expansion  of  the  slab 
is  prevented  by  friction  with  the  subbase  or  by  adjacent 
expanded  slabs. 

Warping  stresses  result  from  uneven  temperature  distri- 
bution over  the  cross  section  of  the  slab.  The  same  effect 
can  occur  as  a  result  of  uneven  distribution  of  moisture 
in  the  slab.  Differential  drying  shrinkage  can,  similarly, 
cause  flexural  stresses,  but  with  magnitudes  smaller  than 
those  due  to  temperature  differential. 


One  of  the  major  causes  of  damage  in  the  pavement  is  the 
passage  of  vehicles.  Different  kinds  of  traffic  loads  are 
applied  to  concrete  pavements  during  their  service   lives. 
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These  loads  vary  in  magnitude  and  in  point  of  application. 
Consideration  of  traffic  loads  for  highways  is  generally 
based  on  the  determination  of  the  equivalent  damage 
effects  of  all  the  vehicle  types  expressed  in  terms  of  the 
number  of  load  repetitions  of  a  standard  vehicle.  The 
relative  amount  of  damage  caused  by  a  particular  type  of 
vehicle  is  related  to  the  damage  caused  by  the  standard 
vehicle  by  using  equivalence  factors    [42], 

Loads  applied  by  traffic  vary  in  location  and  magnitude. 
A  plot  of  the  distribution  of  load  application  with 
respect  to  the  distance  to  the  edge  of  the  pavement,  shows 
that  most  of  the  loadings  occur  within  fairly  well  defined 
areas.  Loads  of  the  same  magnitude  but  applied  at  dif- 
ferent locations  produce  different  effects  in  the  pave- 
ment. Therefore,  critical  load  locations  should  be  con- 
sidered; in  general,  these  are  at  the  corner  of  the  slab, 
at  joints,  or  at  the  outer  edge  of  the  slab. 

Similarly,  the  magnitudes  of  the  applied  loads  are  vari- 
able. Vehicles  of  different  weights  and  axle  configura- 
tions transit  over  highway  pavements.  Therefore,  a  clas- 
sification of  vehicles  in  groups  with  similar  weights  and 
axles  is  necessary. 


When  a  load  is  applied,  the  pavement  responds,  deflecting 
and  transferring  the  pressure  over  its  supporting 
material.  Depending  on  the  flexural  stiffness  of  the   slab 
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and  the  location  of  the  load,  the  slab  acquires  certain 
deflected  shapes,  producing  flexural  stresses  along  its 
cross  section  and  bearing  stresses  on  the  sublayers.  The 
pressure  applied  on  the  slab  is  distributed  through  a 
larger  area  at  the  slab-subbase  interface.  The  pressure  is 
furthermore  reduced  as  it  is  distributed  through  the 
thicknesses  of  the  sublayers.  When  a  vehicle  passes  over  a 
point  in  the  pavement,  the  slab  is  under  flexural  and 
shear  stresses,  while  the  sublayers  are  under  the  pressure 
of  the  deformed  slab.  The  slab  and  sublayers  are 
relieved  of  these  stresses  when  the  vehicle  moves  away. 
Therefore,  with  vehicles  passing,  the  pavement  structure 
at  specific  points  is  subjected  to  repetitive  stressing. 

In  the  slab,  the  cumulative  effect  of  loads  causes  micro- 
cracking  and  reduction  of  strength.  The  concrete  strength 
in  general  and  the  flexural  strength  in  particular  are 
reduced.  As  the  strength  is  reduced,  the  flexural  stiff- 
ness of  the  slab  is  reduced;  thus  a  higher  pressure  is 
transmitted  to  the  sublayers,  producing  larger  deflec- 
t  ions  . 


The  traffic  loads  also  damage  the  sublayers  material.  The 
material  immediately  underneath  the  loaded  slab  is  sub- 
jected to  repetitive  application  of  pressure.  Each  load 
application  produces  a  small  permanent  deformation  in  the 
subbase  and  subgrade.  Even  though  compaction  and  stiffen- 
ning   of  the  material  can  occur,  the  presence  of  excessive 
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water  may  damage  the  structure  of  the  sublayer  material, 
thereby  weakening  the  sublayers.  Additionally,  if  water 
accumulates  in  spaces  between  the  slab  and  sublayers  pro- 
duced by  warping  or  compression,  the  sudden  application  of 
a  load  ejects  the  water.  The  water  moving  at  high  speed 
causes  erosion  and  additional  enlargement  of  the  void. 
This  condition  is  especially  critical  at  the  edge  and 
joints  of  the  slab,  where  pumping  can  occur.  Moreover,  at 
these  locations,  the  openings  allow  the  penetration  of 
water,  and  the  deflections  are  larger  because  of  the 
discontinuity. 

Load  transfer  devices  are  often  provided  at  the  joints  to 
reduce  the  deflections  caused  by  traffic  loads.  Dowel  bars 
are  placed  to  transfer  part  of  the  load  across  the  joint. 
These  bars  cause  the  adjacent  slab  to  deflect  a  propor- 
tion of  the  deflection  of  the  loaded  slab.  These  dowel 
bars,  with  one  end  cast  in  one  slab  and  the  other  end 
embedded  with  a  lubricant  in  the  other  slab,  transfer 
shear  and  moment  actions  across  the  joint.  The  efficiency 
of  the  bar  is  not  100%,  i.e.  deflections  at  both  sides  of 
the  joint  are  not  equal.  Part  of  the  loss  of  load  transfer 
is  due  to  the  deformation  of  the  bar  itself,  and  part  is 
due  to  deformations  in  the  concrete-bar   interfaces. 


In  tied  joints,  the  slabs  are  separated  by  openings  or 
weak  sections  but  tied  together  with  steel  mesh.  This 
mesh  does  not  transfer  a  substantial  part  of  the  load,  but 
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holds  the  slabs  together  so  that  transfer  is  accomplished 
by  aggregate  interlock.  Only  shear  transfer  occurs.  Simi- 
larly, only  shear  transfer  occurs  in  the  keyway  type  of 
joints  . 

The  action  of  repetitive  loading  reduces  the  efficiency  of 
the  joints.  As  load  transfer  is  reduced,  the  slabs 
behave  like  two  separate  entities.  In  doweled  joints,  the 
embedded  extreme  of  the  bar  becomes  loose,  allowing  larger 
deflections.  The  weakening  of  the  concrete-bar  zone  can 
produce  cracking  and  eventual  breaking  of  the  surrounding 
concrete.  In  tied  or  keyway  joints,  the  separation  or  des- 
truction of  the  slab  interface  can  produce  a  decay  in  the 
load  transfer  efficiency. 

3^.3^  Capabilities  and  Limi  tat  ions  of  Available  Methods 
Several  finite  element  computer  methods  for  analysis  of 
rigid  pavements  have  been  developed  in  the  last  twenty 
years.  Herein,  the  following  computer  methods  will  be  dis- 
cused:  1)  ILLISLAB,  2)  JSLAB,  3)  WESLIQID,  and  4) 
WESLAYER. 


In  these  finite  element  methods,  the  pavement  structure  is 
divided  into  three  main  elements:  slab,  sublayers  ,  and 
joints.  All  these  methods  are  based  on  the  fundamental 
expressions  for  thin  plates  developed  by  Zienckiewicz  and 
Cheung.  The  available  methods  of  analysis  can  provide 
stresses  and  deflections  caused  by  a  static  load, 
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representative  of  traffic.  The  designer  chooses  a  cer- 
tain thickness  so  that  for  the  calculated  stress  level, 
the  pavement  is  capable  of  resisting  the  expected  number 
of  load  repetitions.  However,  the  load  considered  in  the 
analysis  is  a  static  load,  in  contrast  to  the  actual  load 
which  is  a  moving  load.  The  response  of  the  pavement  to 
the  actual  load  is  stiffer  than  the  response  obtained  with 
the  assumed  static  load;  the  calculated  stresses  and 
deflections  are  higher  than  the  actual  values.  As  a 
result,  a  conservative,  but  more  expensive  design  is 
obt ai  ne  d . 

The  present  methods  give  results  only  for  the  specific 
properties  of  the  pavement  at  the  moment  considered  in 
the  analysis.  They  cannot  provide  the  stresses  and  deflec- 
tions that  would  occur  after  the  pavement  has  been  in 
service  for  certain  times.  Since  the  state  of  the  pave- 
ment after  several  years  in  service  cannot  be  directly 
determined,  the  maintenance  cost  cannot  be  accurately  cal- 
culated . 


The  mechanical  behavior  of  concrete  has  been  described 
[45]  as  nonlinear  at  low  stresses,  with  expansion  near 
failure.  Nonlinearity  is  caused  by  micr ocracking  due  to 
segregation,  shrinkage,  or  temperature  changes.  The  ini- 
tial cracks  that  exist  in  concrete  develop  and  propagate 
as  loads  are  applied.  Near  failure,  cracks  propagate  and 
join  together,  causing  disruption  in   the   concrete   mass. 
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When  reinforcing  steel  is  provided,  the  cracks  distribute 
in  areas  larger  than  those  under  high  stresses.  Steel 
keeps  the  concrete  areas  together  up  to  stresses  under 
which  disruption  would  otherwise  occur.  The  available 
analysis  methods  assume  a  linear  elastic  behavior  of  con- 
crete at  all  stress  levels.  Prescence  of  reinforcing 
steel  is  not  considered.  Therefore,  the  strains 
corresponding  to  certain  stresses,  especially  high 
stresses,  are  different  from  the  actual  corresponding 
strains  . 

In  laboratory  tests  [45]  the  values  of  strength  and 
modulus  of  elasticity  obtained  for  concrete  depend  upon 
the  rate  with  which  the  load  is  applied.  The  values 
obtained  in  an  impact  test  are  higher  than  those  from  a 
static  slow-rate  test.  In  the  pavement  slab,  loads  are 
more  similar  to  impact  loads  than  to  static  loads.  There- 
fore the  use  of  static  strength  and  static  modulus  of 
elasticity  is  somewhat   unrealistic. 

In  the  sublayers,  the  nonlinear  behavior  is  even  more  pro- 
nounced. Deformations  occur  as  a  result  of  plastic  defor- 
mation of  the  grains  and  aggregations  of  particles.  The 
stiffness  properties  of  the  sublayer  material  are  even 
more  sensitive  to  the  rate  of  load  application.  Therefore, 
the  assumption  of  static  values  with  linear  behavior  is 
farther  away  from  the  actual  behavior. 
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Loads  in  pavements  are  applied  in  a  repetitive  manner. 
This  produces  fatigue  in  the  slab,  the  sublayers,  and  the 
joint  devices.  Fatigue  produces  a  decrease  in  the  pavement 
stiffnesses.  Actually,  a  reduction  in  the  pavement  life  is 
constantly  occurring  due  to  traffic  damage.  Therefore,  in 
the  available  methods  of  analysis,  the  assumption  of  con- 
stant mechanical  properties  for  the  slab,  the  sublayers, 
and  the  joint  devices  is  not  realistic. 

3^.4_  Ideal  Model  for  Analys  is  of  Rigid  Pavement  s 
From  the  design  point  of  view,  the  ideal  analysis  method 
should  be  able  to  compute  the  short  and  long  term 
responses  of  the  pavement.  The  structural  analysis  of  the 
pavement  with  properties  equal  to  those  at  the  beginning 
of  the  service  life  allows  the  design  of  the  geometry  of 
the  pavement.  With  this,  the  initial  cost  can  be 
estimated.  During  service,  the  pavement  deteriorates  and 
its  structural  capacity  decreases.  In  order  to  keep  it  in 
good  condition,  repair  is  performed.  The  cost  of  this 
maintenance  work  depends  on  the  state  of  the  pavement  at 
the  moment  of  repair.  The  analysis  method  should  be  able 
to  predict  the  structural  damage  that  occurs  with  time; 
this  provides  a  guidance  for  the  calculation  of  the 
maintenance  costs. 


From  the  point  of  view  of  the  mechanical  behavior,  the 
ideal  method  of  analysis  should  reproduce  the  actual  per- 
formance of   the   pavement   structure.   The   stress-strain 
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relations  of  the  individual  components  should  be  accu- 
rately represented  in  the  analysis  model.  Concrete  in  the 
pavement  slab  can  be  represented  by  a  nonlinear  plastic 
model.  The  effect  caused  by  the  presence  of  microcracks 
should  be  reflected  in  the  stress-strain  relationships. 
Similarly,  the  presence  of  reinforcing  steel  must  be  taken 
into  account.  In  the  sublayers,  the  material  should  be 
represented  by  a  nonlinear  plastic  model. 

Traffic  loads  have  a  repetitive  nature.  This  causes 
fatigue  in  all  the  stressed  pavement  elements.  The  con- 
crete properties  should  be  dependent  ^upon  the  loading  his- 
tory. Decreases  in  strength  and  modulus  of  elasticity 
should  be  considered;  so  should  be  the  formation  of  cracks 
and  fractures  of  the  slab.  In  the  sublayers,  repetitive 
loading  can  create  areas  of  reduced  support.  Excessive 
water  and  load  repetitions  cause  erosion  and  formation  of 
voids  underneath  the  slab.  This  further  loss  of  support 
should  be  considered  in  the  analysis. 

The  values  considered  for  the  mechanical  properties  should 
be  the  values  obtained  with  loads  similar  in  nature  to 
the  actual  traffic  loads.  Impact  values  for  the  modulus  of 
elasticity  of  concrete  and  for  the  stiffness  of  the  sub- 
layers should  be  used.  These  values  represent  more  accu- 
rately the  performance  of  the  pavement.  Load  transfer 
devices  at  joints  are  also  subjected  to  repetitive  load- 
ing.  Therefore,   the   analysis  method  should  consider  the 
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decay  in  load  transfer  efficiency.  This  decay  can  vary 
from  a  quantitative  variation  of  the  stiffness  values  to  a 
qualitative  variation  of  the  actual  performance.  The 
analysis  method  should  follow  this  kind  of  performance. 
Finally,  the  combination  of  traffic  loading  should  also  be 
included . 

3_.5_  Conclusions 
The  available  finite  element  methods  for  analysis  of  rigid 
pavements  provide  fairly  accurate  results  of  the  stresses 
and  deflections  produced  by  static  load.  With  these 
methods,  any  static-load  distribution  can  be  analyzed  in 
pavements  with  irregular  geometric  characteristics.  How- 
ever, the  analysis  is  carried  out  assuming  linearly  elas- 
tic behavior  for  the  pavement  components.  Moreover,  the 
long-term  effect  of  the  traffic  load  is  not  considered. 
The  presence  of  the  reinforcing  steel  is  not  taken  into 
account  either. 


An  ideal  method  of  analysis  should  reproduce  the  actual 
behavior  of  the  pavement  components.  Nonlinear  stress- 
strain  relationships  are  more  representative  of  the  slab 
and  sublayer  behaviors.  The  ideal  model  should  also  take 
into  account  the  presence  of  reinforcement,  as  well  as  the 
long-term  effect  of  the  traffic  loads.  Like  fatigue  damage 
in  the  concrete  slab  and  pumping  development. 
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CHAPTER  4   NUMERICAL  EVALUATION  OF  EXISTING 
FINITE  ELEMENT  METHODS 

4_.1_  Int  roduc  t  Ion 
In  this  chapter,  results  calculated  using  two  finite 
element  methods  are  compared  with  results  obtained  with 
Wes t ergaa rd '  s  equations,  and  with  the  Pickett  and  Ray 
influence  charts.  Additionally,  deflections  measured  in  a 
laboratory  model  are  used  to  determine  equivalent  compo- 
site values  of  k.  For  the  comparison  of  deflections,  a 
rectangular  plate  on  an  elastic  sublayers  with  full  sup- 
port and  with  interior,  and  edge  loads  is  used.  To 
determine  the  equivalent  composite  values  of  k,  a  small- 
scale  model  of  slab  resting  on  a  soil  sublayers  is  used. 
Numerical  evaluations  of  the  ILLISLAB  and  JSLAB  finite 
element  programs  are  carried  out  by  comparing  the  theoret- 
ical results  obtained  with  these  finite  element  methods 
with  results  obtained  using  closed  form  solutions. 


4_.2_  Cases  S  tudied 
Two  cases  are  considered:  (1)   a  rectangular  concrete  slab 
fully   supported   by   a  uniform  semi-infinite  elastic  sub- 
layers, as  described  in  Figure  4.1;  and  (2)   a  small-scale 
model  consisting  of   micro-concrete  rectangular  slabs 
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INTERIOR  LOAD 


^ -t<LOADED  AREA//pfY; 


EDGE  LOAD 


t8*  254.   mm 
E-  27,579  MPa 
u=»  0.15 

k-  54.2  MPa/m 
P-  44.48  kN 
aa  254.  mm 
b"  254.  mm 


(10  in.) 
(4,000  ksi) 

(200  pel) 
(10,000  lb) 
(10  in.) 
(10  in.) 


Fig.  4.1  Elastic  plate  lying  on  elastic  sublayers. 
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tk~  >;*•.-  *ri  «/*sM--_^ 


slabs;   673  ram  x  406  mm  x  25  ram  (26.5  in.  x  16  in.  x  l.in.) 


Ec=  28,400.  MPa 

u=  0.18 

transverse  joints 

dowel  diameter=   3.5  mm 
dowel  length  =»  51  mm 
dowel  spacing  =  34  mm 
Es=  200,000  MPa 


(4120  ksi) 


(0.15    in.) 
(2.0   in.) 
(1.33in.) 

(29,000  ksi) 


Fig.    4.2    Small-scale    model    of    slab    resting    on    layer 
of    crushed    limestone. 
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resting  on  a  layer  of  crushed  limestone,  as  described  in 
Figure  4.2.  For  the  first  case,  deflections  for  given 
static  loads  are  calculated  using  Westergaard 's  equations, 
Pickett  and  Ray  Influence  Charts,  and  ILLISLAB  and  JSLAB 
finite  element  methods.  For  the  second  case,  deflections 
caused  by  static  and  impact  loads  on  a  laboratory  model, 
are  used  to  determine  values  of  modulus  of  subgrade  reac- 
tion using  a  back-calculation  procedure. 

4.3  Calculation  by  Means  of  Westergaard's  Equations 
The  calculation  of  deflections  in  this  case  is  carried  out 
using  the  following  Westergaard's  Equations: 

For  interior  load: 


\ 


w  =  p/4  \  (3(1:y  >>- 


Et3k 


(4.1) 


-0.275(l-u  )P  ,,  2^,2.,..  2,  2. 
~ —  ((a  +b  )/4+x  +y  ) 

Et 

2 
-0.239(l-w  >p((a2+4ab+b2)/8+(a.b)/(a+b)(-x2+y2) 

Et 


For  ed  ge  load : 


P  ^  2+1  .2u 
\JEt3k 


(l-(0.76+0.4u)b/(JL 


(4.2) 


(l-(0.76+0.4u)y/jL)) 
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Whe  re : 


w=def lection  at  point  (x,y) 


a,b-semiaxes  of  an  ellipse   representing   the   loaded 


area 


x , y=hor izontal  rectangular  coordinates. 

t=thickness  of  slab 

E=modulus  of  elasticity  of  slab. 

y   =  Poisson  ratio  of  slab. 

k=  modulus  of  subgrade  reaction 

J^  =  radius  of  relative  stiffness 

Deflections  are  calculated  assuming  a  circular  area  with 
uniformly  distributed  load.  For  the  cases  shown  in  Fig. 
4.1,  the  calculated  deflections  are: 


max 


(Interior  load)=  0.120  mm  (4.73  mils) 


(Edge  load)     =  0.371  mm  (14.6  mils) 


4.4  Calculation  by  Means  of  Influence  Charts 
Calculation  of   deflections   using   the   Pickett   and   Ray 
Influence   Charts  are  carried  out  with  the  charts  given  in 
reference  [7],   The  charts  developed  for  the  Dense   Liquid 
assumption   were   used.   Scale   adjustment   is  obtained  by 
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determining  the  radius  of  relative  stiffness.   The   number 

of   loaded   blocks  in  the  charts  is  counted  to  the  nearest 

half  block.   For  the  case  shown  in  Fig.  4.1,   the   results 
are  as  follows : 


w      (Interior  load)  =  0.126  mm  (4.99  mils) 
max 


w      (Edge  load)      =  0.360  mm  (14.2  mils) 
max 


4.5  Calculation  by  Means  of  Finite  Element  Methods 
Deflections  for  different  load  conditions  were  calculated 
using  the  JSLAB  and  ILLISLAB  computer  schemes.  For  each 
loading  condition  different  mesh  sizes  and  mesh  configura- 
tions were  assumed.  The  loaded  area  was  assumed  to  be  a 
square  of  25  4  .mm  x  2 5 4. mm  (10.  in.  x  10.  in.),  therefore 
the  applied  pressure  was  689.  kPa  (100  psi)  (Fig.  4.1). 
Two  factors  were  considered  for  both  the  ILLISLAB  and 
JSLAB  computer  schemes:  (1)  effect  of  size  of  slab,  and 
(2)  effect  of  shape  of  elements.  Three  different  slab 
sizes  were  assumed:  5.00  m  x  5.00  m  (197  in.  x  197  in.); 
10.00  m  x  10.00  m  (394  in.  x  394  in.);  and  20.00  m  x  20.00 
m  (787  in.  x  787  in.).  For  the  effect  of  element  shape, 
three  different  side  ratios  were  assumed:  1.0,  0.5,  and 
0.25.  Deflection  results  are  summarized  in  Tables  4.1  and 
4.2,  and  in  Figures  4.3  through  4.6. 
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Table  4.1  Deflections  for  interior  load 


Scheme 


Slab 
Size 
(m) 


5  .00x5  .00 
10  .00x10  .00 
20  .00x20  .00 

5  .00x5  .00 
10  .00x10  .00 
20  .00x20  .00 


Elements   Element 
Size 
(m) 

Effect  of  Slab  Size 


4x4 

8x8 

16x16 

4x4 

8x8 

16x16 


1  .27x1  .27 
1 .27x1  .27 
1 .27x1  .27 

1 .27x1  .27 
1 .27x1  .27 
1  .27x1  .27 


Deflections 
(mm ) 


136 

124 
124 

136 
124 
124 


Effect  of   Element  Shape 


20  .00x20  .00 
20  .00x20  .00 
20  .00x20  .00 

20 .00x20  .00 
20  .00x20.00 
20  .00x20  .00 


16x16 

16x8 

16x4 

16x16 

16x8 

16x4 


1 .27x1  .27 
1 .27x2  .54 
1  .27x5  .00 

1 .27x1  .27 
1 .27x2.52 
1 .27x5  .00 


0.124 
0  .118 
0  .083 

0.124 
0  .118 
0  .084 


I=1LLISLAB 
J=JSLAB 
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Table  4.2  Deflections  for  edge  load 


Scheme      Slab 
Size 
(m) 

Elements   Element    D 
Size 
(m) 

ef lections 
(mm  ) 

ef f  ect 

of  slab  size 

I 
I 
I 

5  .00x5.00 
10  .00x10  .00 
20  .00x20  .00 

4x4 

8x8 

16x16 

1  .27x1  .27 
1 .27x1  .27 
1  .27x1  .27 

0  .392 
0.369 
0.369 

J 
J 
J 

5  .00x5  .00 
10  .00x10  .00 
20  .00x20  .00 

4x8 

8x8 
16x16 

1 .27x1  .27 
1  .27x1  .27 
1 .27x1  .27 

0.392 
0  .369 
0.369 

I 
I 
I 

effect  of 
20.00x20.00   16x16 
20.00x20.00   16x8 
20.00x20.00   16x4 

element  shape 
1 .27x1  .27 
1  .27x2.54 
1  .27x5.00 

0.369 
0.354 
0  .352 

J 
J 

J 

20  .00x20  .00 
20  .00x20  .00 
20  .00x20  .00 

16x16 

16x8 

16x4 

1  .27x1  .27 
1  .27x2.54 
1  .25x5.00 

0  .353 
0  .353 
0.350 

1  = 
J  = 

=ILLISLAB 
=  JSLAB 

64 


150 


ILLISLAB  AND  JSLAB 


00 


SlibSlzi      500x500 
(cm) 

Number  of  4x4 

Elements 


PICKETT  AND  RAY 


1000  xlOOO 
8  x  8 


WFSTERGAARD 


2000x  2000 
16x16 


Fie.  4.3  Deflections  calculated  by  means  of  ILLISLAB    and 
JSLAB  schemes  for  different  slab  sizes  (interior  load). 
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Fig.  A  .A  Deflections  calculated  by  means  of  ILLISLAB  and 
JSLAB  schemes  for  different  element  shapes  (interior  load) 


66 


400 


'o 

3  350 


300 


ILUSLAB  AND  JSLAB 


WESTERGAARD 
□ 


PICKETT  AND  RAy 


Slab  Size    500  x  500 

1  000  x  1  000 

2000  x  2000 

(cm) 

Number  of       4x4 

8x8 

16  x    16 

Elements 

Fig.  4.5  Deflections  calculated  by  means  of  ILLISLAB  and 
JSLAB  schemes  for  different  slab  sizes  (edge  load  ). 
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Fig.  4.6  Deflections  calculated  by  means  of  ILLISLAB  and 
JSLAB  schemes  for  different  element  shapes  (edge  load). 
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4_.6_  Results  from  Labora  tory  Model 
Experimental  deflections  were  obtained  by  testing  a 
laboratory  model  of  a  rigid  pavement  section.  The  1/9 
scale  model  consisted  of  several  slabs  of  673.  mm  x  406. 
mm  x  25.  mm  (26  1/2  in.x  16  in.  x  l.in)  connected  by 
doweled  transverse  joints  and  an  aggregate  interlock  long- 
itudinal joint  (Fig.  4.2).  The  slabs  were  made  out  of 
micro-concrete  of  35,136.  kPa  (5100  psi)  compressive 
strength;  28,400.  MPa  (4120  ksi)  modulus  of  elasticity; 
and  u  =0.18  .  The  slabs  were  connected  at  the 
transverse  joints  by  round-smooth,  steel  bars  3.5  mm  (0.15 
in.)  in  diameter,  51.  mm  (2.  in.)  in  length,  and  34.  mm 
(1.333  in.)  in  spacing. 


Static  deflections  were  obtained  from  static  loads  applied 
gradually  by  means  of  a  mechanical  jack.  For  the  case  of 
impact  load,  deflections  were  registered  while  a  weight 
was  dropped  in  a  small-scale  model  of  falling  weight 
def le ct omet er .  Static  and  impact  deflections  were  measured 
by  means  of  linear  vertical  transducers  connected  to  a 
digital  oscilloscope  and  to  an  oscillographic  recorder. 
Static  loads  were  applied  over  two  rectangular  areas  of  28 
mm  x  17  mm  (1.10  in.  x  0.667  in.)  and  53.6  mm  (2.11  in.) 
at  centers.  Impact  loads  were  applied  over  a  circular  area 
33.8  mm  (1.333  in.)  in  diameter.  In  both  cases,  loading 
areas  were  centered  at  the  middle  of  the  slabs.  With  the 
deflections  obtained  at  each  load  level,  values  of  modulus 


69 

of  subgrade  reaction  were  obtained  by  a  back-calculation 
procedure  using  finite  elements.  This  procedure  consists 
of  determining  the  value  of  k,  such  that  for  a  given 
applied  load,  the  theoretical  deflection  coincides  with 
the  measured  deflection.  An  iteration  scheme  was  imple- 
mented in  the  ILLISLAB  computer  scheme  to  calculate 
automatically  an  average  value  of  k  so  that  the  theoreti- 
cal deflections  are  equal  to  the  measured  deflections. 
Results  are  listed  as  follows: 


Table  4.3  Backcalculat ed  values  of  k  for 
measured  deflections. 


Static  k 

Load 

Def  ] 

.ectlon 

k 

N   (lb) 

ram  ( 
Slab 

>ils) 
No  1 

MPa/m 

(pci) 

329   (74) 

0.050 

(1.99) 

34 

(126) 

507  (114) 

0  .075 

(2.96) 

36 

(133) 

658  (148) 

0.099 

(3.91) 

34 

(129) 

822  (185) 

0  .120 

(4.73) 

37 

(136) 

987  (222) 

0.134 

(5.30) 

40 

(151) 

1152(259) 

0  .157 
Slab 

(6.20) 
No  2 

40 

(150) 

329   (74) 

0  .083 

(3.27) 

16 

(  60) 

507  (114) 

0.104 

(4.11) 

22 

(  81) 

658  (148) 

0  .130 

(5.15) 

23 

(  86) 

822  (185) 

0  .146 

(6.23) 

27 

(101) 

987  (222) 

0  .158 

(6.23) 

32 

(118) 

1152(259) 

0.176 
Slab 

(6.94) 
No  3 

34 

(127) 

329   (74) 

0.036 

(1.44) 

55 

(206) 

507  (114) 

0  .049 

(1.95) 

67 

(249) 

658  (148) 

0.061 

(2.42) 

72 

(267) 

822  (185) 

0  .071 

(2.78) 

82 

(303) 

987  (222) 

0  .077 

(3.03) 

95 

(351) 

1  152(259) 

0  .084 

(3.31) 

105 

(388) 

Table  4.4  Backcalculat ed  values  of  k  for 
measured  deflections. 
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Load 
N    (lb) 


298  (67) 

618  (139) 

916  (206) 

1205  (271) 


293  (66) 

631  (142) 

907  (204) 

1187  (267) 


2  7  1  (61) 

587  (132) 

929  (209) 

1236  (278) 


Impact  k 

Def lection 
mm    (mils ) 

Slab  No  1 
0.016  (0.642) 
0.034  (1.357) 
0.052  (2.071) 
0.071  (2.805) 

Slab  No  2 

0.023  (0.930) 

0.048  (1.903) 

0.069  (2.723) 

0.096  (3.788) 

Slab  No  3 

0.020  (0.786) 

0.042  (1.664) 

0.063  (2.513) 

0.088  (3.467) 


MPa/m  (pci) 


170  (629) 

164  (606) 

157  (581) 

149  (553) 


93  (344) 

101  (373) 

101  (375) 

92  (340) 


105  (391) 

111  (410) 
119  (440) 

112  (415) 


The  average  values  of  the  static  modulus  of  subgrade  reac- 
tion were  obtained  by  means  of  the  least  square  error  pro- 
cedure. The  static  k  for  a  load  of  750.  N  (168  lb)  was  47. 
MPa/m  (173  pci);  while  the  impact  value  of  k  for  the  same 
load  was  122.8  MPa/m  (453  pci).  The  difference  of  the 
averages  is  75.8  MPa/m  (280  pci).  The  least  square  aver- 
ages are  shown  in  Fig.  4.7. 
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Fig.  4.7  Modulus  of  subgrade  reaction  obtained  from  impact 
and  static  deflections  measured  in  laboratory  model. 
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4_.7_  Summary 
Deflections  in  a  square  plate  supported  by  an  elastic 
sublayers  caused  by  interior  and  edge  loads  were  calcu- 
lated using  tbe  ILLISLAB  and  JSLAB  computer  schemes.  The 
results  obtained  by  means  of  these  two  finite  element 
methods  were  compared  with  results  obtained  by  Wester- 
gaard  equations  and  the  Pickett  and  Ray  influence 
charts.  Different  slab  sizes,  number  of  elements,  and 
shape  of  elements  were  used.  For  both  the  interior  and 
edge  loads  deflections  calculated  by  means  of  the  finite 
element  methods  were  closer  to  the  Westergaard  and  Pick- 
ett and  Ray  solutions,  when  larger  slab  sizes  or  more 
elements  were  used.  For  different  element  shapes,  the  two 
finite  element  schemes  give  the  same  results  for  the  case 
of  interior  load.  A  specific  trend  is  not  defined  however. 
For  the  case  of  edge  load,  the  JSLAB  scheme  gives  more 
consistent  results  for  all  element  shapes  than  ILLISLAB 
does  . 


As  reported  by  Thompson  et  al  [44],  the  backca lculat ed 
values  of  k  are  higher  for  an  Impact  load  than  for  a 
static  load.  For  the  static  load,  the  k  values  appear  to 
increase  as  the  load  increases.  This  may  be  due  to  an 
increase  in  the  composite  value  of  k  resulting  from 
adjustments  of  the  slab-sublayer  interface  and  to  an 
increase  of  the  area  of  contact  as  load  Increases.  Initial 
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support  may  be  low  due  to  curling.  The  value  of  k  at  dif- 
ferent load  levels  would  be  expected  to  decrease  as  the 
magnitude  of  the  load  increases.  This  trend  is  not 
observed  very  markedly  in  Fig.  4.7,  but  the  value  of  k 
decreases  as  the  magnitude  of  the  loai  increases. 
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CHAPTER  5   COMPUTER  METHOD  FOR  ANALYSIS  OF  RIGID  PAVEMENTS 

5.«J_  Int  roduc  t  ion 

In  this  chapter,  a  Nonlinear  Finite  Element  Method  for 
analysis  of  rigid  pavements  with  fatigue  is  presented.  It 
is  termed  the  Purdue  Method  for  Analysis  of  Rigid  Pave- 
ments (PMARP).  Fatigue  damage  caused  by  traffic 
deteriorates  the  pavement  components  so  that  their  stiff- 
ness properties  decay  with  time.  The  modulus  of  elasticity 
and  strength  of  concrete  decrease  as  traffic  imposes  load 
repetitions  stressing  and  straining  the  concrete.  Micro- 
cracking  takes  place  reducing  the  strength  and  modulus  of 
elasticity  of  the  concrete  slab.  Voids  are  also  produced 
due  to  pumping  action  created  by  repetitive  traffic  loads. 
The  method  of  analysis  proposed  herein  takes  into  account 
the  damaging  effect  caused  by  load  repetitions  on  the  slab 
as  well  as  on  the  sublayers.  The  damage  is  quantified  as 
decay  in  the  stiffness  of  concrete,  amount  of  cracking, 
decay  in  load  transfer  efficiency,  and  amount  of  damage 
due  to  pumping. 


The  computer   implementation   is   based   on   the   computer 
method   ILLISLAB  previously  developed  at  the  University  of 
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Illinos  [9].  In  the  method  presented  herein,  the  pavement 
components  are  assumed  to  have  a  stress-  or  strain- 
dependent  behavior,  and  the  concrete  slab  is  considered  as 
an  orthotropoic  plate.  Methods  for  including  damage  in 
the  concrete  slab,  damage  in  the  load  transfer  devices, 
and  formation  of  voids  are  described.  In  the  first  part, 
an  outline  of  the  behavior  of  rigid  pavements  under 
repetitive  loading  is  presented.  A  method  for  analysis  of 
rigid  pavements  with  fatigue  is  given  at  the  end  of  the 
chap  ter . 


5_.2_  S  t  rue  tural  Behavior  of  Rigid  Pavements  with  Pumping 
The  pavement  structure,  consisting  of  a  series  of  layers 
made  up  of  different  materials  with  different  stiffnesses, 
provides  a  smooth  surface  to  ride  on,  and  a  way  to 
transfer  the  highly  concentrated  loads  to  the  lower  layers 
at  successively  lower  stresses.  The  concrete  slab 
transfers  the  vertical  stresses  into  a  larger  area  in  the 
subbase  which,  in  turn,  transfers  them  to  even  a  larger 
area  in  the  subgrade  (Fig.  5.1).  The  slab  is  under  the 
action  of  bearing  and,  more  important,  flexural  stresses, 
while  the  subbase  and  subgrade  are  under  vertical  and  con- 
finment  stresses.  The  level  of  stresses  occurring  at  each 
pavement  component  depends  on  the  relative  values  of 
stiffnesses  of  the  individual  structural  elements.  For 
example,  a  relatively  stiff  slab  (  high  concrete  strength 
or  large  thickness)  supports  the  majority  of  the  load, 
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Figure  5.1  Idealized  vertical  cut  of  pavement  structure 
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transferring  relatively  low  stresses  to  the  subbase  and 
subgrade  (the  load  is  distributed  over  a  large  area  in  the 
subbase  and  subgrade). 

In  the  conventional  analysis  of  rigid  pavements,  the  vari- 
ation of  the  pavement  properties  caused  by  repetitive 
application  of  traffic  loads  is  not  taken  into  account. 
The  traditional  premise  underlying  the  design  of  rigid 
pavements  is  that  full  contact  exists  between  the  slab  and 
the  supporting  materials.  However,  several  factors  can 
make  the  concrete  slab  lose  contact  with  the  subbase 
including  warping,  curling,  and  pumping.  In  addition, 
fatigue  damage  caused  by  a  repetitive  load  application 
reduces  the  strength  and  stiffness  of  concrete  pavements. 
Therefore,  the  response  of  a  pavement  to  a  certain  load 
application  depends  on  its  structural  soundness,  which 
varies  with  the  service  life  of  the  pavement. 


Pumping  can  be  analyzed  from  several  viewpoints.  Herein, 
only  the  structural  aspect  of  pumping  is  considered.  The 
main  effects  produced  by  pumping,  from  the  structural 
viewpoint,  are  the  development  of  zones  which  create  low 
bearing  capacity,  or  even  voids.  In  the  traditional 
methods  of  analysis,  the  pavement  structure  is  assumed  to 
be  in  a  full  and  perfect  contact  at  all  points  of  the 
slab-sublayers  interface.  However,  this  condition,  if  it 
ever  occurs,  exists  only  at  the  beginning  of  the  service 
life  of  a  pavement.  After  a  pavement  is  opened  to  service, 
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the  environment  and  traffic  actions  cause  changes  in  its 
structural  conditions. 

The  effect  of  pumping  is  the  ejection  of  water  with  sub- 
layer material  from  under  the  slab  producing  a  void.  Free 
water  can  accumulate  in  an  initial  void,  from  where  it  can 
be  ejected  through  cracks  and  joints  when  the  slab  is 
deflected  down  by  a  passing  load.  Muddy  water  is  ejected 
indicating  erosion  of  the  subbase.  The  initial  voids  or 
spaces  underneath  the  slab  can  be  produced  by  temperature 
curling  of  the  slab  or  by  plastic  deformation  of  the  sub- 
base.  The  continuing  application  of  loads  increases  the 
size  of  the  initial  voids  until  a  relatively  large  space 
is  formed.  The  change  in  the  support  conditions  when 
voids  are  formed  can  augment  the  critical  stresses  in  the 
slab.  Cracking  and  faulting  may  occur  as  a  result  of  this 
increase  of  stresses. 

In  the  slab,  the  strength  of  the  concrete  increases  with 
time  as  the  hydration  process  develops.  At  the  same  time, 
however,  seasonal  and  daily  temperature  changes  produce 
differential  dilatations  in  the  components  of  the  con- 
crete which  causes  micr ocracking .  This  mi cr ocracking 
reduces  the  strength  of  concrete.  It  can  be  assumed  that 
these  two  effects  can,  approximately,  cancel  each  other. 


Additionally,  pavement  structures   are   subjected   to   the 
action  of  traffic.   The  repetitive  application  of   traffic 
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loads  causes  stressing  and  fatigue  damage  in  the  concrete 
slab.  This  is  the  type  of  damage  considered  in  the  method 
presented  herein. 

5_  .3_  Proposed  Method  for  Ana  lysis  of  Rigid  Pavement  s 
In  the  following,  a  Non-linear  Finite  Element  method  for 
analysis  of  rigid  pavements  (PMARP)  considering  the  decay 
of  the  pavement  stiffness  properties  due  to  traffic  damage 
is  presented.  This  method  can  be  used  to  calculate  the 
stresses  and  deflections  of  rigid  pavements  subjected  to 
fatigue  damage  caused  by  traffic.  It  can  also  be  used  to 
estimate  the  current  state  of  damage  of  rigid  pavements 
that  have  been  under  traffic  for  certain  time.  Thus,  the 
method  can  be  applied  in  design  as  well  as  in  maintenance. 


In  the  first  part  of  this  section,  the  bases  for  the 
development  of  the  nonlinear  method  for  structural 
analysis  of  rigid  pavements  are  presented.  The  pavement 
structure  is  divided  into  several  well  defined  components 
whose  behaviors  reflect  the  behavior  of  the  whole  pave- 
ment structure.  Since  the  behaviors  of  these  pavement 
components  are  stress-  or  strain-  dependent,  the  solution 
is  obtained  through  an  iterative  process.  In  the  second 
part,  the  computer  implementation  of  the  method  is 
illustrated. 
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5.3.1  Finite  Element  Discretization 
The  method  presented  herein  is  based  on  a  finite  element 
technique  in  which  the  whole  structure  is  divided  into 
specific  discrete  elements.  Two  different  methods  for 
dividing  the  pavement  structure  have  been  generally  used. 
In  the  first  method,  the  pavement  is  divided  into  three 
basic  elements:  sublayer  elements,  load  transfer  elements, 
and  slab  elements.  In  the  first  case,  the  slab  is  formed 
by  a  combination  of  beam  and  torsion  elements.  (i.e.  Hud- 
son and  Matlock)  [34]  In  the  second  technique,  the  pave- 
ment structure  is  also  divided  into  three  elements,  but 
the  slab  is  divided  into  plate  elements  [33,38]  .  In  the 
method  presented  herein,  the  second  approach  is  used;  the 
slab  elements  are  treated  as  thin  plate  elements.  See  Fig- 
ure 5.2. 

Slab  Elements .-  The  concrete  slab  is  assumed  to  behave 
like  a  thin  plate  with  small  deformations.  Kirchoff 
hypotheses  are  applicable: 

( i ) .-Def lect ions  at  midsurface   are   small   com- 
pared with  thickness  of  plate. 

(ii).-The   raidplane   remains   unstrained   during 
bending. 


(iii). -Plane  sections  initially  normal  to  the 
midsurface  remain  plane  and  normal  to  that  sur- 
face, after  bending. 
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Figure    5.2    Decomposition    of    pavement    section 
finite    elements  . 
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(iv).-The  stress  normal  to  the  midplane,  a^  , 
is  small  compared  with  other  stress  components 
and  therefore  can   be  neglected. 

From  hypothesis  (i),  the  slope  of  the  deflected  surface  is 
assumed  to  be  small  and  the  square  of  the  slope  can  be 
neglected  when  comparing  with  unity.  Hypothesis  (iii) 
implies   that  the  vertical  shear  strains  ,   Txz   and   Yyz 

are  negligible.  The  deflection  of  the  plate  is  therefore 
associated  principally  with  bending  strains.  The  normal 
strain,  e  ,  resulting  from  transverse  loading  may  also 
be  omitted.  Hypothesis  (iv)  is  acceptable  if  transverse 
loads  are  not  highly  concentrated. 

If  an  overlay  or  stabilized  subbase  are  present,  such 
layers  are  also  considered  as  thin-plate  elements.  No 
bond  is  assumed  at  the  interfaces,  thus,  the  stiffness  of 
the  whole  group  of  layers  is  that  of  three  superimposed 
plates.  No  composite  action  is  considered. 

The  finite  element  formulation  of  the  plate  elements  is  as 
given  by  Zienkiewicz  and  Cheung  [33,38].  In  this  formula- 
tion, three  displacements  are  considered  as  nodal  parame- 
ters:  transverse   displacements,    w    ;  rotation  about  x 


axis , 


Xi 


;  and  rotation  about  Y  axis, 


Yi 


,        as       shown 


in    Fig.     5.3.    Each    element    has     twelve    degrees     of     freedom. 
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Figure    5.3    Isotropic    and    orthotropic    thin    plate    elements 
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Corresponding  with  these  displacements,  there   are   three 

nodal   forces:   transverse   force,    f  ;  a  couple  about  x 

axis,   f    ;  and  a  couple  about  y  axis,  f    .  Therefore: 
nodal  displacements 


N 


Xi 
no  dal  forces 


(5.1) 


f         N 
fi 


f,  ~h 
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For  an  element   e  ,  the  nodal  forces  and  displacements  are 
re lated  by: 

ke  ae   =   fe  (5.3) 

For  the  entire  slab  composed  of  n  plate  elements,  we  have: 


[K] 
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(5.4) 


where     [K]     is     the    stiffness    matrix    of     the    whole    slab    ; 


i  =  n 

Z[kf 
1-1 


(5.5) 


{a}  and   {f }  are  the  nodal  displacements  and  nodal   forces 
for  the  whole  structure. 


In  the  present  method,  there  are  two  options  to  represent 
the  concrete  slab:  (1)  as  an  isotropic  slab  ,  and  (2)  as 
an  orthotropic  slab.   The  orthotropic  case  is  to  take  into 


85 

account   the  presence  of  reinforcement  in  both  the  X   and 

Y   directions.  For  the  orthotropic  case,  only  one  plate 

layer   can   be  considered.   The  stiffness  matrices  for  the 
orthotropic  and  isotropic  cases  are  as  follows: 

ke  =  //    BtDBdxdy  (5.6) 

The  result  of  this  integration   is   given   by   Zienkiewicz 
[46]  (p  237), 


ke  =  rn1l     (D  K.+D  K-+D.K.+D   K.  )L 


(5.7) 


60ab  v  x  1  y  2  13  xy  4 
where  K  ,  K_  ,  K_  ,and  K,  have  constant  values  depend- 
ing on  the  dimensions  ,  a  and  b,  of  the  element.  L  is 
also  a  constant  [46]. 


The  terms   D   ,  D   ,  D.  ,  and  D     represent  the   flexural 
x     y     1         xy 

and   torsional  rigidities  of  the  orthotropic  plate  and  are 
calculated  as  follows: 
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and   E    ,  and 
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are  the  moduli   of 


elasticity  and  Poisson  ratios  in  the  X  and  Y  directions. 
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Stresses  are  computed  by: 

E 

o       =   — : (e   +u    e    ) 

x  1"uxwy  x       x    y 

E 

o      =   — -. ^ (e    +y    e    ) 

y  1-wxUy         y      y    x 


(5.9) 
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Where     the    strains     depend    on    the    nodal       displacements,       the 
displacement    function,     and    the    nodal    coordinates; 


ove,   ja| 


(5.10) 


In  the  above,   <a>   is  obtained  from  the  solution   of   the 

F.E.  equations  and  [B]=  QC  .  Q  depends  on  the  dis- 
placement function  and  C  on  the  nodal  coordinates. 
Expressions   for   Q  and  C  can  be  found  in  reference  [46]  p 

236 ,  and   y   ■  y   »  y 
'        x     y 


Foundation  Elements .-  The  subgrade  is  considered  as  a 
Winkler  sublayers  in  which  the  deflection  at  a  point  is 
assumed  to  depend  only  on  the  pressure  acting  at  that 
point.  Therefore,  the  subgrade  can  be  represented  by  a 
series  of  discrete  vertical  elements  with  only  one  degree 
of  freedom. 

In  the  method  presented  herein,  the  vertical  elements  are 
assumed  to  have  a  nonlinear,  deflection-dependent 
behavior.   The   values   of   the   coefficient   of   subgrade 


reaction  are  resilient  values  obtained  from  impact  tests 
[44]  .  These  coefficients  are  the  ratio  between  the 
applied  pressure  and  the  resulting  deflection.  The  verti- 
cal elements  are  placed  at  the  nodal  points  of  the  slab. 
The  displacement  component  at  each  nodal  point  is  a  verti- 
cal displacement,  w  ,  in  the  Z  direction.  The 
corresponding  force  component  is  the  vertical  force  f  . 
The  force  and  displacement  are  related  by: 
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Load  Trans  f er  Element  s  .-  The  concrete  slab  is  divided  into 
small  sections  by  longitudinal  and  transverse  joints.  At 
these  joints,  devices  to  transfer  portions  of  the  load 
among  adjacent  slabs  are  provided.  In  the  present  method, 
joints  can  be  considered  as  doweled  ,  aggregate  interlock, 
and  keyed  joints  [9]  (Fig.  5.4).  Doweled  joints  are 
represented  by  bar  elements  which  are  capable  of  transfer- 
ring shear  and  moment  forces  across  the  opening  of  the 
joint.  These  elements  have  two  degrees  of  freedom:  a 
vertical  displacement 
izontal  axis  transverse  to  the  bar   8 


w   ,   and  a  rotation  about  a  hor- 


Yi 


Corr esp  onding , 


there  are  two  nodal  forces:  a  shear  force   f .  ,  and  a 
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Figure    5.4    Finite    element    representation    of    pavement     joints 


89 


moment  about  a  horizontal  axis  f„.  .  The  displacements 
and  forces  are  related  by  the  stiffness  matrix  of  the  bar. 
Aggregate  interlock  and  keyed  joints  are  represented  by 
vertical  elements  with  only  one  degree  of  freedom.  These 
elements  have  a  vertical  displacement  and  a  correspondent 
vertical  force  (Fig.  5.4).  The  value  of  the  stiffness  of 
this  element  depends  upon  the  type  of  joint  and  the  par- 
ticular conditions  of  the  joint. 


5.3.2  Fatigue  Damage  in  Concrete  Pavements 
Highway  pavements  are  subjected  to  many  repetitions  of 
traffic  loads  during  their  service  lives.  The  importance 
of  fatigue  damage  in  the  performance  and  design  of  rigid 
pavements  is  self-evident.  Each  time  a  vehicle  travels 
over  a  road  its  weight  is  subsequently  supported  by  dif- 
ferent zones  of  the  road.  For  certain  critical  positions 
of  the  load,  its  effects,  like  stresses  and  deflections, 
reach  maximum  values.  Generally,  these  maximum,  or  criti- 
cal, values  are  used  in  the  design  of  the  entire  pavement. 
Since  the  traffic  consists  of  vehicles  of  different  con- 
figurations and  sizes,  the  magnitude  of  the  traffic  load 
depends  on  the  particular  vehicle  traveling  over  the  road 
at  a  given  time.  Heavy  trucks  impose  heavy  loads  with 
corresponding  high  stresses.  A  complete  representation  of 
loads  imposed  on  a  road  can  be  obtained  from  the  traffic 
serviced  by  the  road.  This  traffic  is  quantified  by  the 
number  of  occurrences  of  vehicles  of  each   different   type 
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and  size.  Since  it  would  be  very  complex  to  analyze  a  par- 
ticular pavement  considering  the  loading  of  each  vehicle, 
it  is  necessary  to  introduce  some  simplifications.  First, 
it  is  assumed  that  the  total  traffic  can  be  classified  in 
a  finite  number  of  vehicle  types.  Second,  the  Palraegran- 
Miner  hypothesis  is  assumed  to  be  valid,  that  is:  the  dam- 
age caused  by  several  vehicles  of  different  types  is  cumu- 
lative and  independent  of  the  order  in  which  the  vehicles 
travel  over  the  pavement  [47].  The  supporting  of  loads  is 
carried  out  by  all  the  interconnected  elements  of  the 
structure:  the  concrete  slab,  the  subbase,  and  subgrade 
together  resisting  the  load.  When  a  vehicle  passes  over 
the  pavement  all  the  components  of  the  structure  suffer 
certain  fatigue  damage.  The  fatigue  damage  suffered  by 
each  structural  element  depends  upon  the  relative  magni- 
tude of  stress  or  strain  occurring  in  that  element,  and  on 
the  fatigue  properties  of  the  particular  element.  In  the 
following  ,  the  fatigue  behavior  of  the  pavement  com- 
ponents is  treated  from  the  viewpoint  of  a  pavement  struc- 
ture . 


Fa t  igue  in  Slab .-  The  concrete  slab  constitutes  a  very 
important  structural  element  in  rigid  pavements.  For  the 
concrete  slab,  as  for  any  other  concrete  structure  sub- 
jected to  repetitive  loading,  it  is  important  to  consider 
the  possibility  of  fatigue  failure.  The  fact  that  concrete 
can   suffer   from  fatigue  under  repeated  stresses  has  been 
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recognized  since  the  early  1900's  .  Since  then,  a  great 
deal  of  research  has  been  conducted  on  the  behavior  of 
concrete  under  fatigue  loading.  Yet,  at  the  present  time, 
a  perfect  model  of  the  fatigue  behavior  of  concrete  is  not 
avai lable . 


In  rigid  pavements  the  concrete  slab  is  constantly  sub- 
jected to  stresses  of  different  natures.  Traffic  loads 
produce  bending  stresses  as  well  as  impact  and  horizontal 
stresses.  Environmental  actions,  such  as  temperature  and 
moisture  changes,  are  also  continuously  affecting  the  con- 
crete of  the  slab.  All  these  actions  produce  repetitive 
stresses  and  strains  in  the  structure  of  the  concrete 
causing  damage  and  , eventually,  failure.  In  the  method  of 
analysis  considered  herein,  only  the  damage  produced  by 
repetitive  loading  is  taken  into  account. 

Fatigue  performance  is  usually  expressed  in  terms  of  an 
endurance  curve  in  which  the  mean  value  of  the  number  of 
cycles  to  failure  under  a  particular  loading  condition  is 
plotted  together  with  the  magnitudes  of  the  cyclic  stress 
applied  normalized  with  respect  to  the  flexural  strength 
(Fig.  5.5).  Various  forms  of  fatigue  tests  have  been  used 
to  investigate  the  fatigue  behavior  of  plain  concrete,  for 
example:  direct  compression,  direct  tension,  indirect  ten- 
sion, and  flexure.  For  rigid  pavements,  the  fatigue 
behavior  in  flexure  is  of  particular  interest. 
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Figure    5.5    Typical    flexural    fatigue    performance    of    concrete 
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Several  factors  affect  the  fatigue  performance  of  concrete 
in  flexure:  stress-history,  loading  amplitude,  rate  of 
loading,  rest  periods,  age  of  testing,  and  moisture  condi- 
tions. In  the  following  ,  a  brief  discussion  of  these 
factors  is  presented. 

In  practice,  service  loading  is  not  applied  in  regular 
cyclic  patterns  over  periods  of  hours  or  days  as  in  a  con- 
ventional fatigue  test.  Traffic  loading  is  constantly 
varying  in  amplitude  and  frequency.  There  are  few  pub- 
lished results  of  the  behavior  of  concrete  under  variable 
amplitude  loading,  in  spite  of  the  fact  that  service  load- 
ing histories  on  any  concrete  structure  are  certainly  not 
represented  by  constant  amplitude  loadings.  In  a  real 
structure,  under  a  given  loading  the  stress  distribution 
may  change  with  time  as  the  modulus  of  elasticity  and 
strength  vary  with  age  and  with  moisture  movements.  The 
only  true  form  of  representative  loading  would  be  to 
reproduce  service  loading  histories  in  real  time  together 
with  realistic  moisture  movements  and  temperature  changes. 
This,  however,  would  be  impracticable  and  would  not  con- 
stitute an  accelerated-test  procedure.  In  view  of  this,  a 
simple  constant  amplitude  type  of  loading  Is  used  and  a 
simple  cumulative  fatigue  damage  is  assumed  for  the  esti- 
mation of  the  fatigue  behavior  of  concrete  [48]  . 


Although  it  has  been   known   for   a   long   time   that   the 
compressive   and   flexural   strengths  of  concrete  increase 
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with  the  increase  in  rate  of  loading,  the  effect  on 
fatigue  performance  is  less  marked  [48]  .  In  the  same 
manner,  rest  periods  between  loading  cycles  seem  unlikely 
to  have  a  very  significant  effect  of  fatigue  performance 
of  concrete  [49], 

For  a  given  stress  level,  there  is  a  marked  increase  in 
cycles  to  failure  as  the  age  of  concrete  increases  [49]. 
At  the  same  time  there  is  a  corresponding  increase  in  the 
static  strength  so  that  the  applied  fatigue  loading 
represents  a  diminishing  proportion  of  the  static  strength 
at  the  time  of  testing.  If,  however,  all  the  cycles  to 
failure  are  plotted  against  the  appropriate  flexural 
strengths  at  different  ages,  all  the  points  tend  to  lie 
close  to  a  single  endurance  curve.  However,  the  changes  in 
temperature  and  moisture  produce  differential  stresses  in 
the  concrete  mass.  These  differential  stresses  may  produce 
microcracking  in  the  pavements  and  thus  the  strength  and 
modulus  of  elasticity  of  concrete  may  be  reduced.  For 
simplicity,  the  gain  in  strength  after  28  days  is  assumed 
to  be  reduced  by  microcracking  caused  by  temperature  and 
moisture  changes.  Therefore,  the  variation  of  strength  of 
concrete  due  to  maturity  can  be  neglected  if  no  well 
defined  variation  of  this  strength  is  known. 


Decay  in  Load  Transfer  Efficiency .-  The  occurrence  of 
pumping  is  greatly  affected  by  the  efficiency  of  the  load 
transfer  devices.  Unfortunately,  this   is   an   area   where 
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considerable  research  remains  to  be  done  in  order  to  have 
a  clear  understanding  of  the  phenomena  occurring  at  the 
joints.  What  is  clear  now  is  that  through  the  service 
life  of  a  pavement  the  load  transfer  efficiency  is  modi- 
fied by  the  action  of  traffic.  The  application  of 
extended  repetitive  loading  decreases  the  initial  ability 
of  a  given  system  to  transfer  load. 


From  the  little  research  that  has  been  conducted, 
[21,22,23,24,25]  some  conclusions  can  be  drawn:  the  abil- 
ity to  transfer  loads  depends  on  many  variables,  such  as 
joint  width,  dowel  diameter,  dowel  length,  dowel  spacing, 
initial  dowel  looseness,  and  concrete  strength  and  modulus 
of  elasticity.  Teller  and  Cashell  [28]  found  that  a  defin- 
ite exponential  relation  exists  between  dowel  diameter  and 
load-transfer  capacity,  other  variables  being  constant. 
For  a  given  dowel  diameter  and  condition  of  loading, 
decreasing  the  width  of  the  joint  opening  decreases  the 
bending  stresses  in  the  dowel.  It  also  decreases  the  dowel 
deflections  and  hence  increases  the  percentage  of  load 
transfer  both  initially  and  after  extended  repetitive 
loading.  The  condition  of  dowel  looseness  has  an  important 
effect  on  the  structural  performance  of  the  dowel,  since 
it  can  function  at  full  efficiency  only  after  this  loose- 
ness is  taken  up  by  load  induced  deflection.  This  is  true 
for  both  the  initial  looseness  and  the  subsequent  loose- 
ness developed  during   the  course  of  repetitive    loading 
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Mainly,  the  decrease  in  the  load  transfer  efficiency  is 
due  to  damage  that  occurs  in  the  concrete  surrounding  the 
dowel  bars.  As  reported  in  reference  [28],  loss  in  ini- 
tial capacity  to  transfer  load  depends  on  the  number 
of  load  cycles  and  on  the  magnitude  of  the  load  (Fig. 
5.6).  From  the  data  reported  by  Cashell  and  Teller  [28], 
a  regression  equation  can  be  found  to  relate  the  values  of 
modulus  of  dowel  support  with  the  number  of  load  cycles 
and  the  magnitude  of  the  stress  at  the  bar-concrete  inter- 
face. 


5.3.3  Pumping  Development 
From  the  structural  point  of  view,  pumping  is  important 
whenever  voids  are  formed  under  the  slab  and  the  support 
conditions  of  the  pavement  are  modified.  In  the  method  of 
analysis  presented  in  this  chapter,  formation  of  voids 
resulting  from  pumping  is  included.  Support  conditions  in 
the  pavement  are  modified  in  accordance  with  the  magnitude 
of  the  voids  formed  during  pumping  action.  Estimation  of 
void  sizes  is  achieved  by  means  of  an  analytical  method 
which  has  been  developed  in  part  from  statistical 
analysis  of  the  pumping  indices  reported  in  the  AASHO  Road 
Test  [2],  and  in  part  from  the  structural  analysis  of  the 
pavement  sections  used  in  that  test.  Pumping  index  is 
defined  as  the  volume  of  material  pumped  out  from  under 
the  slab,  per  unit  length  of  pavement. 
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Figure  5.6  Increase  in  dowel  looseness  as  a 
function  of  load  cycles. 
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The  material  pumped  out  from  underneath  the  concrete  slab 
leaves  a  space  which  modifies  the  support  conditions  of 
the  slab.  The  voids  formed  during  the  pumping  action  offer 
an  area  of  low  or  no  support  for  the  slab.  The  identifi- 
cation of  the  shape  and  location  of  the  voids  allows  for  a 
more  accurate  analysis  of  the  pavement  structure.  In  the 
following  paragraphs,  a  model  to  define  the  development  of 
voids  in  concrete  pavements  is  presented. 

The  amount  of  material  pumped  in  a  rigid  pavement  depends 
on  the  following  factors: 

1  .-Structural  properties  of  pavement. 

2. -Pumping  susceptibility  of  subbase  material. 

3. -Magnitude  and  number  of  load  applications 

4. -Climatic  conditions. 

The  analytical  model  presented  herein  consists  of  two 
parts:  (1)  a  qualitative  description  of  location  and  shape 
of  voids,  and  (2)  a  quantitative  description  of  the  size 
of  the  voids. 


Qualitative  Description .  The  qualitative  description  of 
the  location  and  shape  of  voids  is  based  on  field  observa- 
tions and  on  analysis  of  deflection  patterns  of  rigid 
pavements.  It  has  been  observed  that  pumping  occurs  at  the 
edges  of  the  pavement  slabs,   and   at   the   sides   of   the 
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transverse  joints,  with  the  most  significant  pumping 
occurring  at  the  leave  side  of  the  joint  [50,1], 

Pumping  occurs  when  water  is  ejected  at  high  speed  from 
underneath  the  slab  when  it  is  deflected  downward  by  a 
passing  load.  Therefore,  ejection  of  water  takes  place 
where  free  water  accumulates  and  where  high  deflections 
occur.  Water  accumulates  at  those  places  where  initial 
space  exists  between  the  slab  and  the  supporting  material. 
The  initial  spaces,  or  voids,  can  be  caused  by  negative 
temperature  gradient  in  the  slab  (the  top  of  the  slab  has 
a  lower  temperature  with  respect  to  the  bottom  of  the 
slab),  or  by  plastic  deformation  of  the  sublayers.  From 
the  structural  analysis  of  pavement  sections  under  nega- 
tive temperature  differential,  it  can  be  inferred  that  the 
slabs  adopt  a  bowl  shape,  leaving  spaces  at  the  corners  of 
the  slabs.  At  the  corners  also,  specially  at  the  exterior 
corner,  the  maximum  deflections  take  place. 


For  the  pumping  model  described  here,  it  is  assumed  that 
pumping  occurs  mainly  at  two  places:  at  the  joints,  and 
along  the  exterior  edge  of  the  pavements.  It  is  assumed 
that  voids  start  to  develop  at  the  exterior  corner  of  the 
slab  and  propagate  towards  the  centerline  of  the  pavement. 
The  shapes  of  the  areas  of  the  voids  are  assumed  to  be 
rectangular  at  the  edge  and  triangular  at  the  corner.  For 
purposes  of  structural  analysis  and  for  allocation  of 
voids  with  no  support  in  the  finite  element   system,   void 
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areas  will  be  provided  to  the  nodes  with  the  highest 
deflections,  as  described  later  in  the  computer  implemen- 
tation. 

Quantitative  Description.  A  mathematical  model  to  predict 
Che  amount  of  pumping  has  been  developed  from  the  AASHO 
Road  Test  data.  This  model  considers  only  the  structural 
features  of  the  pavement  and  the  amount  of  traffic.  With 
this  model  a  so  called  "pumping  potential"  can  be  deter- 
mined for  a  certain  pavement  section  and  volume  of  traffic 
given  in  ESAL  .  The  pumping  potential  is  the  pumping  index 
that  would  occur  for  the  conditions  which  were  present  in 
the  AASHO  Road  Test;  mainly,  a  non-treated  subbase  and  the 
rainfall  conditions  of  that  location.  The  pumping  poten- 
tial can  be  modified  by  multiplying  it  by  factors  which 
depend  on  the  erosion  susceptibility  of  the  subbase  (  type 
of  subbase  treatment),  and  on  the  amount  of  rainfall. 

An  equation  to  predict  pumping  potential  was  obtained  from 
a  regression  analysis  of  the  pumping  indices  of  the  AASHO 
Road  Test  and  of  the  amount  of  energy  imposed  on  the  pave- 
ment by  traffic.  The  amounts  of  deformation  energy  imposed 
by  one  application  of  a  80.1  KN  (18K)  -single  axle  load  on 
a  limited  deflection  basin  were  determined  for  all  the 
sections  of  the  AASHO  Road  Test.  These  energies  were  com- 
puted using  a  finite  element  technique.  The  AASHO  pumping 
indices  were  normalized  to  account  for  the  fact  that  dif- 
ferent  slab   lengths   were  used.  This  was  accomplished  by 
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dividing  the  pumping  indices  by  the  number  of  joints  per 
30.5  m  (100  feet)  ;  this  gives  the  normalized  pumping 
index.  The  deformation  energy  applied  on  a  deflection 
basin,   with  deflections  larger  than  20  mils,  was  computed 


whe  re 


energy  =   I  A    k    w 
i=l 


A    =  nodal  area 

k    =  nodal  subgrade  modulus 


=  nodal  deflection 


(5.13) 


n  =  number  of  nodes  with  deflections  larger  than   20 


mi  Is 


The  values  of  deformation  energy  for  different  thicknesses 
we  re  as  f o 1 lows : 


thickness 


deformation  energy 


89.  mm  (3.5  in..) 

127.  mm  (5.0  in.) 

165.  mm  (6.5  in.) 

203.  mm  (8.0  in.) 

24  1.  mm  (9.5  in.) 

279  .  ram(ll  .0  in.) 

318.  mm(12.5  in.) 


27.370  j  (2A2.25  in. -lb) 

10.072  j  (  89.15  in. -lb) 

2.728  j  (  24.15  in. -lb) 

0.986  j  (  8.73  in. -lb) 

0.309  j  (  2.74  in. -lb) 

0.095  j  (  0.84  in. -lb) 

0.028  j  (  0.25  in. -lb) 


102 

With  these  values  and  using  a  least  square  approximation, 
a  regression  equation  between  the  values  of  log(total 
energy)  and  the  values  of  ln(NPI)  was  obtained  for  202 
cases.  The  resulting  coefficient  of  correlation  was  R  ■ 
0.78.  The  correlation  equation  can  be  expressed  as  fol- 
lows : 

NPI  =  exp  [-2.884  +  1.652  log ( tot al / 10000  )  ]   (5.14) 
whe  r e : 

2 

NP I=norma lized  pumping  index  (in    ) 

total  =  (energy  per  application)  (ESAL)   (in. -lb) 

This  equation  is  implemented  in  the  computer  method  of 
analysis,  to  determine  the  value  of  pumping  index  for  a 
pavement  section   given  certain  traffic  volume. 


5.3.4  Computer  Implementation 
In  this  section,  the  computer  implementation  of  the  non- 
linear method  for  analysis  of  rigid  pavements  is 
presented.  The  computer  implementation  is  illustrated 
with  the  flow  chart  shown  in  Fig.  5.7.  In  general,  the 
computer  scheme  is  an  Iterative  process  in  which  the  pro- 
perties of  the  pavement  section  are  modified  during  each 
iteration,  in  accordance  with  the  level  of  stresses  and 
number  of  load  repetitions  that  occur  at  the  critical  sec- 
tion of  the  pavement.  In  the  first  iteration,  stresses, 
strains   and   deflections  are  calculated  for  all   elements 
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assuming  the  initial  properties  of  the  pavement.  In  subse- 
quent iterations,  the  properties  of  each  individual  ele- 
ment are  modified  given  the  stresses  and  deflections  that 
have  taken  place  in  that  particular  element  and  in  accor- 
dance with  the  fatigue  behavior  of  the  element. 

The  computer  scheme  presented  herein,  called  PMARP,  is 
based  on  modifications  of  the  computer  code  ILLISLAB 
developed  previously  at  the  University  of  Illinois  [9], 
The  computer  method  presented  herein,  is  a  Fortran  pro- 
gram with  a  main  subroutine  for  memory  allocation,  and 
eight  subroutines  with  the  following  functions: 

Subroutine  Slab.  This  subroutine  is  to  read  all  the 
pavement  properties  and  geometry,  to  call  all  the  other 
subroutines  and  to  carry  out  the  iterative  process. 

Subrout  ine  Stiff;  calculation  of  all  stiffness 
matrices  given  the  updated  pavement  properties. 

Subroutine  Up t  ri;  inversion  of  stiffness  matrices. 

Subroutine  Load ;  allocation  of  boundary  conditions 
and  applied  forces. 

Subroutine  Pi  sp ;  solution  of  finite  element  equations 
and  calculation  of  strains  and  displacements. 


Subroutine  Stress;  calculation  of  stresses. 
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suDr ouhnes 


C  st*st  ) 


'NPUT    PAVEMENT 

GEOMETRY 


INPUT     ALL 
INITIAL    PROPERTIES 


COMPUTE 

STIFFNESS    PROPERTIES 


INVERT 
STIFFNESS    MATRIX 


SPECIFY 

BOUNDARY    CONDITIONS 
AND    APPLY    FORCES 


SOLVE  F.E.  EQUATIONS 
AND  OBTAIN  STRAINS 
AND    DISPLACEMENTS 


PUMPING    DETERMINATION 


COMPUTE     STRESSES 


OBTAIN  FATIGUE    DAMAGE 

AND    MODIFY   PAVEMENT 

PROPERTIES 


(       ST0"      ) 


Figure  5.7  Flow  chart  for  computer  scheme  of  nonlinear 
finite  element  method  for  analysis  of  rigid  pavements. 


105 

Subrou t  ine  Check;  modification  of  pavement  proper- 
ties. 

Subrout  ine  Pump;  calculation  of  pumping  index  and 
allocation  of  support  conditions.  A  detailed  description 
of  the  method  of  finite  elements  can  be  found  elsewhere 
[46,33,38],  In  this  section,  detailed  descriptions  of  the 
iterative  process,  the  modification  of  pavement  proper- 
ties, modification  of  load  transfer  efficiency,  develop- 
ment of  voids,  and  final  output  are  presented. 

I terat ive  Process  .  Referring  to  the  flow  chart  in  Fig. 
5.7,  the  computation  process  consists  of  a  specified 
number  of  iterations  during  which  the  following  operations 
are  executed: 

a.  Determination  of   stiffness   matrices   using   updated 
pavement  properties 

b.  Inversion  of  updated  stiffness  matrices. 

c.  Specification  of  boundary  conditions  and   application 
of  nodal  forces. 

d.  Solution  of  F.E.  equations  and  calculation  of  strains 
and  nodal  displacements. 


e.   Determination  of  pumping  index  in  third  iteration  and 
allocation  of  support  conditions. 
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f.  Calculation  of  stresses. 

g.  Calculation  of  amount  of  damage  to  pavement,  decay  in 
pavement  properties,  and  determination  of  updated 
pavement  properties. 

h.   Print  results  and  stop. 

The  nonlinear  analysis  is  carried  out  by  means  of  an 
iterative  process.  The  stress-  and  strain-dependent  pro- 
perties of  structural  elements  are  updated  in  each  itera- 
tion. In  the  computer  scheme,  the  total  number  of  load 
repetitions  is  applied  by  steps.  For  each  step,  several 
iterations  are  performed  to  obtain  convergence  between  the 
assumed  values  of  stiffnesses  and  the  values  calculated 
after  the  fatigue  damage.  The  total  number  of  iterations, 
is  equal  to  the  number  of  loading  steps  multiplied  by  the 
number  of  iterations  per  step.  During  each  iteration, 
stresses,  strains,  and  deflections  are  calculated  assuming 
the  current  stiffness  properties  and  the  current  number  of 
load  repetitions.  The  stiffness  matrices  of  the  slab,  load 
transfer  elements,  and  sublayer  elements  are  updated 
correspondingly. 


Modif  ica t  ion  of  Concrete  Slab  Properties.-  The  concrete 
slab  is  divided  into  rectangular  plate  elements.  Two 
options  to  analyze  these  elements  are  provided:  (1)  iso- 
tropic  and    (2)  orthotropic  elastic  thin  plate  elements. 
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Only  the  orthotropic  case  is  explained  herein;  the  iso- 
tropic case  is  a  special  case  of  the  orthotropic  one 
for  which  the  mechanical  properties  are  the  same  in  all 
the  horizontal  directions. 

The  capability  to  analyze  orthotropic  plates  has  been 
developed  to  take  into  account  the  presence  of  different 
amounts  of  reinforcement  in  the  X  and  Y  directions.  The 
stiffness  matrix  of  each  plate  element  is  calculated  as 
shown  in  equation  (5.7)  ,  which  is  a  function  of  dimen- 
sions and  mechanical  properties  of  the  element.  The 
mechanical  properties  are  defined  in  terms  of  E  and 
E   .  These  values  are  the  equivalent  values  of  the  moduli 

y 

of  elasticity  of  the  composite  section.  If  the  sections 
subjected  to  bending  in  X  and  Y  directions  are  composed 
of  concrete  and  steel  bars,  the  equivalent  moduli  are 
obtained  as  follows: 

Equivalent  Modulus  of  Elasticity: 


E  =  E  +E   =■ 
c   si 


(5.15) 


For  the  X  and  Y  directions 


E  -  E  +0.75  p   (  -  )2   E 

x     c  rx  t       s 

E  =  E  +0.75  p   (  -  )2   E 

y     c  *y  t       s 

Whe  re 


E   ■  Modulus  of  elasticity  in  X  direction 
x 

E   -  Modulus  of  elasticity  in  Y  direction 

y 
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E   = 


E   = 


Xl  = 

l2    ' 
d  = 

t  = 


Modulus  of  elasticity  of  concrete 
Modulus  of  elasticity  of  reinforcing  bars 
Percentage  of  reinforcement  in  X  direction 
Percentage  of  reinforcement  in  Y  direction 
Moment  of  inertia  of  the  concrete  section 
Moment  of  inertia  of  one  rei nf or ci ngbba r 
Diameter  of  rei nf or cingbba r 
Thickness  of  slab 


The  values  of  moduli  of  elasticity  of  concrete  are  modi- 
fied as  the  loads  are  imposed.  The  modification  is  carried 
out  in  accordance  with  the  level  of  stresses  and  the 
number  of  load  repetitions,  making  the  following  assump- 
tions : 

(i).-The  moduli  of  elasticity  of  slab  elements 
in  the  X  and  Y  directions  are  equal  to  the 
values  of  the  composite  sections  in  the  X  and  Y 
directions  as  given  in  equations  5.8. 

(ii). -Decay  in  strength  and  modulus  of  elasti- 
city in  the  composite  section  takes  place  in  the 
concre t e  on ly  . 

(iii).-A  flexural  endurance  curve  is  known  for 
the  particular,  or  similar,  concrete. 


iv.-The  Palmgren-Mine r  linear  cumulative   damage 
hypothesis  is  applicable. 
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(v).-Rest  periods  between  loading  cycles  have  no 
effect   on   the   fatigue   behavior  of   concrete. 
(vi).-The  relations  among  the  modulus  of  elasti- 
city, the  compressive  strength,  and  the  flexural 
strength  of  concrete  are  known. 
Under  these  assumptions,  the  strength  and  modulus  of  elas- 
ticity  of   a   slab  element  subjected  to  certain  number  of 
stress  repetitions  can  be  calculated   by   determining   the 
remaining   fatigue  life  of  the  concrete  .  Three  conditions 
may  occur: 


if  f   >  1  .0  »>  cracking  takes  place 

if  f  .   <  f   <  1  .0  =>  fatigue  damage  takes  place 
min     r 


if  f    <  f  .   =>  no  damage  occurs 

r      min  " 


Whe  re 


f   =  relative  stress 

r 


MR 


o     ■  maximum  principal  stress 
max 

MR  =  modulus  of  rupture  of  concrete 

f  .   ■  minimum  stress  necessary-to-produce  damage 
min 


Depending  on  the  value  of   f    ,   the   properties   of   the 
slab   element   are   modified   because   cracking   occurs 
(f   >  1.0)     ,    or    because    of    fatigue    decay 

(f  .   <  f   <  1.0)    .   When  cracking  occurs,  the  thickness 

*  min     r  ° 

of  the  composite  section  is  modified  so  that  moment 
equilibrium  in  the  cracked  section  is  satisfied.  This  is 
accomplished  in  the  computer  program  by  means  of  a   trial- 
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and-error  process  to  find  the  position  of  the  neutral  axis 
in  the  composite  section.  The  number  of  iterations 
required  depends  upon  the  approximation  error  desired. 
After  the  neutral  axis  is  found,  the  new  thickness  of  the 
slab  element  is  calculated  as  the  initial  thickness  minus 
the  depth  of  the  crack,  which  is  assumed  to  extend  through 
the  entire  tension  zone. 

When  only  fatigue  damage  occurs  the  properties  of  the  con- 
crete slab  ,  at  the  particular  point,  are  calculated  using 
the  endurance  curve.  For  example,  taking  a  concrete  ele- 
ment with  initial  MR  -  3447  kPa   (  500  psi)   and  subjected 

to  1000,000  load  repetitions   a  relative  stress   f    -  0.7 

r 


( 


o  =2413kPa       (350psi)  )       is       produced.         If       the 

max  e  r 


endurance  curve  of  this  concrete  is  the  one  shown  in  Fig. 
5.5,  the  number  of  load  repetitions  to  failure  is 
500,000.  Therefore,  the  remaining  fatigue  life  is  400,000 
repetitions,  .which  corresponds  to  f  ■  0.75  Thus,  the 
updated  MR  is  2413. /0. 75  -  3217  kPa  (467  psi).  With 
this  value  of  the  modulus  of  rupture,  the  strength  and  the 
modulus  of  elasticity  are  recalculated  knowing  the  rela- 
tions among  these  parameters.  Then,  the  properties  of  the 
slab  at  the  particular  point  are  recalculated  using  the 
composite  section  as  given  in  equations  (5.8). 


Modif  ica t ion  of  Load  Transfer  Ef  f lcency  .-  Unfortunately, 
very  little  experimental  research  has  been  conducted  to 
investigate  the  behavior  of  dowel  bars  under  the  action  of 
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repetitive  loading.  Therefore,  the  method  proposed  herein 
to  modify  the  load  transfer  efficiency  as  a  function  of 
the  number  of  load  cycles,  is  a  simplified  but  practical 
approach  based  on  the  available  experimental  data.  It  is 
intended  that  the  model  proposed  herein  shall  be  con- 
sidered as  a  guide  which  can  be  modified  or  adjusted 
later,  especially  if  experimental  data  are  available  for 
a  particular  case. 

The  method  for  modifying  the  values  of  load  transfer  effi- 
ciency is  based  on  the  experimental  results  presented  by 
Teller  and  Cashell  [28].  These  results  are  used  to  obtain 
a  trend  of  behavior  for  the  variation  of  load  transfer 
efficiency  as  a  function  of  load  repetitions  and  level  of 
stresses.  In  the  computer  implementation,  the  value  of 
load  transfer  efficiency  is  calculated  as  a  function  of 
the  geometry  of  the  dowels  and  the  concrete  properties  by 
means  of  a  parameter  called  Dowel-Concrete  Interaction 
Factor.  The  value  of  this  factor  is  given  in  [22] 


DCI  = 


G0.75    2.5 
d 


0.041D,0*75  +  0.0004G0,25  w. 
d  J 


(5.16) 


whe  re : 


G=  modulus  of  dowel  support 


D,  =  dowel  diameter 
d 

w.  =  joint  width 

J 
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The  initial  value  of  DCI  is  changed  by  multiplying  it  by  a 
reduction  fator,  RF  ,  obtained  from  a  regression  equation 
obtained  from  data  reported  in  reference  [28], The  value 
of  RF  is  calculated  as: 


R 


F  =  0.268  -  0.0457  log(n)  +  1.123  f  . 

"  rl 


(5.17) 


whe  re 


RF  =  reduction  factor 

n  ■  number  of  load  repetitions 

f  ,  =  relative  load  acting  on  dowel 
rb 


rb     P 
c 

Pn    ■    concentrated    load    acting    on    dowel 

Dd    U~    °d    } 
Pc    =    crack    load    =    3/2    (1+w/L) f  t 

t=thickness    of     slab 

L    =    embedded    length    of    dowel 

1/2 
f       =     10    (    f       )1/Z 

t  c 

The  expression  for  crack  load  is  given  by  Marcus  [27]. 
The  values  of  DCI  are  calculated  in  each  iteration  in  the 
computer  process  depending  on  the  number  of  load  cycles 
and  the  level  of  loads  applied  to  the  dowels.  If  more 
experimental  data  become  available,  a  more  accurate 
expression  for  the  value  of  RF  can  be  implemented  in  the 
computer  scheme. 


Pumping  Development .  Pumping  development  and  allocation  of 
support    conditions    are   included   in   this   method   of 
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analysis.  Pumping  index,  determined  by  means  of  equation 
5.14,  is  used  to  determine  the  volume  of  pumped  material 
and  to  determine  the  area  of  voids.  The  area  of  voids  is 
allocated  in  the  pavement  structure  by  removing  the  sup- 
port at  several  nodes  to  cover  the  calculated  void  area. 
Allocation  of  nodes  with  zero  support  is  provided  consid- 
ering the  nodes  with  the  highest  deflections.  The  pro- 
cedure is  as  follows: 

1.  Compute  amount  of  deformation  energy  per  one  ESAL 
application. 


whe  re : 


i"n 
energy  =    I   A   k   w 
i  =  l 


A    =  nodal  area 


k    =  nodal  subgrade  modulus 


w    =  nodal  deflection 


(5.18) 


2.  Compute  total  applied  energy. 

total=(energy  per  one  ESAL)x(ESAL) 

3.  Compute  normalized  pumping  index. 

NPI  -  exp(-2. 8847  +  1  .6521  log( total/ 10000 ) ) 


4.  Compute  total  volume  of  pumped  material   per  slab 
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vol=(NPI ) (slab  length) (number  of  joints  per  100  ft) 

5.  Obtain  void  area. 

a rea=vol/ (a ve rage  void  depth) 

6.  Allocate  nodes  with  zero  support  to  cover  void  area. 

Allocation  of  nodes  with  zero  support  is  done  by 
sorting  and  arranging  the  nodes  in  descending  order  in 
accordance  with  their  deflections.  Given  the  void  area, 
nodes  with  zero  support,  with  their  corresponding  area, 
are  assigned  until  the  total  area  covered  by  the  nodes  is 
equal  to  that  of  the  void.  After  zero  support  conditions 
are  allocated  to  certain  nodes,  these  nodes  remain  with 
no  support  through  all  the  iterative  procedure. 


The  determination  of  the  area  without  support  (area  of 
void)  is  accomplished  by  dividing  the  volume  of  the  void 
by  an  average  void  depth.  The  size  of  the  area  with  no 
support  is  also  a  critical  parameter  in  the  analysis. 
Unfortunately,  there  is  only  minimal  detailed  experimental 
evidence  concerning  the  actual  size  and  configuration  of 
the  voids.  Therefore,  the  assumption  of  an  average  void 
depth  represents  an  attempt  not  to  describe  the  geometri- 
cal properties  of  a  void,  but  to  provide  a  parameter  to 
obtain  an  unsupported  area,  which  is  in  accordance  with 
measured  values  obtained  by  means  of  deflection-based 
experimental   methods.   In  the  computer  code,  the  value  of 
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the  average  void  depth  is  a  necessary  part  of  the  input. 
In  the  absence  of  experimental  evidence,  the  analysis  must 
necessarily  be  based  on  a  reasonable  assumption  for  this 
parameter.  More  experimental  evidence  of  the  size  of 
voids,  or  of  functionally  equivalent  sizes,  is  urgently 
needed . 

Resilient  Modulus  of  Subgrade  Reaction.-  The  structural 
analysis  of  rigid  pavement  generally  has  been  carried  out 
assuming  a  static  constant  value  for  the  stiffness  of  the 
subbase.  This  type  of  analysis  assumes  a  linearly  elastic 
behavior  of  the  subbase.  In  the  finite  element  analysis  of 
rigid  pavements,  this  assumption  is  represented  by  a 
Winkler  sublayers  consisting  of  a  series  of  spring  ele- 
ments. However,  even  though  this  assumption  provides 
fairly  accurate  results,  the  actual  behavior  of  the  sub- 
base  is  complex  and  the  nature  of  the  traffic  load  is  not 
static.  A  better  representation  can  be  obtained  by  using 
an  impact  resilient  modulus  of  subgrade  reaction,  k 
[44]  In  accordance  with  this  assumption,  the  value  for  the 
stiffness  at  a  specific  point  in  the  sublayer  depends  on 
the  value  of  the  deflection  at  that  point.  This  results 
in  a  non-linear,  load-deflection  relationship. 


The  method  of  analysis  proposed  herein  has  the  capa- 
bility for  considering  a  composite  resilient  value  for  the 
stiffness  of  the  sublayers.  In  Figure  5.8,  a  typical 
curve   relating  the  values  of  resilient  modulus  versus  the 
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deflections  Is  shown  [44],  In  the  computer  implementa- 
tion, first,  the  values  of  k  are  introduced  as  part  of 
the  input,  by  specifying  the  k  values  corresponding  to 
several  deflections.  In  the  iterative  process,  the  value 
of  the  modulus  of  subgrade  reaction  assumed  at  each  point 
is  calculated  from  the  deflection  occurring  at  that  point. 
In  the  iterations,  the  initially  assumed  value  of  k  is 
adjusted,  as  the  iterations  proceed  in  accordance  with  the 
resulting  deflection.  Thus,  a  nonlinear  pressure- 
deflection  ratio  is  obtained. 


5_.4_  Conclus  ions  and  Re  commend  at  ions 

A  Non-Linear  Finite   Element   Method,   called   PMARP,  for 
analysis   of   rigid   pavements  subjected  to  fatigue  damage 

caused  by  traffic   has   been   presented.    The   method  is 

implemented   in   an  iterative  computer  scheme  in  which  the 

pavement  properties  are  modified  in   accordance   with  the 
fatigue  damage  produced  by  repetitive  loading. 
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Plate  Deflections,  mil* 


Figure    5.8    Deflection-dependent    behavior    of    resilient 
composite    modulus    of    subgrade    reaction. 
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In  a  rigid  pavement,  most  of  the  structural  support  is 
provided  by  the  concrete  slab.  This  slab,  then,  under  the 
action  of  traffic  is  subjected  to  many  number  of  load 
repetitions  along  the  service  life  of  the  pavement.  The 
load  repetitions  produce  gradual  damage  in  all  the  pave- 
ment components  until  eventually  cracking  and  faulting 
occurs  producing  the  failure  of  the  structure.  In  jointed 
concrete  pavements,  the  devices  which  transfer  the  load 
across  the  joints,  specially  transverse  joints,  also  lose 
their  efficiency  along  the  life  of  the  pavement.  Fatigue 
damage  imposed  to  the  concrete  slab  reduces  its  strength 
and  modulus  of  elasticity  because  of  the  occurrence  of 
mi  croc  rack s  .  Similarly,  fatigue  damage  caused  by  repeti- 
tive loading  reduces  the  load  transfer  efficiency  of  the 
dowel  bars  by  increasing  the  degree  of  looseness  of  the 
bars.  These  decays  in  the  properties  of  the  pavement  com- 
ponents modify  the  response  of  the  pavement  to  a  given 
load.  Stresses  and  deflections  calculated  assuming  the 
initial  constant  stiffness  values  differ  from  those 
obtained  by  using  reduced  values  which  account  for  fatigue 
damage.  Assuming  fatigue  damage,  results  closer  to  actual 
field  conditions  are  obtained. 


In  the  computer  implementation  of  the  method  of  analysis 
presented  herein,  the  modifications  to  the  properties  of 
the  concrete  slab  are  based  on  the  flexural  behavior  of  a 
particular   concrete   specified  as  part  of  the  input  data. 
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Similarly,  the  modification  of  the  values  of  the  resilient 
modulus  of  subgrade  reaction  is  made  on  the  basis  of 
values  of  k  specified  in  the  input  data.  For  the 
modification  of  the  load  transfer  efficiency,  however, 
only  the  initial  value  can  be  specified  in  the  input  data. 
The  fatigue  behavior  is  controlled  by  an  expression  built 
into  the  computer  scheme.  Nevertheless,  this  expression 
can  be  modified  by  changing  the  values  of  the  coefficients 
used  in  the  computer  implementation.  In  any  event,  accu- 
racy of  the  absolute  values  obtained  with  the  method  pro- 
posed herein,  depends  on  the  accuracy  of  the  information 
provided  in  the  input;  especially  the  information  concern- 
ing the  fatigue  behavior  of  the  different  pavement  com- 
ponents. The  relative  results,  however,  are  more  accurate 
than  those  obtained  by  less  comprehensive  methods,  and  are 
especially  useful  in  comparative  studies. 


The  support  conditions  are  modified  as  a  consequence  of 
pumping.  Pumping  of  material  from  under  the  concrete  slab 
produces  areas  of  low  or  no  support.  A  pumping  model 
developed  using  the  AASHO  Road  Test  data  has  been  imple- 
mented in  the  method  of  analysis  presented  here.  By  means 
of  this  pumping  model,  the  support  conditions,  of  a  pave- 
ment section  being  analyzed,  are  modified  automatically  by 
means  of  a  subroutine  called  "PUMP".  The  modification  of 
the  support  conditions  depends  upon  the  amount  of  pumping 
which  is  also  calculated  by  means  of  this  subroutine. 
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CHAPTER  6   NUMERICAL  STUDY 

§_•!_   Int  roduc  t  Ion 

In  this  chapter,  numerical  results  obtained  by  means  of 
the  method  of  analysis  described  in  the  previous  chapter 
are  presented.  These  numerical  results  are  analyzed  for 
the  following  purposes:  (1)  to  verify  the  sensitivity  of 
the  proposed  method  to  the  variations  of  different  pave- 
ment parameters,  and  (2)  to  provide  results  for  estima- 
tion of  pavement  performance. 

First,  a  sensitivity  study  is  presented  in  this  chapter  by 
evaluating,  for  the  pavement  section  shown  in  Fig.  6.1, 
deflections  and  amount  of  fatigue  damage  caused  by  varia- 
tions in  the  following  parameters: 

1  .-Slab  thickness  . 

2 .-Amount  of  traffic. 

3 .-Amount  of  reinforcement. 

4  .-Type  of  sublayer. 

5. -Endurance  behavior  of  concrete. 

In  the  second  part,  results  in  terms  of  amount  of  pumping, 
amount  of  fatigue  damage  in  the  slab,  size  of  damaged 
area,  and  maximum  deflections  are  obtained  for  different 
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Figure    6.1    Two-slab    pavement    section    with    doweled    joint 


122 


values  of  the  following  pavement  parameters: 
l.-Slab  thickness. 
2 .-Amount  of  traffic. 
3  .-Type  of  sublayer. 

In  this  chapter,  the  terra  sublayer  refers  to  the  layers 
(subbase  and  subgrade)  below  the  concrete  slab.  The  compo- 
site resilient  modulus  of  subgrade  reaction  is  an 
equivalent  value  which  corresponds  to  the  stiffness  of  the 
sublayers  (subbase  and  subgrade).  Therefore,  the  term 
type  of  sublayer  refers  to  the  sublayers  whose  composite 
stiffnesses  are  given  in  Fig.  6.2. 

6_.2_  Sensitivity  Study 
The  effects  of  several  pavement  parameters  on  the  results 
obtained  using  the  method  of  analysis  presented  in  Chapter 
5  are  shown  here.  These  effects  are  quantified  in  terms  of 
maximum  deflection  and  amount  of  damage.  In  this  case, 
full  contact  is  assumed  to  remain  at  all  load  repetitions. 
The  pavement  section  shown  in  Fig.  6.1  is  analyzed  assum- 
ing the  following  pavement  parameters: 

l.-Slab  thickness  (t). 

t=150.  mm  (6.  in.) 

t=178.  mm  (7.  in.) 

t=203.  mm  (8.  in.) 

t=254.  mm  (10.  in.) 
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Figure  6.2  Composite,  deflection  dependent  resilient 
moduli    of  subgrade  reaction . 
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Figure  6.3  Endurance  curves  for  concrete  under  flexure. 
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2  .-Equi va le nt  single-axle  load  applications  (ESAL). 
ESAL=100 ,000 
ESAL=1 ,000,000 
ESAL=10 ,000 ,000 

3  .-Amount  of  reinforcement  (p). 
p  =  0  .% 
p  =  0  .05% 
p=0.25% 

A  .-Type  of  sublayer  (Fig.  6.2). 
Soft. 
Medium . 
Stiff. 

5. -Endurance  behavior  of  concrete  (Fig.  6.3). 

Type  1. 

Type  2. 


6_.3^  Pavement  Performance 
Results  in  terms  of  pumping  index,  amount  of  fatigue   dam- 
age  in   slab,   decay   in   load   transfer   efficiency,  and 
deflections  are  obtained  for  the  following  parameters: 

1  .-Slab  thickness  (t) 

t=150  mm  (6  in.) 

t=203  mm  (8  in.) 

t=254  mm  (10  in.) 

t=305  mm  (12  in.) 
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2  .-Equivalent  single-axle  load  applications  (ESAL) 
ESAL=100 ,000 . 
ESAL=1 ,000,000. 
ESAL=10,000,000  . 


3  .-Type  of  sublayer 

Soft 

Me  diura 

Stiff 

4. -Average  void  depth 


=  6.3  mm  (0.25  in.) 


6_.4_  Results 
Results  in  terms  of  amount  of  damage  and  deflections  are 
shown  in  Tables  6.1  to  6.10.  These  results  were  obtained 
for  the  case  in  which  fatigue  damage  in  the  load  transfer 
devices  and  pumping  are  not  considered.  Deflections  are 
considered  at  corner  of  the  loaded  slab  (  node  57  in  fin- 
ite element  mesh)  produced  by  a  177.9  kN  (40-K)  single 
axle  load  applied  after  the  pavement  has  been  under  the 
given  ESAL.  Amount  of  damage  in  the  pavement  is  quanti- 
fied in  terras  of  total  damaged  area,  and  in  terms  of 
fatigue  damage  calculated  as  follows: 


i  =  n 


Damage  =   I  A    D 
i-1  1  X 


(6.1) 
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where,  A  is  the  partial  damaged  area,  and  D  is  the 
decay  of  either  the  modulus  of  elasticity  or  the  effective 
thickness  of  the  slab  given  by: 


Ec   -  Ec, 
o f_ 

Ec 


(6.2) 


or 


Di  = 


o     f 


(6.3) 


Ec    and   Ec,   are  the  initial  and  final   values   of   the 
o  r 

modulus  of  elasticity  of  the  slab.  Similarly,  t  and 
t  are  the  initial  and  final  effective  slab  thicknesses. 
Tables  6.11  to  6.16  show  the  results  in  terms  of  amount  of 
fatigue  damage  in  the  slab  and  deflections  for  a  177.9  KN 
(40-K)single  axle  load.  In  this  case,  fatigue  damage  in 
the  load  transfer  devices  is  included.  Additionally,  the 
change  in  load  transfer  efficiency  is  also  presented. 
Load  transfer  efficiency  is  quantified  as: 


1   i=n  62 
LTE  =  -    I  -t=- 

n   i=l  61 


(6.4) 


where:  LTE  is  the  value  of  load  transfer  efficiency,  and 
6.  and  6„  are  the  values  of  the  deflections  at  the 
loaded  and  opposite  slabs  respectively.  The  quantity  n  is 
the  number  of  dowels  considered.  The  LTE  results  shown  in 
Tables  6.11  to  6.16  correspond  to  nodes  57  and  65. 
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Table  6.1   Effects  of  thickness  and  ESAL  on  damaged  area 

and  amount  of  fatigue  damage  (p=0.05%,  soft  sublayer, 

and  concrete  type  1). 


ESAL 

THICKNESS 

ISO.    mm 

178.    mm 

20  3.    mm 

25  4 .    mm 

100,000. 

area    m2(in2 ) 
daaage 

2.70    (4200) 
1006.63 

1.80    (2800) 
13.57 

0.90    (1400) 
2.69 

0. 
0. 

1,000,000. 

area    a2(in2) 
dama  ge 

2.70    (4200) 
1009.10 

1  .93    (3000) 
75.76 

0.90    (1400) 
26.70 

0. 
0. 

10,000,000. 

area    a2  (  in2  ) 
dacage 

2.70    (4200) 
1270.16 

1.93    (3000) 
597.20 

0.90    (1400) 
257  .00 

0. 

0. 

...     J 

Table  6.2   Effects  of  thickness  and  ESAL  on  deflections 
(p=0.05%,  soft  sublayer,  and  concrete  type  1) 


ESAL 

THICKNESS 

150.    mm 

178.    mm 

203.    mm 

254 .     mm 

100,000. 

2.69 
(  106  .0) 

1  .38 
(54.4) 

1  .16 
(45.8) 

0.88 
(35.0) 

1  ,000,000. 

2.70 
(  106.4) 

1  .38 
(54.5) 

(45.9) 

0.88 
(35.0) 

10,000  ,000. 

2.76 
( 109 .0) 

1.41 
(55.7) 

1  .18 
(46.4) 

0.88 
(35.0) 

de  f lee 1 1 ons 

In    mm    (mils) 

Table  6.3   Effects  of   ESAL  and  percentage  of 
reinforcement  (p)  on  damaged  area   and   amount 
of  fatigue  damage  (t=l78  mm,  soft  sublayer, 
and  type  1  concrete) 
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P* 

ESAL 

100,000 

1  ,000,000 

10,000,000 

oz 

area    mz ( In2     ) 
damage 

1.80    (2800) 
13.58 

1.93    (3000) 
75.88 

1.93  (3000) 
598.1 

0.05Z 

a  rea    m2 ( in'  ) 
damage 

1.80    (2800) 
13.57 

1.93    (3000) 
75.76 

1.93  (3000) 
597.2 

0.25Z 

area    m    ( In2 ) 

damage 

1.80    (2800) 
13.52 

1.93    (3000) 
75.46 

1.93  (3000) 
597  .1 

Table  6.4   Effects  of  ESAL  and  percentage  of 
reinforcement  (p)  on  deflections  (t=178  mm, 
soft  sublayer,  and  type  1  concrete) 


pz 

ESAL 

100,000 

1,000,000 

10  ,000,000 

oz 

0.05Z 
0.25Z 

1  .38 
(54.4) 

1.38 
(54.5) 

1  .41 
(55.7) 

1.38 
(54.4) 

1.38 

(54.5) 

1  .41 
.    (55.7) 

1  .38 

(54.4) 

1.38                      1.41 
(54.5)                  (55.6) 

deflections    In    mm    (mil) 
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Table  6.5   Effects  of  ESAL  and  percentage  of  reinforcement 
(p)  on  damaged  area  and  amount   of    fatigue     damage 
( t  =  1 5  2  mm,  soft  sublayer,  and  type  1  concrete). 


p: 

ESAL 

100,000 

1,000,000 

10 ,000,000 

oz 

area    m2(in2     ) 
damage 

3.22    (5000) 
2755.6 

3.22    (5000) 
2840.00 

3.48    (5400) 
3428.0 

0.05Z 

area    m2 (in2  ) 
damage 

2.70    (4200) 
1006  .6 

2.70    (4200) 
1009.1 

2.70    (4200) 
1270.10 

0.25Z 

a  rea    m2 ( In2  ) 
damage 

2.70    (4200) 
991  .6 

2.7    (4200) 
1003.3 

2.7    (4200) 
1073.7 

Table  6.6   Effects  of  ESAL  and  percentage  of  reinforcement 
(p)  on  deflections  (t=152  mm,  soft  sublayer, 
and  type  1  concrete) 


P* 

ESAL 

100,000 

1,000,000 

10,000,000 

oz 

4.31 
(169.7) 

4.30 

(169.4) 

4.51 
(177.7) 

-      ■ 
0.05Z 

2.69 
(106.0) 

2.70 
(  106.4) 

2.77 
(109.0) 

0.25Z 

6.69 
( 106  .0) 

2.70 

( 106.4) 

2.71 
(106.7) 

def  lee 

Clons     In    mm    (ml 

1) 
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Table  6.7  Effects  of  ESAL  and  type  of  sublayer  on  damaged 
area  and  amount  of  fatigue  damage  ( t  =  1 7 8  mm,  p=0.05%,  and 
type  1  concrete). 


I 

1 

ESAL 

1 
Sublayer 

100,000 

1  ,000,000 

10,000,000 

[STIFF 

■a  rea    a2(  In'     ) 
da  na ge 

1.03    (1600) 
13.46 

1  .03    ( 1600) 

43.63 

1  .03    (  1600) 

243.3 

; MEDIUM 

are*    a2  ( In    ) 
dama  ge 

1.29    (2000) 
13  .50 

1.29    (2000) 
47  .34 

1.41     (2200) 
474.2 

1 

|  SOFT 

i 

area    a2  (In2  ) 
dama  ge 

1.80    (2800) 
13.57 

1  .93    (3000) 
75.76 

1.93     (3000) 
597  .2 

Table  6.8   Effects  of  type  of  sublayer  and  ESAL 
on  deflections  (t=178  mm,  p=0.05%,  and 
type  1  concrete  )  . 


ESAL 

Sublayer 

100,000 

1 

,000,000 

10,000,000 

STIFF 

0.82 
(32.5) 

i 
i 

0.82 
(32.5) 

0.83 
(33.0) 

MEDIUM 

1.11 
(43.8) 

j 

1.11 
(43.8) 

1 

1  .13 
(44.6) 

SOFT 

1  .38 
(54.4) 

i 
i 

| 

1  .38 
(54.5) 

1  .41 
(55.7) 

de  f le  ct  lo 

is     In    om    (all) 

i 
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Table  6.9   Effects  of  ESAL  and  type  of  concrete  on 

damaged  area  and  amount  of  fatigue  damage 

C  t" 17  8  mm,  p-0.05%,  and  soft  sublayer). 


ESAL 

100,000 

1,000,300 

10,000,000 

TYPE    OF 

CONCRETE 

area    m2 ( In2    ) 
Type    1 

1.80    (2800) 

1.93    (3000) 

1.93    (3000) 

damage 

13.57 

75.76 

597.2 

a  rea    m2 ( in2  ) 
Type    2 

1.80    (2800) 

1.80     (2800) 

1.80    (2800) 

damage 

7.97 

18.63 

105.7 

Table  6.10   Effects  of  type  of  concrete  and  ESAL   on 
deflections  (t=178  mm,  p=0.05Z,  and  soft  sublayer). 


TYPE    OF 

CONCRETE 

ESAL 

100,000 

1,000,000 

10,000,000 

1.38 
TYP»    1           1      (54.4) 

1 

1.38 
(54.5) 

1.41 
(55.7) 

!         1.38 
TyP«    2           1      (54.4) 

1.38 
(54.4) 

1.38 
(54.5) 

deflection*    In    mm    (mil) 
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The  amount  of  damage  is  very  sensitive  to  the  variation  of 
thickness  and  ESAL  (Table  6.1).  A  change  of  fatigue  damage 
from  an  average  of  1095  to  an  average  of  95  occurs  when 
the  thickness  varies  from  150  mm  (6  in.)  to  20  3  mm  (8 
in.).  The  variation  of  amount  of  fatigue  damage  as  a 
function  of  ESAL  is  more  pronounced  for  a  thickness  of  203 
mm  (8  in.)  than  for  150  mm  (6  in.).  The  values  of  fatigue 
damage  for  t=150  mm  (6  in.)  are  much  higher  however.  The 
deflection  values,  as  expected,  are  very  sensitive  to 
variations  in  thickness.  Deflections  remain  almost  con- 
stant for  different  ESAL  (Table  6.2). 


The  flexural  resistance  of  the  slab  is  controlled  in  the 
computer  scheme  by  specifying  the  percentage  of  reinforce- 
ment, and  the  fatigue  flexural  behavior  of  the  particular 
concrete.  The  effects  of  percentage  of  reinforcement  and 
flexural  behavior  of  concrete  are  given  on  Tables  6.3, 
6.5,  and  6.9.  The  percentage  of  reinforcement  has  a  small 
effect  on  the  amount  of  fatigue  damage,  while  the  flexural 
resistance  has  a  severe  effect.  Therefore,  when  reinforce- 
ment is  provided,  not  only  should  the  amount  of  reinforce- 
ment be  specified  as  input  in  the  computer  method  of 
analysis,  but  also  the  fatigue  flexural  behavior  should  be 
modified.  A  concrete  with  higher  flexural  resistance 
should  be  assumed  for  higher  reinforcement  percentages 
even  for  the  same  concrete  strength.  The  percentage  of 
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Table  6.11   Amount  of  fatigue  damage  and  damaged  area 
caused  in  pavement  section  with  decay  in  load  transfer 
efficiency  (p=0.05%,  soft  sublayer,  concrete  type  1). 


ESAL 

Thickness    (nam) 

178 

203 

228                       254 

100,000. 

1  r  ,     2  \ 

area    a  '  (  i  «■•■  ) 

1  .68(2600) 

0.90< 1400) 

0.26(400)       0.13(200) 

dama  ge    ( In2  ) 

16  .35 

2.34 

0.83                          0. 

1 ,000  ,000  . 

a  re  a    ta2  ( In2  ). 

1  .68( 2600) 

0.90( 1400) 

0.64( 1000)    0.13( 200) 

damage     ( in2 ) 

155.01 

47.05 

19.06                  3.78 

10,000,000. 

2  /  i     2  \ 
area    a-  (  i  n*  ) 

1 .68(2600) 

1  .03(  1600) 

0.64(1000)    0.13(200) 

daaage     ( in2 ) 

551 .44 

314.52 

1  18 .34               31.65 

Table  6.12   Amount  of  fatigue  damage  and  damaged  area 
caused  in  pavement  section  with  decay  in  load  transfer 
efficiency  (p=0.05%,  medium  sublayer,  concrete  type  1). 


Thickness    (mm) 

ESAL 

178                    203                         228                   254 

100,000. 
1 ,000  ,000  . 
10 ,000 ,000  . 

area    a2  (  In2  ) 
damage     ( in2 ) 

area    a2(U2) 
daoa  ge     ( In2 ) 

a  rea    a2 ( in2  ) 

daaa  ge     ( 1 n2  ) 

1.03(1600)     0.64(1000)     0.26(400)    0.13(200 
13.19                  1.51                  0.                            0. 

1.16(1800)     0.64(1000)     0.38(600)    0.13(200 
142.02               30.58                12.49                     3.84 

1.29(2000)     0.77(1200)     0.38(600)    0.13(200 
4'80.45             246.90               127.04                  32.27 
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Table  6.13  Amount  of  fatigue  damage  and  damaged  area 
caused  In  pavement  section  with  decay  in  load  transfer 
efficiency  (p=0.05%,  stiff  sublayer,  concrete  type  1). 


ESAL 


Thickness      (mm) 


178 


:o3 


25' 


100  ,000  . 


1 .000,000. 


10 ,000 ,000  . 


area    a2;in2)  0.77(1200)  0.64(1000)  0.26C400)  0.13(200) 

damage     'In2)  13.45  1.32                    0.35  0.35 

area    a2'in2)  0.77(1200)  0.64(1000)  0.26(400)  3.13(200) 

damage     (In2)  23.50  26.70                    9.12  3.84 

area    ::i'.n2)  0.90(1400)  0.64(1000)  0.39(500)  0.13(200) 

damage     :in2)  341.22  186.50               111.45  32.43 


Table  6.14   Effect  of  thickness  and  ESAL  on  deflections 
and  on  load  transfer  efficiency  (p-0.05%,  soft 
sublayer,  concrete  type  1). 


ESAL 

Thickness 

(mm) 

173 

203 

228 

254 

100,000. 

1  .-1(55) 
36.12 

1  .20(47) 
83.92 

1  .04(41) 
81  .32 

0.92( 36) 

79  .92 

1 ,000  ,000  . 

1.-4(56) 
33.32 

1  .22(48) 
81.32 

1.06(42) 
78.32 

0.94(37) 
75.82 

10,000,000. 

1.-9(59) 
31  .32 

1  .26(49) 
77.82 

1.10(43) 

73.12 

0.99(39) 
63.32 

deflections     i 

i    an    (mi  Is  ) 
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Table  6.15   Effect  of  thickness  and  ESAL  on  deflections 
and  load  transfer  efficiency  (p-0.05%,  medium  sublayer, 

concrete  type  1 ) . 


ESAL 

Thickness 

(mm) 

178                     20  3 

228 

254 

100,000. 

1.04(41)       0.89(35) 
82.02               79.22 

0.78( 30) 

76.4: 

0.69( 27) 
"3.92 

1 ,000  ,000  . 

1.06(42)         0.91(36) 
79.52               75.92 

0.79(31 ) 
72.3: 

0.71( 28) 
u8.92 

10 ,000  ,000  . 

1.10(43)         0.94(37) 
76.52               71.72 

0.83(33) 
66.2: 

0.75( 29) 
60.92 

def  lee  t  i  oris     1 

n    mm     (  mi  Is  ) 

Table  6.16   Effect  of  thickness  and  ESAL  on  deflections 
and  on  load  transfer  efficiency  (p-0.05%,  stiff  sublayer, 

concret  e  type  1 )  . 


ESAL 

Thickness      (mm) 

178 

203 

228 

254 

100,000. 

0.85(33) 
78  .92 

0.73(29) 
74.62 

0.64(25) 

72.42 

0.58(23) 
69  .52 

1  ,000  ,000  . 

0.87(34) 
7  5.82 

0.75(30) 
7  1  .82 

0.66(26) 
67.82 

0.60(24) 
64.02 

10,000,000. 

0.90(35) 

72.52 

0.78(31) 
66.92 

0.69(27) 
61.12 

0.63(25) 
55.22 

deflect!  one 

In    mm    (alls! 
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reinforcement  has  a  greater  effect  on  deflections  when 
thinner  pavement  slabs  are  considered  (Tables  6.4  and 
6.6). 

The  sensitivity  of  the  amount  of  damage  to  type  of  sub- 
layer is  higher  for  higher  values  of  ESAL  (Table  6.7). 
Moreover,  the  variation  of  the  amount  of  damage  as  a  func- 
tion of  ESAL  is  higher  for  lower  values  of  subgrade 
modulus.  The  values  of  the  deflections  are  not  very  sensi- 
tive to  the  ESAL  for  each  value  of  subgrade  modulus,  but, 
as  expected,  deflection  values  are  very  sensitive  to 
values  of  subgrade  modulus  (Table  6.8).  The  small  varia- 
tion of  deflections  as  a  function  of  ESAL  is  due  to  the 
fact  that  in  this  case  full  support  was  assumed  to  exist 
at  all  levels  of  ESAL. 


In  general,  higher  variations  of  the  amount  of  damage 
occur  when  decay  in  load  transfer  efficiency  is  considered 
(Tables  6.1  and  6.11).  When  this  decay  is  included,  the 
type  of  sublayer  affects  the  deflection  values  at  all  lev- 
els of  ESAL  (Table  6.14,  6.15,  and  6.16).  The  values  of 
ESAL  have  a  noticeable  effect  on  the  load  transfer  effi- 
ciency. The  load  transfer  efficiency  decreases  as  the 
number  of  ESAL  increases.  For  the  same  number  of  ESAL, 
the  thicker  the  slab,  and  the  stiffer  the  sublayers,  the 
lower  the  values  of  load  transfer  efficiency. 
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Numerical  results  were  obtained  assuming  several  pavement 
parameters  ,  to  evaluate  the  performance  of  the  pavement 
section  shown  in  Fig.  6.1  when  pumping  development  is 
included.  As  described  before,  four  different  thicknesses, 
three  different  traffic  levels,  and  three  different  values 
of  subgrade  modulus  were  used. 


First,  pumping  indices  were  determined  for  the  twelve  dif- 
ferent  cases.   The   results   are   shown  in  Tables  6.17  to 

6.20.  The  same  results  are  plotted  in   Fig.   6.4   to   6.6. 

3 

These       pumping    indices    are    given    in       in       /in       ,     as    in    the 

3 

AASHO  Road  Test,  and  in   dm   /dm   .   The   areas   with   no 

support  were  obtained  assuming  an  average  void  depth  of 
6.3  mm  (0.25  in.).  The  curves  were  drawn  by  graphically 
fitting  lines  through  the  pumping  index  values  obtained 
from  the  analysis.  As  shown  in  the  Figures,  the  pumping 
index  varies  seriously  with  the  number  of  ESAL,  the  thick- 
ness of  the  slab,  and  the  stiffness  of  the  sublayer. 
Since  the  state  of  cracking  of  the  concrete  slab  varies 
from  element  to  element,  the  stiffness  properties  of  the 
slab  vary  in  a  discrete  manner.  Therefore,  the  variations 
of  deflections,  amount  of  damage,  and  damaged  area  are 
also  discrete.  To  compensate  for  this,  regression  equa- 
tions are  obtained  to  fit  the  numerical  values. 
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Table    6.17       Pumping    index    for    different     levels     of 
ESAL    and    different    types    of    sublayers,     t-152    mm. 


Sublayer   type 

ESAL 

100000  . 

type    1                       type    2 

ty  pe    3 

7.560                        4.200 

3.000 

48.774                       27.097 

19  .355 

1000000  . 

39.760                      22.020 

15  .750 

256.516                    142.064 

101.613 

10000000  . 

211.660                    134.500 

83.390 

1365.546                    867.740 

537  .999 

sq  .  in. 

s  q  .  cm 

Table    6.18       Pumping    index    for    different     levels    of 
ESAL    and    different    types    of    sublayers,     t-203    mm. 


Sublaver 

type 

ESAL 
100000. 

type    1 

type    2 

type    3 

2.900 

1  .210 

0.620 

18.710 

7  .806 

4  .000 

1000000. 

15.180 

6.360 

3.260 

97  .935 

41  .032 

21  .032 

10000000. 

81  .580 

34.600 

18.260 

526.322 

223.225 

1 17  .806 

sq . in. 

sq  .  cm 
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Table    6.19    Pumping    index    for    different     levels    of 
ESAL    and    different    types    of    sublayers.     t=254    mm. 


Sublayer 

type 

ESAL 
100000  . 

type    1 

type    2 

type    3 

1  .180 

0.310 

0.060 

7  .613 

2  .000 

0.387 

1000000  . 

6.200 

1.670 

0.350 

40.000 

10.774 

2  .258 

10000000. 

32.500 

8.950 

2.000 

209.677 

57.742 

12.903 

sq  .  in. 

sq  .  cm 

Table    6.20       Pumping    index    for    different     levels    of 
ESAL    and    different    types    of.    sublayer.     t»305    mm. 


Sublayer 

type 

ESAL 
100000  . 

type    1 

type    2 

type    3 

0.410 

0.040 

0. 

2.645 

0.258 

0. 

1000000. 

2.170 

0.220 

0. 

14 .000 

1  .419 

0. 

10000000  . 

11  .570 

1.300 

0. 

74  .645 

8.387 

0. 

s q  .In. 

sq  .  cm 
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Figure  6.4  Pumping  indices  for  different  slab 
thicknesses  and  number  of  load  applications. 
Sublayer  type  1  . 
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75 
(4.84) 


Sublayer  Typ«  2 


Figure  6.5  Pumping  indices  for  different  slab 

thicknesses  and  number  of  load  applications. 

Sublayer  type  2. 
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10" 

log(ESAL) 


Figure  6.6  Pumping  indices  for  different  slab 
thicknesses  and  number  of  load  applications. 
Sublayer  type  3. 
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Deflections  at  node  57  were  determined  for  the  twelve 
cases.  From  these  values,  a  regression  equation,  in  terms 
of  ESAL,  thickness  (t  in.),  and  modulus  of  subgrade  reac- 
tion (k  pci)  was  obtained  as  follows: 


w5?  =  exp(al+a2  log(ESAL)  a3  t   +   a4  k  )      (6.5) 


w57   =  deflection  at  node  57  (mil) 


coef  f  icient 
al=  -0.6714 
a2  =   0.339 
a3=  -0.422 
a4=  -0.0022 


95%  confidence  interval 
-1.556         0.213 
0.216         0.463 
-0.467        -0.377 
-0.00283      -0.00160 


=  0.93 


The  numerical  results  are  given  in  Tables  6.21  to  6.24, 
and  graphical  representations  of  the  regression  equation 
are  given  in  Figures  6.7  to  6.  9.  These  Figures  show 
the  deflections  caused  by  an  80.1  kN  (18-K)  single  load. 
This  load  is  applied  after  the  given  number  of  ESAL 
repetitions  has   been  imposed  on  the  pavement. 
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Table    6.21       Deflections    produced    by    a    18-K    single 

axle     load,     for    different     levels     of    ESAL    and 

different     types     of     sublayers.     t=152    mm. 


Sublayer 

tyoe 

ESAL 

100000  . 

type     1 

type    2 

type    3 

203.900 

102  .300 

74.210 

5.179 

2  .598 

1  .885 

1000000  . 

215  .000 

134  .900 

H5  .400 

5  .461 

3  .426 

2  .931 

10000000  . 

261  .400 

168  .900 

129.100 

6  .640 

4  .290 

3  .279 

ml  Is 

ram 

Table    6.22       Deflections    produced    by    18-K    single 

axle    load,     for    different     levels     of    ESAL    and 

different    types    of    sublayer.     t=203    mm. 


Sublayer 

type 

ESAL 
100000  . 

type    1 

type    2 

type    3 

36.050 

22  .430 

16  .470 

0  .916 

0.570 

0  .418 

1000000  . 

114  .300 

71.010 

24  .390 

2.903 

1  .804 

0  .620 

10000000. 

124.600 

97.790 

71  .210 

3.165 

2  .484 

1  .809 

mils 

mm 
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Table    6.23       Deflections    produced    by    18-K    single    axle    load, 
for       different       levels     of    SSAL    and    different     types     of     sub-, 
laye  r .     t  =  254    mm . 


Sublayer 

type 

ESAL 
100000. 

type    1 

type    2 

type    3 

20  .620 

13.540 

10  .210 

0.524 

0.344 

0.259 

1000000  . 

36  .680 

17  .270 

11  .390 

0.932 

0.439 

0.289 

10000000  . 

38.190 

28.320 

14  .480 

0.970 

0.719 

0.368 

mi  Is 

mm 

Table  6.24  Deflections  produced  by  18-K  single  axle  load, 
for  different  levels  of  ESAL  and  different  types  of  sub- 
layer.    t*305    mm  . 


^ubla^pr   tvpe 

ESAL 
100000  . 

type    1 

type    2 

type    3 

14  .360 

9.960 

8.290 

0.365 

0.253 

0.211 

1000000  . 

18.360 

10.280 

8.360 

0.466 

0  .261 

0.212 

10000000  . 

29  .150 

12  .570 

8.430 

0.740 

0.319 

0.214 

ml  Is 

mm 
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Figure    6.7    Deflections     caused    by    a    80.1     kN       (18-K)       single 
axle    load,     for    different     levels    of    ESA.L    repetitions    and    for 
different    thicknesses.    Sublayer    type    1 
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Figure  6.8  Deflections  caused  by  a   80.1  KN  (18-K)  single 
axle  load,  for  different  levels  of  ESAL  repetitions  and  for 
different  thicknesses.  Sublayer  type  2. 
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Figure  6.9  Deflections  caused  by  a  80.1  KN   (18-K)  single 
axle  load,  for  different  levels  of  ESAL  repetitions  and  for 
different  thicknesses.  Sublayers  type  3. 
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A  better  fitting  of  the  regression-calculated  deflections 
with  deflections  calculated  by  means  of  the  structural 
analysis  could  have  been  obtained  if  individual  regression 
lines  were  obtained  for  particular  cases  of  thickness  or 
subgrade  modulus.  The  values  of  subgrade  modulus  for  the 
regression  equation  are  the  static  modulus  of  subgrade 
reaction.  A  conversion  factor  of  0.6  was  used  to  obtain 
these  static  values  given  the  average  value  of  resilient 
modulus  for  a  range  of  deflections  from  0.  to  40.  mils. 
These  values  were  as  follows: 


type  of  sublayer 

type  1  (soft) 

type  2  (medium) 

type  3  (stiff) 


static  k 

54.1  MPa/m  (200.  pci) 

108.2  MPa/m  (400.  pci) 

162.3  MPa/m  (600.  pci) 


The  amount  of  damage,  quantified  as  explained  in  Chapter  5 

2  2 

and  expressed  in   dm     (   in.    ),  was  determined  for  the 

twelve  cases.  The  numerical  results  obtained  from  the  ana- 
lyses (shown  in  Tables  6.25  to  6.28)  were  used  to  calcu- 
late a  regression  equation  for  amount  of  damage  in  terms 
of  ESAL,  thickness,  and  type  of  sublayer.  This  equation 
is  as  f ol lows : 
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D  =  exp((ATAN(al+a21og(ESAL)+a3  t+a4  k)5.5)    (6.6) 

2 
D=fatigue  damage  (   in.    ) 

t=slab  thickness  (in.) 

k=static  subgrade  modulus  (pci) 


coefficient 
al-  11.985 
a2=  4.718 
a3=  -3.183 
a4  =   -0.0214 


95%  confidence  interval 

-15.1093       39.0794 

1.1204        8.3190 

-4.9828       -1.3835 

-0.0394       -0.0034 


=  0.54 
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Table    6.25       Amount    of     damage     caused    at     different     levels 
of    ESAL    and    for    different     types    of     sublayer.     t=152    mm. 


SuW.it 

•er   cype 

ESAL 
100000  . 

type     1 

cype    2 

type    3 

3481  .000 

179  1  .000 

1531  .300 

224.559 

115.537 

98.784 

1000000  . 

3906.000 

2827.000 

1937  .000 

251  .976 

182.370 

124  .956 

10000000  . 

5083.000 

3780.000 

3010.000 

327  .904 

243.848 

194  .175 

sq  .  in. 

3  q  .dm 

Table    6.26       Amount    of    damage    caused    at    different    levels 
of    ESAL    and    for    different    types    of    sublayer.     t*203    mm. 


Sublayer 

type 

ESAL 
100000. 

type     1 

type    2 

cype    3 

13.400 

5.900 

2.810 

0.864 

0.381 

0.181 

1000000  . 

2589  .500 

1462.800 

89  .610 

167  .049 

94.365 

5  .781 

10000000  . 

3677  .000 

2523.000 

2129  .000 

237.203 

162.759 

137  .342 

sq  .  in. 

3q  .  dm 
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Table    6.27       Amount    of     damage     caused    at     different     levels 
of    ESAL    and    for     different     types     of     sublayer,     t-254    mm. 


Sublayer 

type 

ESAL 

type     1 

type    2 

type    3 

100000  . 

0.700 
0.045 

0  .400 
0  .026 

0.300 
0.0  19 

1000000  . 

73.200 
4  .722 

30  .000 
1  .935 

6.500 
0  .4  19 

10000000  . 

1  180.000 

76.122 

620.700 
40  .04  1 

53  .000 
3.4  19 

sq .in. 

sq  .  dm 

Table    6.28       Amount    of    damage    caused    at    different    levels 
of    ESAL    and    for    different    types    of    sublayer,     t-305    mm. 


Sublayer   type 

ESAL 
100000  . 

type 

1                     type    2 

type    3 

0. 

0. 

0. 

0. 

0. 

0. 

1000000  . 

0. 

0. 

0. 

0. 

0. 

0. 

10000000. 

0. 

0. 

0. 

0. 

0. 

0. 

9q  .  in. 

sq  .  dm 
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Figure    6.10    Amount    of    damage    caused    by    different     levels 
of    ESAL    repetitions    and    for    different    thicknesses. 

Sublayer    type    1 . 


154 


5000 
(323) 


3750 
(242) 


2500 
(162) 


a 

«-      1250 

°       (8  1) 


0 
(0) 


Sublayer  Type  2 

t=  1  50  «in/ 

"    / 

t  =  205  mil/ 

,       .     .    -*e<.  i 

1             i l_*— ..J-l- 

t»254  mm/ 

10' 


logCESAL 


Figure    6.11    Amount    of    damage    caused    by    different    levels 

of    ESAL    repetitions    and    for    different    thicknesses. 

Sublayer    type    2. 
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Figure  6.12  Amount  of  damage  caused  by  different  levels 

of  ESAL  repetitions  and  for  different  thicknesses. 

Sublayer  type  3. 
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The  areas  damaged  under  fatigue  in  the  concrete  slab  were 
determined  for  the  twelve  cases  (Tables  6.29  to  6.32). 
With  the  numerical  results  obtained  from  the  computer 
method  of  analysis,  a  regression  line  was  obtained  to  fit 
all  the  cases.  This  regression  equation  is  as  follows: 

DA  =  exp(ATAN(al+a2  log(ESAL)+a3  t  +a4  k)6.)   (6.7) 
whe  re : 


DA=damaged  area  (   in.    ) 
t=slab  thickness  (in.) 
k=static  subgrade  modulus  (pci) 


coeff  icient 
al=  39.006 
a2=   3.941 
a3=  -4  .387 
a4=  -0.036 


95%  confidence  interval 

15.928      63.271 

0.346       6.636 

-5.959      -2.815 

-0.521      -0.021 


R    =0.73 


Graphical  representation  of  this  equation  is  given  in  Fig- 
ures 6.13  through  6.15. 
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Table  6.29   Size  of  area  damaged  by  different  numbers  of 
ESAL  applications  and  for  different  types  of  sublayers. 

t" 1 5  2  mm. 


Sublayer   type 

ESAL 

100000 . 

type    1 

type    2 

type    3 

10400  .000 

5400.000 

4000  .000 

670.904 

348.354 

258  .040 

1000000. 

10600.000 

8200.000 

6400.000 

683.806 

528.982 

412.864 

10000000. 

10800.000 

8800.000 

7800.000 

696.708 

567  .688 

503.178 

Bq  .in. 

sq  .  dm 

Table  6.30   Size  of  area  damaged  by  different  numbers  of 
ESAL  applications  and  for  different  types  of  sublayers. 

t=»203  mm. 


Sublayer 

type 

ESAL 

100000  . 

type    1 

type    2 

type    3 

5000.000 

2600.000 

1400.000 

322.550 

167  .726 

90.314 

1000000  . 

9000.000 

6000.000 

3600.000 

580.590 

387  .060 

232.236 

10000000  . 

9600.000 

6800.000 

5800.000 

6 19  .296 

438.668 

374.158 

sq    In. 

sq.   da- 
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Table    6.31       Size    of    area    damaged    by    different     numbers     of 
ESAL    and    for    different     types     of     sublayers. 
t*254    mm. 


s^laver 

cvpe 

ESAL 

100000  . 

type     1 

type    2 

type    3 

1000.000 

200.000 

200  .000 

64  .5  10 

12  .902 

12  .902 

1000000  . 

3400  .000 

1000  .000 

400.000 

219  .334 

64  .510 

25  .804 

10000000. 

4800.000 

2600.000 

500.000 

309  .648 

167  .726 

32  .255 

sq    In. 

sq    dm 

Table    6.32       Size    of    area    damaged    by    different    numbers    of 
ESAL    applications    and    for    different    types     of    sublayers. 

t»305    mm. 


SuM  iver   type 

ESAL 

100000  . 

type 

1                      type    2 

type    3 

0. 

0. 

0. 

0. 

0. 

0. 

1000000  . 

0. 

0. 

0. 

0. 

0. 

0. 

10000000  . 

0. 

0. 

0. 

0. 

0. 

0. 

sq  .  in. 

sq  .  dm 
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Figure    6.13       Area    of    slab    damaged    under    flexural    fatigue; 
for    different     levels    of-ESAL    and    different     slab    thicknesses 
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Figure  6.14   Area  of  slab  damaged  under  flexural  fatigue; 
for  different  levels  of  ESAL  and  different  slab  thicknesses 

Sublayer  type  2. 
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Figure  6.15   Area  of  slab  damaged  under  flexural  fatigue; 
for  different  levels  of  ESAL  and  different  slab  thicknesses 

Sublayer  type  3. 
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6_.5_  Summary 
A  given  pavement  section  was  used  to  obatain  numerical 
results  of  pavement  response  using  the  method  of  analysis 
presented  in  Chapter  5.  These  results  were  used  to  study 
and  demostrate  the  sensitivity  of  the  method  of  analysis 
to  variations  in  sevaral  pavement  parameters.  Addition- 
ally, numerical  results  were  also  obtained  to  determine 
the  response  of  a  given  pavement  section  to  different 
traffic  levels.  Pavement  response  was  quantified  by 
means  of  amount  of  fatigue  damage  caused  in  the  concrete 
slab,  decay  in  load  transfer  efficiency,  and  maximum 
deflections.  The  parameters  which  cause  the  largest  varia- 
tions in  pavement  response  are:  slab  thickness,  number  of 
equivalent-single-axle-load  applications,  modulus  of 
subgrade  reaction,  and  flexural  behavior  of  concrete.  The 
numerical  values  are  given  in  Tables  6.1  to  6.10  for  the 
case  in  which  only  fatigue  damage  in  the  concrete  slab  is 
considered.  Tables  6.11  to  6.16  show  the  same  type  of 
results  but  for  the  case  in  which  decay  in  load  transfer 
efficiency  is  also  considered.  Higher  variations  of  the 
amount  of  damage  occur  when  decay  in  load  transfer  effi- 
ciency is  included.  The  number  of  load  repetitions  has  a 
noticeable  effect  on  the  load  transfer  efficiency.  Pave- 
ments with  thicker  slab  or  stiffer  sublayers  have  lower 
values  of  load  transfer  efficiency  after  a  given  number  of 
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load  repetitions.  Pavement  performance  was  determined 
for  different  traffic  levels  assuming  that  pumping  takes 
place.  In  this  case,  in  addition  to  considering  fatigue 
damage  in  the  slab  and  decay  in  the  load  transfer  effi- 
ciency, pumping  of  the  material  under  the  slab  was  also 
considered.  Results  of  pavement  performance  were 
presented  in  Tables  6.17  to  6.32.  Additionally,  to  com- 
pensate the  discrete  variation  of  these  results,  regres- 
sion equations  were  obtained'  for  values  of  deflections, 
amount  of  fatigue  damage,  and  damaged  area,  as  a  function 
of  slab  thickness,  number  of  ESAL,  and  type  of  sublayer. 


The  computer  scheme  presented  in  Chapter  5  was  used  to 
obtain  the  numerical  results  presented  in  this  Chapter. 
The  pavement  geometries  as  well  as  the  pavement  properties 
presented  at  the  beginning  of  this  Chapter  were  used  as 
input  data.  Loads  were  specified  in  terms  of  number  of 
18-K  ESAL  applications.  Additionally,  a  specified  number 
of  iterations  of  15  to  achieve  convergence  was  used.  A 
complete  listing  of  the  computer  scheme  is  given  in  the 
Appendix.  Additionally,  in  this  Appendix  an  example  with 
the  input  and  output  data  is  also  provided. 


164 


CHAPTER  7   SUMMARY   CONCLUSIONS  AND 
RECOMMENDATIONS  FOR  FURTHER  STUDY: 


Since  the  1940's,  pumping  has  been  considered  as  a  major 
contributor  to  failure  in  concrete  pavements.  It  has  been 
reported  [2,3,4,51]  that  most  of  the  failures  in  concrete 
pavements  have  been  preceded  by  pumping  of  the  material 
underneath  the  concrete  slab.  Traditional  design  pro- 
cedures, based  on  a  stress-fatigue  criteria,  have  been 
unsuccessful  in  providing  pavements  which  accomplish  their 
intended  service  lives.  Premature  failures  due  to  exces- 
sive pumping,  erosion  of  sublayers  and  joint  faulting  have 
occurred.  Therefore,  in  addition  to  the  stress-  fatigue 
criterion,  an  erosion  criterion  must  be  considered. 

The  conventional  methods  for  structural  analysis  of  con- 
crete pavements  do  not  consider  the  variation  in  the  pave- 
ment properties  caused  by  repetitive  application  of  loads. 
In  these  methods,  it  is  assumed  that  the  geometric  and 
mechanical  properties  retain  the  same  values  through  the 
life  of  the  pavement. 

The  major  limitation  in  the  analysis  method  developed  by 
Westergaard,  or  any  method  based  on  his  equations,  is  that 
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it  is  applicable  only  if  the  analyzed  structure  is  con- 
tinuous. A  rigid  pavement  is  composed  of  several  elements 
with  different  mechanical  behaviors.  In  general,  the  pave- 
ment can  be  divided  into  a  concrete  slab,  sublayers  and 
load  transfer  devices.  Each  one  of  these  components 
behaves  in  a  completely  different  manner,  therefore,  it  is 
impossible  to  use  a  general  differential  equation  to 
represent  the  behavior  of  the  whole  structure.  Addition- 
ally, the  existence  of  joints,  cracks,  and  gaps  between 
slab  and  sublayers  are  very  common,  which  makes  the 
analysis  more  complex  and  the  solution  in  a  closed  form 
impossible.  Therefore,  these  methods  are  used  only  to 
check  results  obtained  with  other  methods  using  simple, 
i  dea lized  cases. 


One  of  the  main  causes  of  pavement  failures  is  the 
deterioration  and  gradual  weakening  of  the  pavement  com- 
ponents. Repetitive  loading  produces  micr ocracking  in  the 
slab  as  well  as  erosion  and  voids  in  the  sublayers.  This 
reduces  the  resistance  of  the  pavement  to  the  loads 
applied  by  traffic.  The  response  of  the  pavement  to  a 
specific  load  is  different  throughout  its  entire  service 
life.  Cumulative  damage  reduces  the  general  stiffness  of 
the  structure.  Therefore,  in  the  analysis  of  rigid  pave- 
ments, especially  if  erosion  and  pumping  are  expected,  the 
pavement  components  should  be  considered  as  having  a 
stress   dependent   behavior.    The  response  of  the  slab  as 


166 

well  as  that  of  the  sublayers  should  be  dependent  upon  the 
previous  stress  history  applied  on  the  pavement.  Pavements 
properties  should  be  revised  every  few  years  to  take  into 
account  the  reduction  of  stiffness  caused  by  exposure  to 
traffic  loads.  In  this  way,  a  more  exact  calculation  of 
stresses  and  deflections  in  the  pavement  at  different 
6tages  of  its  life  i6  accomplished.  Therefore,  a  more  pre- 
cise estimation  of  the  number  of  allowable  load  repeti- 
tions (and  thereby  of  service  life)  can  be  achieved.  Then, 
the  designer  has  more  and  better  information  to  carry  out 
an  economical  analysis. 

The  present  methods  only  give  results  for  specific  proper- 
ties that  the  pavement  has  at  the  moment  considered  in  the 
analysis.  They  cannot  provide  the  stresses  and  deflections 
that  occur  after  the  pavement  has  been  in  service  a  cer- 
tain time.  The  state  of  the  pavement  after  several  years 
in  service  cannot  be  determined,  and  thus,  the  maintenance 
cost  cannot  be  accurately  estimated. 


The  mechanical  behavior  of  concrete  has  been  described 
[45]  as  nonlinear  at  low  stresses,  and  with  expansion  near 
failure.  Nonlinearity  is  caused  by  mi crocra eking  due  to 
segregation,  shrinkage,  or  temperature  changes.  The  ini- 
tial cracks  that  exist  in  concrete  develop  and  propagate 
as  stresses  are  applied.  Near  failure,  cracks  propagate 
and  bridge  together  causing  disruption  in  the  concrete 
mass.   When   reinforcing   steel   is   provided,   the  cracks 
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distribute  in  areas  larger  than  those  under  high  stresses. 
Steel  keeps  the  concrete  areas  together  up  to  stresses 
under  which  disruption  would  otherwise  occur.  The  avail- 
able analysis  methods  assume  a  linear  elastic  behavior  of 
concrete  at  all  stress  levels.  Presence  of  reinforcing 
steel  is  not  considered.  Therefore,  the  strains 
corresponding  to  certain  stresses,  especially  high 
stresses,  are  different  from  the  actual  corresponding 
strains. 

In  laboratory  tests  [45]  the  values  of  strength  and 
modulus  of  elasticity  obtained  for  concrete  depend  upon 
the  rate  with  which  the  load  is  applied.  The  values 
obtained  in  an  impact  test  are  higher  than  those  from  a 
static  slow-rate  test.  In  the  pavement  slab,  loads  are 
more  similar  to  impact  loads  than  to  static  loads.  There- 
fore the  use  of  static  strength  and  static  modulus  of 
elasticity  is  unrealistic. 

In  the  sublayers,  the  nonlinear  behavior  of  the  layers  is 
even  more  pronounced.  Deformations  in  the  sublayers  occur 
as  a  result  of  relative  movements  of  aggregates  and 
grains.  Consolidation  can  occur.  The  stiffness  properties 
of  the  sublayers  may  be  even  more  sensitive  to  the  rate  of 
load  application.  Therefore,  the  assumption  of  static 
values  with  linear  behavior  may  be  far  from  reality. 


Loads  in  pavements  are  applied   in   a   repetitive   manner. 
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This  produces  fatigue  in  the  slab,  the  sublayers,  and  the 
joint  devices.  Fatigue  produces  a  decrease  in  the  pavement 
stiffnesses.  Therefore,  in  the  available  methods  of 
analysis,  the  assumption  of  constant  mechanical  properties 
for  the  slab,  the  sublayers,  and  the  joint  devices  is  not 
realistic. 


From  the  point  of  view  of  the  mechanical  behavior,  an 
ideal  method  of  analysis  should  reproduce  the  actual  per- 
formance of  the  pavement  structure.  The  stress-strain 
relations  of  the  individual  components  should  be  accu- 
rately represented  in  the  analytical  model.  Concrete  in 
the  pavement  slab  can  be  represented  by  a  nonlinear  plas- 
tic model.  The  effect  caused  by  the  presence  of  micro- 
cracks  should  be  reflected  in  the  stress-strain  relation- 
ships. Similarly,  the  prescence  of  reinforcing  steel  must 
be  taken  into  account.  In  the  sublayers,  the  sublayers 
should  be  represented  by  a  nonlinear  plastic  model. 

Traffic  loads  are  repetitive  in  nature.  This  causes 
fatigue  of  all  the  stressed  pavement  elements.  The  con- 
crete properties  should  be  dependent  upon  the  loading  his- 
tory. Decreases  in  strength  and  modulus  of  elasticity 
should  be  considered,  so  should  the  formation  of  cracks 
and  break-ups  in  the  slab.  In  the  sublayers,  repetitive 
loading,  and  excess  water  may  create  erosion  and  pumping 
leading  to  the  formation  of  voids  underneath  the  slab. 
This  loss  of  support  should  be  considered  in  the  analysis. 
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The  values  considered  for  the  mechanical  properties  should 
be  the  values  obtained  with  loads  similar  in  nature  to  the 
actual  traffic  loads.  Impact  values  for  the  modulus  of 
elasticity  of  concrete  and  for  the  stiffness  of  the  sub- 
layers should  be  used.  These  values  more  accurately 
represent  the  performance  of  the  pavement. 

Load  transfer  devices  at  joints  are  also  subjected  to 
repetitive  loading.  Therefore,  the  method  of  analysis 
should  consider  the  decay  in  load  transfer  efficiency. 
This  decay  can  vary  from  a  quantitative  variation  of  the 
stiffness  values  to  a  qualitative  variation  of  the  actual 
performance.  The  analysis  method  should  follow  this  kind 
of  performance. 

A  non-linear  finite  element  method  for  analysis  of  rigid 
pavements  subjected  to  fatigue  damage  caused  by  traffic 
has  been  presented.  The  method  is  implemented  in  an 
iterative  computer  scheme  called  PMARP  in  which  the  pave- 
ment properties  are  modified  in  accordance  with  the 
fatigue  damage  produced  by  repetitive  loading. 

In  a  rigid  pavement,  most  of  the  structural  support  is 
provided  by  the  concrete  slab.  This  slab,  then,  under  the 
action  of  traffic  is  subjected  to  many  load  repetitions 
through  the  service  life  of  the  pavement.  The  load 
repetitions  produce  gradual  damage  in  all  the  pavement 
components    until,    eventually,   cracking   and   faulting 
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occurs,  producing  failure  of  the  structure.  In  jointed 
concrete  pavements,  the  devices  needed  to  transfer  the 
load  across  the  joints,  especially  transverse  joints,  also 
lose  their  efficiency  along  the  life  of  the  pavement. 
Fatigue  damage  imposed  to  the  concrete  slab  reduces  its 
strength  and  modulus  of  elasticity  because  of  the 
occurrence  of  microcracks.  Similarly,  fatigue  damage 
caused  by  repetitive  loading  reduces  the  load  transfer 
efficiency  of  the  dowel  bars  by  increasing  the  degree  of 
looseness  of  the  bars.  These  decays  in  the  properties  of 
the  pavement  components  modify  the  response  of  the  pave- 
ment to  a  given  load.  Stresses  and  deflections  calculated 
assuming  the  initial 

stiffness  values  are  different  from  those  obtained  assum- 
ing the  values  reduced  in  accordance  with  fatigue  damage. 
Assuming  fatigue  damage,  results  closer  to  actual  field 
conditions  are  obtained. 


In  the  computer  implementation  of  the  method  of  analysis 
presented  herein,  the  modifications  of  the  properties  of 
the  concrete  slab  are  based  on  the  flexural  behavior  of  a 
particular  concrete,  specified  as  part  of  the  input  data. 
Similarly,  the  modification  of  the  values  of  the  resilient 
modulus  of  subgrade  reaction  is  made  on  the  basis  of 
values  of  k  specified  in  the  input  data.  For  the 
modification  of  the  load  transfer  efficiency,  however, 
only  the  initial  value  can  be  specified  in  the  input  data. 
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The  fatigue  behavior  is  controlled  by  an  expression  built 
into  the  computer  scheme.  Nevertheless,  this  expression 
can  be  modified  by  changing  the  values  of  the  coefficients 
used  in  the  computer  implementation.  In  any  event,  accu- 
racy of  the  absolute  values  obtained  with  the  method  pro- 
posed herein,  depends  on  the  accuracy  of  the  information 
provided  in  the  input;  especially  the  information  concern- 
ing the  fatigue  behavior  of  the  different  pavement  com- 
ponents. The  relative  results  ,  however,  are  more  accu- 
rate and  especially  useful  in  comparative  studies. 

Additionally,  a  pumping  model  has  been  developed  based  on 
data  of  the  AASHO  Road  Test.  This  pumping  model  has  been 
implemented  in  the  method  of  analysis  presented.  With  this 
pumping  model,  implemented  in  a  subroutine  called  "pump", 
the  pumping  index,  as  well  as  the  support  conditions 
occurring  in  a  pavement  section  can  be  determined  for  any 
number  of  load  applications. 

The  method  of  analysis  presented  herein  has  the  following 
features: 

1.  Non-linear  finite  element  method  for   analysis   of 
cracked  and  jointed  concrete  pavements. 

2.  Finite  element  discretization: 


2.1  Slab;   orthotropic   cracked   and   uncracked   thin 
plates,  with  or  without  reinforcement 
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2.2  Sublayers:   Winkler   sublayers   with    resilient 
subgrade  modulus 

2.3  Joints:  load  transfer  devices  as  beam  elements 
3.  Long  term  effect  of  load  applications. 


3.1  Fatigue  damage  and  cracking  can  occur  in  concrete 


s  lab 


3.2  Fatigue  decay   in   load   transfer   efficiency   of 
dowel  bars  can  occur 

3.3  Pumping  development  and  loss  of  support  can  occur 
4.  Results  from  analysis: 

4.1  Deflections  at  all  nodal  points 

4.2  Final  values  of  subgrade  moduli  at  all  nodes 

4.3  State  of  stresses  at  all  nodal  points 

4.4  Amount  of  damage  caused  to  concrete  slab 


4.5  Damaged  area  in  concrete  slab 


4.6  Total  decay  in  load  transfer  efficiency 


4.7  Potential  of  pumping 


4.8  Area  with  no  support  caused  by  pumping 
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The  computer  implementation  of  the  proposed  method  for 
analysis  of  rigid  pavements  considering  fatigue  is  not  a 
very  critical  task.  A  computer  scheme  of  the  method  of 
analysis  presented  in  Chapter  5  as  well  as  an  example  are 
given  in  the  Appendix.  It  is  more  difficult  to  obtain  the 
appropriate  experimental  information  to  use  the  analytical 
model.  Experimental  and  field  results  are  required  to 
develop  the  following  tasks:  (1)  accurate  representation 
of  behavior  of  plain  and  reinforced  concrete  in  flexural 
fatigue;  (2)  Behavior  of  load  transfer  devices  under 
repetitive  loading;  and  (3)  accurate  pumping  model  which 
considers  the  following  factors:  structural,  climatic, 
traffic,  and  erosion  resistance. 


For  a  typical  pavement  section,  the  parameters  which  pro- 
vide the  largest  effects  on  pavement  response  are:  slab 
thickness;  number  of  ESAL  repetitions,  modulus  of  subgrade 
reaction,  and  flexural  behavior  of  concrete.  The  numerical 
results  obtained  by  means  of  the  method  of  analysis 
presented  here  vary  in  a  discrete  manner  due  to  the  fact 
that  cracking  and  void  formation  are  allocated  in  such  a 
manner.  This  discrete  variation  can  be  solved  in  two  ways: 
by  refining  the  finite  element  mesh  so  that  the  discrete 
changes  are  small,  or  by  using  regression  equations  to 
represent  the  numerical  values  obtained  from  the  analysis. 
Refining  the  finite  element  mesh  would  require  large  com- 
puter memory  which  in  some  cases  could   not   be   provided. 
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The  results  obtained  in  Chapter  6  are  represented  by 
regression  equations.  These  regression  equations  are 
plotted  to  provide  an  easy  way  to  determine  the  pavement 
performance  given  the  pavement  features  and  the  traffic 
conditions.  The  results  provided  in  Chapter  6  are 
obtained  for  the  assumptions  explained  in  Chapter  5.  Pave- 
ment performance  is  very  sensitive  to  the  support  condi- 
tions which  are  defined  by  the  degree  of  pumping.  The 
pumping  model  provided  in  Chapter  5  considers  only  the 
structural  features  of  the  pavement  and  the  volume  of 
traffic.  Any  information  concerning  other  factors  ,  such 
as  amount  of  rainfall  ,  erosion  susceptibility  of  the  sub- 
base,  or  joint  sealing  ,  should  be  included.  This  informa- 
tion can  be  included  by  using  factors  which  modify  the 
potential  volume  of  pumped  material  to  obtain  the  actual 
volume  of  pumped    material. 

Further  research  to  develop  an  appropriate  method  of 
analysis  of  rigid  pavements  is  required  in  the  following 
a  reas  : 

1  Development  of  an  appropriate  nondestructive  method 
to  obtain  the  configuration  of  voids  under  concrete  pave- 
ments, so  that  a  pumping  model  can  be  verified  by  means  of 
actual  field  data. 


2  Development  of  accurate  methods  to   predict   amount 
of  pumping  based  on:  material  erodability,  drainage  condi 
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tions,  amount  of  rainfall,  structural  conditions  of   pave- 
ment, and  volume  of  traffic. 

3  Development  of  procedures  to  relate   pumping   index 
with  actual  loss  of  support  of  concrete  slab. 

4  Development  of  a  more   accurate   representation   of 
behavior  of  reinforced  concrete  under  flexural  stresses. 


5  Development  of  a  more  accurate  model   to   represent 
the  behavior  of  dowel  bars  under  repetitive  loading. 
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INPUT  GUIDE  AND  EXAMPLE 
Input  Guide 

The  "PMARP"  iterative,  nonlinear  finite  element  program 
can  be  used  to  analyze  reinf or ced/unreinf or ced  concrete 
pavements  with  joints.  Different  amounts  of  reinforcement 
may  be  used  in  the  longitudinal  and  transverse  directions. 
Joint  types  may  be:  doweled,  aggregate  interlock,  or 
keyed.  Only  one  layer,  constituted  by  the  concrete  slab, 
can  be  handled.  The  stiffness  of  the  sublayers  can  be 
represented  by  a  composite  resilient  impact  modulus. 
Loads  may  consist  of  one  application  with  any  configura- 
tion, or  of  several  repetitions  with  constant  configura- 
tion. Stresses  and  deflections  at  all  nodal  points,  after 
any  number  of  ESAL  load  repetitions,  can  be  obtained. 
Additionally,  fatigue  damage  in  the  concrete  slab,  and  in 
the  doweled  load  transfer  devices  can  be  quantified. 
Pumping  potential  can  also  be  determined  as  well  as  lose 
of  support  due  to  pumping. 

,  The  pavement  section  may  consist  of  1,  2,  3,  4,  or  6 
slabs.  with  nodes  numbered  from  left  to  right  and  from 
bottom  to  top  (Fig.  A.l). 
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